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Abstract 

One  of  the  biggest  challenges  for  Brazilian  military  logisticians  is  to  support 
effective  jungle  warfare  for  both  real  and  training  operations  carried  out  by  their  combat 
forces  in  the  Amazonian  region.  The  jungle's  heat,  humidity,  and  dense  vegetation  put 
significant  demands  on  the  supply  chain.  Further,  because  of  the  difficulties  of  land  or  air 
transportation,  water  transport  is  the  most  reasonable  transportation  option  to  sustain 
these  deployed  forces.  Planners  must  select  from  among  the  available  watercourses  those 
whose  surroundings  provide  stealthy  navigation  to  the  combat  force  location  where  the 
requested  supplies  can  be  safely  unloaded.  We  seek  a  method  of  determining  a  path 
through  a  river  network  that  blends  short  transit  times  with  maximal  shade  coverage  from 
forest  growth  along  the  riverbanks.  We  combine  an  astronomical  algorithm  for 
computing  shadow  coverage  with  Dijkstra's  shortest  path  algorithm  to  determine  the  start 
time  and  routing  information  necessary  for  a  supply  boat  to  travel  from  a  depot  to  a 
resupply  point  that  minimizes  weighted  risk,  which  is  defined  as  the  product  of  shade 
coverage  and  arc  transit  time. 
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STEALTHY  RIVER  NAVIGATION  IN  JUNGLE  COMBAT  CONDITIONS 


I.  Introduction 


Background 

One  of  the  biggest  challenges  for  Brazilian  military  logisticians  is  to  support 
effective  jungle  warfare  for  both  real  and  training  operations  carried  out  by  their 
combined  forces  (Army,  Navy,  and  Air  Force),  in  the  Amazonian  region  (figures  1  to  4). 
To  support  anti-guerilla  actions,  whether  riverine  or  air  interdiction  operations,  against 
the  drug  traffickers  operating  along  the  Brazil-Colombian  border  or  to  repel  any 
aggressive  international  invading  forces,  in  a  defensive  posture,  the  military  supply 
system  needs  to  adapt  to  the  nature  of  jungle  warfare  (Spencer,  2004). 


Figure  1.  Brazilian  Army  commando 
troops  during  jungle  operational  training  in 
Amazon  rainforest,  simulating  a  real 
tactical  infiltration  mission. 


Figure  2.  Embraer  Super  Tucano  Brazilian 
aircraft,  designed  and  employed  in 
operations  carried  out  in  jungle  warfare 
scenarios. 
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That  means  that  resupply  missions  cannot  disturb  the  unconventional,  low- 
intensity,  guerrilla-style  operations  taking  place  or  expose  the  location  of  hidden  military 
camps  (camouflaged  depots  and  airfields)  during  load/unload  operations.  The  impact  of 
both  the  challenges  of  a  harsh  environment  and  the  hostile  enemy  actions  against  troops 
and  logistics  resources  should  be  considered.  In  order  words,  the  peculiarities  of 
nonlinearity,  reduced  frontages,  and  short  fields  of  observation  found  in  the  jungle 
combat  zones  need  to  be  carefully  weighed  by  the  military  logisticians  while  they  are 
setting  up  the  logistical  capabilities  that  best  fit  the  operations  that  will  be  performed 
(DA,  1993:3-26). 

The  North  Vietnamese  constantly  repaired  and  extended  these  routes,  in  spite  of 
American  bombing.  Over  the  elaborate  trail  and  road  network,  enemy  troops,  fuel, 
and  munitions  flowed  southward.  Carried  most  of  the  way  in  a  series  of  short 
hauls,  with  repeated  changes  of  vehicles,  each  truck,  or  group  of  them,  continually 
shuttled  different  loads  over  the  same  short  stretch  of  road,  almost  always 
traveling  by  night.  At  various  key  points,  troops  could  rest  in  hidden  camps  and 
supplies  could  be  stored  in  carefully  camouflaged  depots.  (Cosmas  and  Murray, 
1986:279) 


Figure  3.  Brazilian  amphibious  com¬ 
mandos  during  riverine  warfare  training 
program. 


Figure  4.  Brazilian  Air  Force’s  Infantry  - 
Special  Force  troop  (dagger)  in  guerrilla 
tactics  training. 
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The  obstacles  and  opportunities  inherent  to  jungle  warfare  environment, 
extensively  evidenced  in  both  the  China-Burrna-India  (CBI)  campaigns  and  Pacific 
Theater  of  Operations,  during  the  Second  World  War,  as  well  as  the  Vietnam  War,  are 
quite  different  from  those  encountered  in  other  types  of  warfare.  As  a  result,  all  aspects 
of  jungle  warfare  are  highly  specialized  -  as,  indeed,  are  those  of  the  desert,  arctic, 
mountain  and  urban  varieties.  Inside  the  forests,  the  phases  of  war,  tactics,  training, 
logistical  support  and  administration  have  to  be  modified  because  of  trees,  which 
provide  multiple  constraints  to  military  operations  (Anderson,  1993;  Cross,  2007:11). 
For  example,  in  the  tropical  environment  of  canopy  forests,  swamps,  or  densely  forested 
mountains,  tanks,  and  heavy  artillery  are  of  little  use.  Similarly,  trucks  (figures  5  and  6), 
and  large  aircraft  are  inefficient  in  providing  regular  logistical  support  because  the 
jungle  limits  the  movement  of  personnel  and  material  among  the  sites  (Scott  and  others, 
2000:91-92). 


Figure  5.  “In  morass  of  mud,  one  jeep  la¬ 
bors  to  pull  out  another  mired  up  to  the 
axles.  The  road  passes  through  Hukawng- 
Mogaung  Valley  (Burma),  where  the  jungle 
is  thickest”.  (LIFE,  1944) 


Figure  6.  Endless  mudded  conditions  faced 
by  the  E  Battery,  82nd  Artillery,  during  the 
monsoon  seasons,  in  the  Vietnam  War. 
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The  dense  vegetation  and  general  lack  of  infrastructure,  along  with  reduced 
visibility  and  engagement  ranges,  also  make  it  extremely  difficult  to  support  deployed 
military  units,  or  any  other  kind  of  forces  on  a  large  scale  in  the  jungle  (Asprey, 
2002:838;  Thompson,  J.,  1994).  The  extract  below  sums  up  why  logistics  in  a  jungle 
environment  is  challenging  work. 

Many  of  the  techniques  and  assumptions  which  were  accepted  as  valid  in 
conventional  warfare  did  not  apply  in  the  harsh,  primitive,  jungle  environment 
and  the  isolated  support  enclaves.  Even  so,  Vietnam  is  a  story  of  remarkable 
logistics  achievement.  At  no  time  was  logistic  support  a  constraint  on  a  major 
tactical  operation.  This  record  was  made  despite  the  conditions  which  imposed  a 
fantastic  strain  on  logistics  operations  and  which  offered  an  enormous  challenge 
to  all  logisticians.  (Heiser,  1974:4) 

Key  Points: 

•  Monsoon  winds  alternately  encourage  and  discourage  most  military  operations 
in  South  and  Southeast  Asia. 

•  Geographical  circumstances  affect  supply,  maintenance,  transportation, 
medical,  and  other  logistical  requirements  at  least  as  much  as  combat 
operations. 

•  Logistical  problems  multiply  and  intensify  in  direct  proportion  to  the  distance 
between  support  bases  and  supported  forces. 

•  Jungle-covered  mountains  reduce  the  advantages  of  airmobile  forces  in  open 
terrain. 

•  Parachute  delivery  systems  and  helicopters  can  sustain  small,  isolated  units  in 
jungles,  but  large  formations  need  main  supply  routes  with  much  greater 
capacities.  (Collins,  1997/1998:126) 

The  lack  of  an  extensive  all-weather  transportation  network  in  many  jungle  areas 
makes  the  mission  of  support  units  more  difficult  (Scott  and  others,  2000:133-136). 
Transportation  difficulties  may  dictate  that  maneuver  units  be  resupplied  by  air,  pack 
animals,  or  human  portage  (figures  7  to  9).  For  example,  the  rapid  movement  of  military 
personnel  and  supplies  in  jungle  regions  is  severely  limited  by  both  the  dense  foliage  and 
the  uneven  terrain,  with  even  helicopters  being  only  of  limited  use  (Heiser,  1974; 
O'Sullivan,  1983:67;  Vongsavanh,  1981:4-5).  Suitable  natural  gaps  in  the  canopy  are 
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scarce  and  it  will  almost  certainly  be  necessary  to  clear  a  landing  site  by  hand  if  aircraft 
are  to  land. 


Figure  7.  Filipino  volunteers  carry  sup-  Figure  8.  A  long  line  of  porters  carry  supplies 
plies  into  the  mountains  to  reach  1st  along  the  Ho  Chi  Minh  Trail  for  the  North  Vi- 
Cavalry  Division  troops  in  the  Prov-  etnamese,  during  the  Vietnam  War.  (Correll, 
ince  of  Leyte  (Philippine  Campaign),  2005:62) 
during  the  Second  World  War.  (An¬ 
derson,  1993:24) 


Figure  9.  Native  carriers  bringing  supplies  through  the  jungles  into  the  hills,  in  the 
Guadalcanal  Island,  during  the  American  Pacific  Offensive  in  1942. 
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There  is  an  extensive  literature  documenting  the  burdens  faced  by  military 

commanders  regarding  how  to  conduct  and  support  their  troops  into  jungle  areas.  Below 

are  the  two  passages  which  briefly  summarize  these  concerns. 

Jungle  imposed  three  major  limitations  upon  the  infantry  that  fought  there  in 
1941-1945.  The  first  -  lack  of  visibility  and  fields  of  fire  -  necessitated  drastic 
revisions  of  small  unit  tactics.  The  second  -  remoteness,  and  obstruction  to 
transport  -  not  only  influenced  the  way  troops  fought  in  battle,  in  ambush  and  on 
patrol,  but  imposed  sometimes  insurmountable  problems  on  supply  and  the 
movement  of  heavy  equipment.  Getting  the  simplest  things  to  the  soldier  could  be 
extremely  difficult:  Gen  Slim  cited  specific  problems  with  milk,  rotting  cartons 
and  rusting  tins.  (...)  The  third  governing  condition  was  the  tropical  climate  itself; 
humidity  combined  with  insect  life  and  stagnant  water  to  produce  a  bewildering 
variety  of  diseases,  many  of  them  unknown  in  Europe.  Whilst  cholera  and  typhus 
were  killers,  it  was  debilitating  fevers,  and  above  all  malaria,  that  cut  a  swath 
through  the  ranks  of  those  who  fought  in  the  jungle.  (Bull,  2007:6) 

The  terrain  was  a  commander's  nightmare  because  it  fragmented  the  deployment 
of  large  formations.  On  the  north  shore  a  tangled  morass  of  large  mangrove 
swamps  slowed  overland  movement.  Monsoon  rains  of  eight  or  ten  inches  a  day 
turned  torpid  streams  into  impassable  rivers.  There  were  no  roads  or  railways,  and 
supply  lines  were  often  native  tracks,  usually  dirt  trail  a  yard  or  so  wide  tramped 
out  over  the  centuries  through  the  jungle  growth.  Downpours  quickly  dissolved 
such  footpaths  into  calf-deep  mud  that  reduced  soldiers  to  exhausted  automatons 
stumbling  over  the  glue-like  ground.  (Drea,  2008:3) 

The  dropping  of  supplies  by  parachute  is  also  not  very  effective  in  the  jungle. 
Supplies  can  easily  fall  into  the  hands  of  the  enemy  or  land  in  a  zone  attacked  by  fire. 
Supplies  can  land  among  the  trees,  steep  hills,  or  sink  in  the  streams  (figures  10  and  1 1). 
This  challenge  was  faced  by  the  American  troops  during  the  Papua  Campaign,  in  the 
Second  World  War,  between  23rd  July  1942  and  23rd  January  1943,  as  shown  in  the 
citation  below. 

The  efficiency  of  the  air  arm  in  direct  support  of  ground  troops  is  strictly 
curtailed.  The  complete  leaf  canopy  prevents  pilots  from  seeing  troops  on  the 
ground;  the  troops  are  often  unable  to  catch 

more  than  fleeting  glimpses  of  the  planes.  Pilots  cannot  se  panels  laid  out  on  the 
ground;  it  is  often  impossible  for  the]  to  see  panels  or  strips  displayed  in  the  tops 


6 


of  trees.  Colored  smoke  pots  placed  on  the  ground  have  been  used  to  indicate  the 
position  of  ground  troops  but  the  rapid  diffusion  of  the,  smoke  renders  them 
impractical.  Smoke  pots  or  smoke  grenades  to  which  lines  are  attached  may  be 
thrown  into  tree:  There  is  now  available  a  “tree  top”  smoke  grenade  which 
believed  to  be  satisfactory.  (DN,  1989:6) 

The  1st  October  plan  was  marked  by  the  innovation  which  would  characterize 
Mac  Arthur's  leadership  throughout  the  Pacific  War:  resupply  by  air.  Once  units 
entered  the  jungled  mountains,  resupply  became  a  major  problem.  The  Australian 
practice  of  relying  on  the  strong  backs  of  New  Guineans  did  not  solve  the 
problem,  since  the  bearers  usually  deserted  when  they  suspected  enemy  presence. 
The  Allies  settled  on  the  airdrop.  Expanding  its  range  as  fast  as  new  airfields 
could  be  constructed,  the  Fifth  Air  Force  proved  invaluable  in  overcoming  the 
obstacles  of  sea  distance  and  rugged  terrain.  Crates  of  food  and  supplies  were 
pushed  out  the  hatches  of  low-flying  C-7  over  breaks  in  the  jungle  ceiling. 
Though  not  perfect  -  hungry,  diseased  troops  sometimes  saw  crates  of  food, 
medicine,  and  ammunition  fall  down  mountainsides  just  out  of  reach  -  the 
airdrops  continued  and  improved  as  aircrews  gained  experience.  (Anderson, 
1993:8) 


Figure  10.  “Parachute  resupply  in  the 
Malayan  jungle  -  unlike  many  others, 
this  one  did  not  get  caught  up  in  the 
trees”.  (Cross,  2007:147) 
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Figure  1 1 .  Most  of  successful  air  resupply 
occurred  in  open  areas  instead  of  inside 
the  dense  wood.  (Thompson,  J., 1991:98) 


In  the  case  of  land  transport,  wheeled  vehicles  are  normally  restricted  to  roads  and 
wider  trails,  and  even  these  may  prove  impassable  during  heavy  rains  (Correll, 
1987/1988:65;  Schweitzer  and  Armstrong,  1966).  Sometimes,  goods  must  be  transported 
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by  cross  loading  from  wheeled  to  tracked  vehicles.  For  instance,  large  wheeled  vehicles 
move  the  supplies  as  far  forward  as  possible,  where  they  are  transloaded  to  tracked 
vehicles  that  move  them  cross-country.  In  rugged  terrain,  the  supplies  may  have  to  be 
further  transloaded  to  pack  animals  or  human  supply  bearers  (figures  12  and  13)  (Astor, 
2004;  Moreman,  2002:126). 


Figure  12.  “Muletrain”:  A  Chin- 
dit  column  moving  through  Bur¬ 
mese  Jungle. 


Figure  13.  Mars  Task  Force  mule  skinners  (2d  Bn., 
475th  Inf.  Regt.)  lead  mules  through  the  swift  river 
that  impeded  their  progress  to  Bhamo,  Burma,  17 
November  1944. 


During  the  Burma  campaign  from  December  1941  until  August  1945,  the  longest 
campaign  in  which  the  British  Army  fought  during  the  Second  World  War,  mules 
were  used  extensively,  exhaustively  and  often  hazardously.  Under  the  command 
of  General  Sir  William  Slim,  the  Fourteenth  Army  comprised  men  of  a  number  of 
diverse  nationalities  -  including  Burmese,  Chinese,  Gurkhas,  as  well  as  troops 
from  East  and  West  Africa.  In  March  1942  Rangoon  was  lost  and  the  Burma  army 
began  its  600-mile  retreat  to  the  Indian  frontier.  There  was  a  shortage  of  mules, 
particularly  the  larger  type  of  equipment-carrying  animals,  which  made  the 
decision  to  leave  the  roads  and  move  through  the  jungle  and  mountain  paths  hard 
to  implement.  Rail  and  road  communications  were  disrupted,  so  sick  animals 
could  not  be  moved  and  the  long  and  arduous  marches,  shortage  of  forage,  saddle 
galls  and  injuries  from  sharp  bamboo  stumps,  poisonous  plants  and  bombing  raids 
took  a  heavy  toll.  In  addition,  leeches,  mosquitoes  and  other  insects  made 
conditions  very  unpleasant  and  when  the  men  and  animals  reached  the  frontier 
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monsoon  conditions  and  the  spread  of  tropical  diseases  such  as  malaria,  dysentery 
and  surra  meant  that  in  1943  the  Fourteenth  Army  suffered  a  total  of  12,130 
animal  casualties.  (Gardiner,  2006:76) 

Mules  were  critical  to  the  operation  because  of  the  versatility  they  added.  The 
mules  increased  the  hauling  capability  so  that  the  soldiers  could  operate  longer 
without  resupply.  If  air  drops  or  local  food  procurement  failed,  the  Chindits  would 
eat  the  mules  in  emergency  situations.  Mules  also  posed  certain  problems.  They 
required  food,  got  sick,  injured,  or  died.  If  something  happened  to  a  mule,  the 
equipment  was  either  shifted  to  other  mules  and  the  soldiers,  or  left  behind. 
(Gregory,  1987:6) 

The  jungle  habitat  is  also  highly  characterized  by  an  extremely  hot  and  humid 
environment  that  adversely  affects  the  equipment  as  well  as  the  physical  condition  of  the 
soldiers  (Cross,  2007).  The  hot  and  damp  tropical  climate  shortens  the  life  of  material 
objects  and  lessens  the  efficiency  of  machines.  Untreated  meat  has  to  be  eaten  the  same 
day  that  it  is  prepared.  The  stitching  in  clothes  rots,  rust  appears  overnight  and,  in 
monsoons,  fungus  grows  on  leather  between  dusk  and  dawn.  Prickly  heat  rashes  make  the 
body  uncomfortable,  tempers  short  and  the  wearing  of  equipment  a  burden,  if  not  an 
impossibility.  Everything  is  affected  in  the  tropics:  food,  matches,  cigarettes,  batteries, 
electrical  equipment  and  munitions.  For  example,  in  the  jungle,  weapons  tend  to  rust 
quickly,  and  must  be  cleaned  and  oiled  more  frequently  than  in  most  other  areas,  canvas 
items  rot  and  rubber  deteriorates  much  faster  than  in  more  temperate  areas:  battery  life  is 
shorter  than  is  normal,  electrical  connections  corrode  quickly,  lenses  and  dials  become 
quickly  fogged  with  internal  moisture,  and  troops  drink  more  water,  requiring  greater 
water  purification  and  means  of  transport.  The  following  lines  exemplify  these 
quandaries. 

The  war  was  fraught  with  problems  from  the  beginning.  The  climate  was  hostile 
and  extremely  difficult  to  deal  with.  It  caused  health  problems  for  many  people 
and  added  to  the  medical  service  requirements.  Further,  the  climate  raised  havoc 
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with  supplies.  Rain  was  often  overwhelmingly  heavy.  Humidity  was  constantly 
high.  Corrosion,  mildew,  and  rot  were  principal  problems  for  most  all  classes  of 
supply.  Mold  attacked  everything  and,  in  some  instances,  totally  destroyed  the 
basic  item.  For  example,  combat  boots  and  field  shoes  seemed  to  be  eaten  rapidly 
by  mold.  Insects  and  rodents  were  constant  annoyances  and  caused  considerable 
storage  problems.  It  was  necessary  for  the  United  States  to  construct  huge 
quantities  of  covered  warehousing  to  provide  minimal  supply  protection.  Several 
large  depot  facilities  were  ultimately  constructed  by  the  Services  to  do  higher 
level  maintenance  and  provide  adequate  supply  support.  Even  the  best  of  storage 
seemed  unable  to  defeat  the  climate  and  the  bugs.  Supply  losses  to  the 
environment  continued  through  our  stay  in  the  country.  (Scott  and  other, 
2000:335) 

New  to  jungle  warfare,  the  division  lacked  even  the  basics  for  survival,  prompting 
one  military  historian  to  label  the  soldiers  of  the  32nd  the  “guinea  pigs”  of  the 
South  Pacific.  Men  were  issued  any  of  the  specialized  clothing  that  later  became 
de  rigueur  for  the  war  in  the  South  Pacific.  For  camouflage,  their  combat  fatigues 
were  hastily  dyed  before  left  Australia.  In  the  rain  and  extreme  humidity,  the  dye 
ran  and  clogged  the  cloth,  causing  men  to  develop  horrible  skin  ulcers.  Soldiers 
were  forced  to  wade  through  vines,  creepers,  brush,  dense  stands  of  razor-sharp 
kunai  grass,  and  elephant  grass  as  high  a  basketball  rim  without  the  aid  of 
machetes.  They  did  not  even  have  insect  repellent  -  astonishing  when  one 
considers  that  they  were  fighting  in  a  bug-ridden  place.  They  were  not  equipped 
with  waterproof  containers  either.  Matches  were  often  unusable.  Quinine  and 
vitamin  pills,  salt  and  chlorination  tablets  got  wet  and  crumbled  in  their  pants 
pockets.  Yet,  in  New  Guinea,  the  32nd  Division  was  asked  to  do  the 
extraordinary.  (Campbell,  2008 :XV) 

As  a  result,  to  maintain  men  in  fighting  condition,  the  meals  served  during  the 
jungle  operations  include  a  source  of  fresh  water,  proteins  and  carbohydrates  because 
there  is  a  tremendous  loss  of  energy  in  that  environment  (Bull,  2007:44)  (figures  14  and 
15).  In  general,  beyond  medications,  a  good  supply  of  vitamins  and  energizers  are  carried 
to  camp  as  well  as  a  large  amount  of  salt  (Thompson,  L.,  1994:24).  These  items  are 
required  after  exhausting  marches  or  long  periods  without  sleep.  Jungle  warfare  is  a 
challenge  to  the  military  logistician  because  of  storage,  transport  and  life  time  constraints 
of  these  items  and  after  supplies. 
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Every  one  of  us  was  sick,  with  malaria,  diarrhea,  or  both,  and  we  were  all 
physically  weak  from  lack  of  food  and  sleep.  Back  in  those  days,  the  marines 
didn’t  have  “talking  doctors”  or  psychiatrists,  and  we  didn’t  know  the  meaning  of 
depression.  Everybody  in  every  one  of  our  platoons  was  in  the  same  physical  and 
mental  condition.  We  lived  on  several  tablespoons  of  captured  Jap  rice  per  day  - 
that  is,  if  anyone  brought  it  out  to  us.  Our  water  situation  was  about  as  bad.  Not 
enough  of  it,  and  sometimes  it  wasn’t  available  at  all.  (Marion,  2004: 140) 


Figure  14:  “Battle  Fatigue:  a  British  soldier 
is  inter-viewed  by  a  doctor  of  154th  Field 
Ambulance,  RAMC,  in  Burma.”  (Bull, 
2007:5) 


Figure  15:  “A  May  1969  photo  showed  a 
wounded  United  States  paratrooper 
waiting  for  medical  evacuation  at  base 
camp  in  the  A  Shau  Valley  near  the  Laos 
border  in  South  Vietnam  during  the 
Vietnam  War”. 


From  an  operational  perspective,  short  fields  of  observation  and  fire  -  and  thick 
vegetation  make  engaging  the  enemy  difficult.  Concerning  that  issue,  Clausewitz  made 
one  single  and  important  statement:  “What  is  more,  in  the  depth  of  the  forest  he  will 
hardly  be  in  a  position  to  impress  the  omnipresent  enemy  with  the  superior  weight  of  his 
numbers.  This  is  without  doubt  one  of  the  worst  situations  in  which  an  attacker  can  find 
himself’  (Clausewitz,  1984:154).  The  same  factors  reduce  the  effectiveness  of  indirect 
fire  and  make  jungle  combat  primarily  a  fight  between  infantry  forces  (Asprey, 
2002:844).  As  a  result,  combat  in  the  jungle  normally  takes  place  at  a  close  range  (Bull, 
2007:5)  and  demands  a  high  consumption  rate  of  munitions.  Troops  waste  more  shots 
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than  normal  to  hit  enemy  forces.  In  addition,  in  jungle  warfare,  small  groups  of 

combatants  often  use  unconventional  mobile  tactics,  such  as  infiltration,  ambush,  raids, 

and  guerrilla  operations,  and  their  profound  knowledge  of  the  terrain  to  combat  a  less 

mobile  formal  army.  These  tactics  include  a  large  variety  of  very  aggressive  surprise 

attacks  on  transportation  routes,  installations  and  fixed  structures  (Bull,  2007:18;  DA, 

2009:133;  Moreman,  2002:202;  Sandler,  2002:346-348).  The  American  Army  soldiers 

confronted  this  same  scenario  in  the  South-East  Asian  theatre  of  operations,  during  the 

Vietnam  War,  as  evidenced  in  the  passage  below. 

In  the  jungle,  the  advantage  lies  with  the  force  that  waits.  A  force  moving  at  night, 
especially  a  sizeable  one,  cannot  avoid  making  noise  regardless  of  the  stealth 
attempted.  Therefore,  as  demonstrated  by  the  Americans  in  Vietnam,  ambushes 
are  probably  the  most  effective  means  of  interdicting  an  enemy  who  uses  the 
night  to  transport  supplies  or  infiltrate  to  and  from  an  area.  (Bennett,  1993:42) 

Those  unconventional  tactical  specificities  of  jungle  warfare,  which  require  skills 

and  stealth  techniques,  as  well  as  nimbleness,  swiftness,  and  readiness  on  the  part  of  the 

jungle  combatants,  restrict  severely  the  amount  of  provisions  that  they  can  carry, 

accentuating  the  importance  of  a  continuous  supply  flow  from  military  sources.  The 

extract  below  may  offer  a  better  comprehension  of  these  circumstances. 

Each  item  of  clothing  and  equipment  must  be  considered  in  terms  of  its  necessity, 
and  serviceability  in  jungle  environment.  Lightness  of  weight  is  essential  because 
of  the  difficulties  of  transport.  Serviceability  is  essential,  because  of  the  problems 
of  resupply.  Every  effort  should  be  made  to  reduce  to  the  minimum  the  amount  of 
equipment  to  be  used,  but  care  should  be  exercised  that  no  essential  items  are 
omitted.  During  training,  men  should  be  required  to  use  only  items  of  clothing 
and  equipment  which  will  be  taken  into  combat  areas;  not  only  will  this  teach 
them  to  live  with  a  minimum  of  essentials,  but  also  it  may  indicate  the  non¬ 
essentiality  of  some  things  originally  thought  to  be  necessary  (WD,  1944:35-36). 

Under  the  environmental  constraints  as  well  the  operational  peculiarities  of  jungle 
warfare  detailed  above,  and  considering  that  60%  of  the  total  Amazon  rainforest  is  in 
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Brazilian  territorial  area,  combat  logistics  planners  should  consider  the  use  of  water 
transport  to  support  deployed  troops  during  both  real  and  training  operations.  That 
perspective  emphasizes  taking  advantage  of  the  rich  regional  natural  resources  -  the 
Amazon  basin  -  to  move  troops  and  supplies.  The  Brazilian  Doctrine  for  Joint  Operations 
recognizes  the  fundamental  and  beneficial  effects  of  the  Amazon  basin  for  unity  of  effort, 
synchronization  and  integration  of  military  operations  in  time  and  space  to  successful 
military  campaigns. 

The  large  Amazon  watercourses  -  brooks,  rivers  and  creeks  (figures  16  and  17)  - 
are  a  very  important  part  of  a  transport  system,  especially  in  remote  areas.  Large  rivers 
often  allow  small  ships  and  large  boats  to  penetrate  several  hundred  miles  inland 
(McClain  and  others,  2001:17).  Where  smaller  streams  branch  out,  a  unit  may  establish 
transfer  points  for  the  transloading  of  cargo  into  smaller  watercraft.  A  river  that  is 
normally  very  shallow  during  the  dry  season  will  be  deeper  during  the  rainy  season, 
permitting  travel  by  larger  craft. 
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Because  of  both  the  large  number  of  waterways  in  the  area,  and  the  characteristics 
of  cargo-capacity,  accessibility,  reduced  noise  and  exposure  level,  particularly  during 
night-operations,  that  are  provided  by  river  vessels,  water  transport  is  the  most  reasonable 
option  to  sustain  those  deployed  forces. 

However,  besides  the  multitudes  of  natural  and  artificial  obstructions  in 
watercourses  (topographic  features,  crossings,  etc),  and  the  hydrological  characteristics  of 
rivers  and  streams  (flooding  cycle,  flow  rate  and  direction,  etc),  all  of  which  affect  river 
navigability,  selecting  optimal  routes  and  schedules  that  reduce  enemy  detection  (aircraft, 
helicopters,  and  unmanned  aerial  vehicles),  while  providing  sufficient  support  to  fighting 
forces,  makes  logistics  to  support  jungle  warfare  operations  a  challenging  job  (figures  18 
and  19). 


Figure  18.  Brazilian  UH-60  Blackhawks 
over  the  Amazon  River  in  operational 
training. 


Figure  19.  Brazilian  Navy’s  fast  patrol 
craft  employed  to  fast  deployment  during 
riverine  operations. 


Problem  Statement 

A  small  deployed  combat  force  (Special  Force  unit,  command  troop,  camouflaged 
airfield,  etc)  requires  time-dependent  resupply  in  a  hostile  jungle  environment,  under  the 
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constraints  of  exposure  to  enemy  fire  (aircraft,  helicopters,  and  unmanned  aerial 
vehicles).  Planners  must  select  among  the  available  watercourses  those  whose 
surroundings  provide  stealthy  navigation  to  the  combat  force  location  where  the 
requested  supplies  are  safety  unloaded. 

Research  Focus 

The  context  of  the  research  is  related  to  combat  logistics  (supplying  under  combat 
conditions),  involving  both  scheduling  and  routing  methods,  in  theater-level  modeling. 

The  boundary  of  the  research  is  related  to  supporting  jungle  warfare  operations 
developed  in  tropical  and  subtropical  moist  broadleaf  forests.  This  includes  ones  that 
occur  in  a  belt  around  the  equator  such  as  the  largest  areas  in  the  Amazon  basin  of  South 
America,  the  Congo  basin  of  central  Africa,  Indonesia,  and  New  Guinea.  This  is  due  to 
the  fact  that  these  locations,  between  10°N  and  10°S  latitude  at  elevations  below  3,000 
feet,  share  the  following  both  climatic  and  topographic  characteristics: 

a)  Their  annual  rainfall  rates  of  more  than  3,600  milliliter  (McClain  and  others,  2001:22), 
over  220  inches,  as  well  as  their  temperature  ranges  fairly  constantly  between  72°F  by 
night  and  80°F  by  day  (Louis  and  Draffen  2005:64),  resulting  in  a  daily  flooding 
cycles,  which  change  the  major  river  attributes  (speed,  depth  and  width); 

b)  Their  riverine  forest,  whose  average  tree  height  extends  from  30  feet  to  180  feet 
(Asner  and  others,  2002:486),  can,  depending  on  the  time  of  day,  produce  large 
shadow  areas,  wide  enough  to  hide,  completely  or  partially,  small  boats  traveling  on 
theses  rivers  (figures  20  and  21). 
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This  research  can  be  extended  with  restrictions  to  military  operations  carried  out 
in  parts  of  the  Afrotropic  (equatorial  Africa),  Indomalaya  (parts  of  the  Indian 
subcontinent  and  Southeast  Asia),  the  Neo tropic  (northern  South  America  and  Central 
Asia),  Australasia  (eastern  Indonesia,  New  Guinea,  northern  and  eastern  Australia),  and 
Oceania  (the  tropical  islands  of  the  Pacific  Ocean),  which  share  the  same  biome. 


Research  Objectives 

Develop  a  planning  methodology  that,  given  a  time-window  period  (complete 
timeframe  data),  determines  route  and  schedule  (start  time)  to  enable  a  water  transport  of 
provisions  to  a  fighting  unit  deployed  in  jungle,  in  a  way  that  reduces  its  chance  of 
detection  by  aerial  enemy  searchers.  In  other  words,  plan  the  most  stealthy  water  route 
from  origin  to  destination,  by  taking  advantage  of  the  environment. 


16 


Theoretical  Lens 


The  fundamental  of  theory  of  war  is  used  as  the  perspective,  or  dominant 

paradigm,  to  guide  this  research.  Combat  logistics  must  to  be  consistent  with  the  nature, 

scope  and  structure  of  war,  as  well  the  realities  of  the  battlefield  (Kress,  2002:10).  This 

assertion  is  attested  in  the  following  passage. 

Even  if  logistics  affected  merely  the  ability  of  armies  to  seize  advantages  in 
individual  battles,  that  fact  alone  would  make  this  topic  crucial  to  the  dynamics  of 
combat.  However,  the  effects  of  initiative  reach  far  higher.  Warfare  is  s 
phenomenon  which  can  manifest  itself  in  many  different  ways.  There  are  wars  of 
attrition,  wars  of  maneuver,  wars  of  position,  and  guerrilla  wars,  to  name  but  four. 
If  one  side  consistently  can  choose  the  manner  of  fighting  in  individual  battles, 
that  side  is  well  on  its  way  to  choosing  the  overall  character  of  the  war.  (...) 
Different  nations  prefer  to  wage  different  types  of  wars.  Indeed,  such 
‘preferences’  are  often  a  matter  of  survival.  The  style  and  tempo  of  combat 
determine  what  warfare  will  demand  from  the  material  resources  of  a  nation,  the 
political  base  of  a  government,  the  technological  sophistication  of  an  army,  de 
diplomatic  position  of  a  state,  the  tactical  imagination  of  an  officer  corps  and  the 
moral  sensibilities  of  a  people.  One  might  note,  in  passing,  that  different  styles  of 
warfare  go  hand  in  hand  with  different  styles  of  logistics.  (Kane,  2001 :9) 

Although  it’s  almost  impossible  to  numerically  assess  the  absolute  influence  of 
the  military  supply  system  on  combat  effectiveness  (Dupuy,  1985:38),  it  definitely 
accounts  for  all  activities  in  war  that  are  pre-conditional  to  the  use  of  the  fighting  forces 
(Haldi,  2002;  Newell,  1991:100;  Thorpe,  1982:67).  As  a  result,  logistics  emerges  as  a 
condition  possibility  to  conduct  war.  It  is  a  tactical  and  strategic  concern  for  engagement 
of  war  (Eccles,  1981:12-22;  Newell,  1991:101-120;  Scott  and  others,  2000:32;  Thorpe, 
1982:10-11).  The  purpose  and  scope  of  logistics  is  defined  by  the  needs  of  tactics  and 
strategy,  which,  by  inference,  connects  both  logistics  and  the  theory  of  war,  implying  still 
that  they  share  the  same  attributes  (Coffin,  1993:41;  Guarino,  1992). 
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To  sum  up,  what  is  asked  for  logistics  and  what  it  is  able  to  provide  are 
profoundly  influenced  by  the  particular  circumstances  of  a  war,  campaign,  or  battle 
(O’Hanlon,  2009:146-147;  Rainer  and  others,  2003;  Scott  and  others,  2000;  Thompson, 
J.,  1994:133-219;  Tucker  and  Roberts,  2005:708). 

Methodology 

Mathematical  programming  is  the  methodology  applied  in  this  study.  Search  and 
detection  theory  (light  and  sound  military  doctrine)  is  the  theoretical  framework  for  the 
study. 

Assumptions 

Assumptions  in  this  research  are: 

a)  The  battlefield  can  be  described  as  a  network  of  a  finite  collection  of  nodes,  which 
represent  known  locations  of  a  deployed  combat  unit; 

b)  Arcs  represent  possible  paths  between  the  nodes; 

c)  Combat  unit  for  support  is  known  since  it  is  derived  from  the  basic  battle  plan; 

d)  Support  unit  capacity  is  unlimited; 

e)  The  effects  of  weather,  equipment  failures,  enemy  interdiction  on  support  unit,  and  the 
“fog  of  war”  are  not  considered; 

f)  The  resupply  entity  doesn’t  stop  at  any  point  along  the  path  other  than  the  unloading 
location.  It  doesn’t  wait  at  that  place  after  delivering  its  cargo.  It  moves  in  a  fairly 
unknown  manner,  since  it  must  to  keep  on  the  watercourse  (network),  without, 
exhibiting  any  evasive  action; 
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g)  The  key  characteristics  of  the  routes  (watercourses)  are  risk  (exposure  to  enemy  fire 
level)  and  flow  speed  and  depth  can  change  independently  of  one  another  through  the 
day  but  are  kept  constant  once  movement  begins. 

Limitations 

Limitations  of  this  study  are: 

a)  The  dynamic  nature  of  troop  movements  inside  the  jungle  during  real  operations  as 
well  as  the  radio  silence  imposed  by  wartime  orders  can  create  difficulty  in  finding  the 
precise  force  locations.  Therefore,  a  precise  location  of  all  combat  units  is  improbable; 

b)  Unlimited  support  unit  capacity  is  unrealistic; 

c)  Combat  units  consume  food,  water,  drugs,  and  ammunition  during  each  time  period. 
Quantities  may  vary  depending  on  the  intensity  of  conflict  and  other  concerns. 
Therefore,  the  exact  consumption  rates  under  combat  conditions  are  difficult  to  obtain; 

d)  Chaos  (chance  and  uncertainty)  and  total  friction  (resistance  and  antagonism)  are 
essential  elements  of  the  battlefield,  but  aren’t  modeled; 

e)  Watercourse  features,  like  depth,  width,  and  flow  speed  are  not  constant  through  the 
day,  and  could  change  slowly  or  abruptly  depending  on  water  system  behavior; 

f)  The  risk  level  calculation  which  considers  an  omnipresent  enemy  exacerbates  the 
potential  sources  of  mission  non-fulfillment  creating  a  pessimistic  scenario,  which  is, 
by  definition,  very  improbable. 
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Implications 

Research  implications  are: 

a)  This  research  will  offer  a  jungle  campaign  or  theater-level  model  which  will  support 
military  planning  decisions  regarding  allocation  of  materiel,  positioning  of  combat 
units,  dimensioning  of  support  installations,  scheduling  of  sorties,  and  operational  and 
tactical  needs-derived  decisions; 

b)  This  research  will  increase  the  application  of  mathematical  programming  methods  in 
defense  analysis,  particularly  regarding  the  impact  of  both  combat  and  environmental 
restrictions  in  the  logistics  effectiveness  of  sustaining  military  operations  in  the  jungle, 
in  different  potential  scenarios; 

c)  The  model  developed  in  this  research  can  provide  insights  into  new  and  distinct 
organizational,  doctrine,  training  and  materiel  designs  to  deal  with  the  specificity  of 
supporting  jungle  warfare  units,  in  both  real  and  operational  exercises; 

d)  The  model  developed  in  this  research  will  contribute  relevant  data  concerning 
compatibility,  feasibility  and  operational  cost  of  tactical  movements,  intensity  of 
combat  actions,  as  well  as  estimating  the  readiness  level  of  logistics  systems,  thus 
helping  tacticians  to  develop  their  deployment  plans. 
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II.  Literature  Review 


Introduction 

Theoretical  perspectives  and  previous  research  findings  about  “minimum- 
exposure  moving  process  analysis”  and  taking  advantage  of  surroundings  are  fairly  rare 
in  military  scientific  literature.  Excluding  the  large  amount  of  studies  published  about  the 
Ho  Chi  Minh  trail,  the  main  military  approach  is  focused  on  both  development  and  use  of 
technologically  advanced  materials  that  provide  very  weak  radar  return  or,  more 
precisely,  ones  that  reduce  radar  cross  section.  In  other  words,  it’s  centered  on  military 
technologies  that  are  intended  to  make  vehicles  nearly  invisible  to  enemy  radar  or  other 
electronic  detection.  It’s  also  quite  considered  a  sub-discipline  of  military  tactics  as  well  a 
passive  electronic  countermeasure  (Rao  and  Mahulikar,  2002:629-641). 

However,  in  most  of  the  English  lexica,  the  term  “stealth”  is  perceived  as  an  act, 
behavior  or  procedure  which  intends  to  make  something  imperceptible,  instead  of  an 
engineering  system.  For  instance,  the  Cambridge  Advanced  Learned  Dictionary  defines 
such  expression  as  “the  movement  which  is  quiet  and  careful  in  order  not  to  be  seen  or 
heard,  or  a  secret  action”  (Cambridge,  2008),  while  the  Merriam-Webster’s  Collegiate 
Dictionary  characterizes  it  as  “the  act  or  action  of  proceeding  furtively,  secretly,  or 
imperceptibly”  (Merriam-Webster,  2003).  Keeping  the  same  viewpoint,  the  Oxford 
Advanced  Learner’s  Dictionary  of  Current  English  describes  the  label  stealth  as  “the  fact 
of  doing  something  in  a  quiet  or  secret  way”  (Hornby,  2007).  As  well,  to  the  Chambers 
Dictionary,  defines  stealth  as  “a  secret  or  unobtrusive  going  or  passage  -  procedure  or 
manner  -  furtiveness  -  (being  able)  to  approach  (an  enemy)  without  being  detected” 
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(Chambers,  2008).  Based  on  the  explanations  presented  above,  it’s  possible  to  affirm  that 
vocable  stealth  can  be  understood  as  the  ability  to  move  with  a  minimum  exposure  from 
an  observer.  This  research  employs  these  lexical  approaches  in  the  military  supply 
movement  analysis  in  jungle  combat  conditions. 

It’s  common  sense  that  the  ability  to  operate  and  maneuver  safely  in  a  battlefield 
environment  using  low-detection  footprint  technique  greatly  increases  the  chance  of 
success  of  the  operational  and  tactic  objectives.  However,  using  shadows  to  conceal 
movement  is  uncommon  in  the  current  scientific  military  writings,  being  found  more 
frequently  on  the  Computer  Science  domain,  especially  in  robotics.  In  the  available 
theoretical  papers  analyzed,  the  emphasized  focal  point  is  to  minimize  the  non-detection 
probability  of  detecting  a  move,  instead  of  deceiving  a  fixed  object  (Dobbie,  1973:907). 

In  this  research,  the  word  “detection”  is  defined  as  “the  act  to  gathering 
information  pertaining  to  the  object  being  sought,  the  sifting  out  of  what  is  important 
information  and  relaying  on  that  information  in  some  efficient  form  to  decision  maker” 
(Moore,  1970:3).  Therefore,  it  infers  the  existence  of  observer(s)  capable  of  identifying 
and  distinguishing,  among  wide-ranging  physical  entities,  patterns  which  can  produce 
behavior  changes  in  that  viewer.  In  the  strictly  military  perspective,  that  approach  implies 
that  enemy  detection  can  induce  hostile  course  of  actions  such  as  damage  and/or 
destructive  acts,  which  would  not  occur  if  any  kind  of  severe  secret  behavior  was 
conducted.  In  most  of  the  research  from  civilian  sectors,  the  major  concern  is  optimizing 
performance  instead  of  minimizing  impact  of  physical  (and  human)  loss.  The  following 
studies  identify  proven  benefits  of  concealment  during  operational  missions. 
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Military  Perspectival  Studies 

Among  the  rare  military  writings  about  stealth,  the  analysis  of  the  Ho  Chi  Minh 
trail  are  the  most  prominent.  The  Ho  Chi  Minh  trail  was  a  logistical  system  ran  from  the 
Democratic  Republic  of  Vietnam  into  the  Republic  of  Vietnam,  through  the  neighboring 
kingdoms  of  Laos  and  Cambodia  (Thompson,  1990:134).  It  was  not  so  much  a  single 
route  but  a  network,  a  honeycomb  of  routes,  passing  through  country  that  was  alternately 
limestone  karst,  triplecanopy  jungle,  and  grassland.  In  general,  it  started  with  the 
transportation  arteries  in  North  Vietnam,  swung  west  into  Laos,  south  to  the  South 
Vietnamese  border,  and  at  various  points  crossed  back  to  the  east  and  into  South 
Vietnam.  The  Laotian  part  of  the  system  continued  further  south  into  Cambodia  and 
intersected  with  a  network  there  which  was  known  as  the  “Sihanouk  Trail”  (Vongsavanh, 
1981:14). 

The  importance  the  Ho  Chi  Minh  Trail  is  its  stealth  structural  features  as  well  as 
its  congruent  internal  management.  Its  singular  characteristics  included  the  use  of 
“camouflaged  vehicles”  (figure  22),  the  “invisible  roads”  dug  out  of  riverbeds  during  the 
dry  season  and  used  at  night  so  as  not  to  leave  convoy  tread  marks,  “submarine  bridges” 
built  just  under  the  waterline  of  a  river  to  be  invisible  from  the  air,  over  which  convoys 
crossed  with  foot  guides,  as  well  as  the  evolution  of  the  trail  from  a  footpath  (figure  23) 
to  an  eventual  logistical  highway  designed  with  truck  stops. 
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Figure  22.  A  camouflaged  truck  towing  a  122mm  gun  in 
February  1973  (Nalty,  2005:290). 


Figure  23.  Bicycles  laden 
with  supplies  in  the  Ho 
Chi  Minh  Trail. 


The  routes  consisted  mostly  of  small  trails,  practically  undetectable  from  the  air, 
for  personnel  movement  and  roads  for  vehicles.  Personnel  generally  walked  or  pushed 
bikes  along  narrow  foot-paths  (figures  24  and  25).  These  personnel  were  both  porters 
carrying  war  supplies  and  soldiers  going  south  to  fight.  The  network  (for  both  vehicle  and 
foot  traffic),  on  the  other  hand,  consisted  of  the  trails,  themselves  connected  by  a  series  of 
small  rest-points,  larger  storage  sites,  and  a  few  major  base  areas  manned  on  a  permanent 
basis  (figure  25).  The  network,  therefore,  needed  not  only  soldiers  to  operate  the  routes, 
but  also  personnel  to  feed  and  support  the  workers. 

In  terms  of  furtiveness  and  unobtrusiveness  aspects,  the  Ho  Chi  Minh  Trail  is 
perceived  by  the  historians  as  the  most  secret  military  logistics  system  in  the  history  of 
armed  conflict.  The  following  passage  confirms  the  success  of  the  Ho  Chi  Minh  trail 
enterprise. 

In  the  world  of  logistics,  there  are  a  few  brand  names  to  match  of  the  Ho  Chi 
Minh  Trail,  the  secret,  shifting,  piecemeal  network  of  jungle  roadways  that  helped 
the  North  win  the  Vietnam  War.  From  the  air  the  Ho  Chi  Minh  was  impossible  to 
be  identified  and  although  the  United  States  Air  Force  tried  to  destroy  this  vital 
supply  line  by  heavy  bombing,  they  were  unable  to  stop  the  constant  flow  of  men 
and  logistical  supplies.  The  miracle  of  the  Ho  Chi  Minh  Trail  “logistic  highway” 
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was  that  it  enabled  the  “impossible”  to  be  accomplished.  A  military  victory  is  not 
determined  by  how  many  nuclear  weapons  can  be  built,  but  by  how  much 
necessary  materiel  can  be  manufactured  and  delivered  to  the  battlefield.  The  Ho 
Chi  Minh  Trail  enabled  the  steady,  and  almost  uninterrupted,  flow  of  logistics 
supplies  to  be  moved  to  where  it  was  needed  to  ultimately  defeat  the  enemy. 
(Taylor,  2007:1-7) 


Figure  24.  “A  group  of  walking  Figure  25.  A  Serial  of  camouflaged  supply  depots 
North  Vietnamese  on  the  Ho  Chi  set  up  through-out  the  Ho  Chi  Minh  trail.  (Nalty, 
Minh  trail”  (Nalty,  2005:284).  2005:285). 


Civilian  Perspectival  Studies 


On  the  first  civilian  essay  studied,  a  group  of  researchers  from  the  University  of 


Southern  California  developed  an  algorithm  using  the  stealth  strategy  to  improve 


performance  tasks  in  multi-robot  assignments.  The  key  idea  was  to  reduce  both  resource 


conflict  and  physical  interference  created  among  the  machines,  since  beyond  some 


critical  limit,  any  additional  robot  decreases  overall  performance.  Keeping  as  low  a 


profile  as  possible,  the  robots  could  carry  out  their  tasks  in  the  presence  of  single  or 


multiple  observers,  which  included  locating,  corralling,  surrounding,  and  navigation 
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(figures  26  and  27).  In  that  study  the  researchers  assumed  an  unlimited  omnidirectional 
sensing  observer.  The  robots  carry  out  their  traverses  one  at  a  time,  and  the  environment 
consists  of  objects  that  can  occlude  the  robots  from  the  observer's  sensors.  Also,  the  team 
initially  had  no  map  of  the  environment  but  the  positions  of  the  observer  and  the  goal  are 
known. 


According  to  the  final  experimental  results,  conducted  both  in  simulation  and  in  a 

real  outdoor  environment,  the  researchers  demonstrated  the  algorithm’s  versatility  in 

taking  advantage  of  an  environment  that  changes  between  robot  traverses.  Their 

outcomes  also  revealed  that  surrounding  configurations  allow  the  robustness  of  the 

approach  to  be  evaluated  in  terms  of  low-visibility  path  selection,  advantage  of  sharing 

information  to  reduce  cumulative  visibility,  repeatability,  and  reactivity  to  a  changing 

environment.  The  method  can  be  briefly  summarized  in  the  following  statements. 

In  our  approach,  the  robots  carry  out  their  traverses  one  at  a  time,  sequentially, 
and  generate  an  occupancy  grid  representation  of  the  environment  en  route.  The 
occupancy  grid  is  modeled  by  potential  fields,  along  with  task  specific 
information  such  as  observer  and  goal  position.  The  combination  of  the  fields 
forms  an  abstract  view  of  the  environment  from  which  navigation  waypoints  are 
extracted.  To  take  advantage  of  multiplicity,  each  robot  commencing  a  traverse  is 
provided  with  the  occupancy  grid  and  filtered  path  information  from  the  previous 
robot.  The  filtered  path  is  a  waypoint  list  generated  from  events  that  occurred 
along  the  robot’s  path.  The  successor  robot  uses  this  information  to  make 
decisions  about  waypoint  selection  and  overwrites  the  provided  occupancy  grid 
with  its  ego-centric  sensor  data.  By  sharing  information,  each  robot  follows  a 
lower-visibility  path  than  its  predecessor,  and  in  the  case  of  static  environments, 
the  paths  are  traversed  more  efficiently.  (Tews  and  others,  2004) 
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Figure  26.  “An  example  of  stealthy  and 
non-stealthy  traverses  shown  on  the  global 
potential  field.  The  irregular  white  objects 
are  the  barricades.  Dark  areas  represent 
shadows  behind  the  barricades  and  there¬ 
fore  attractive  locations  for  the  robot.  The 
grey  area  surrounding  the  stealthy  path  is 
the  accumulated  effect  of  a  local  high¬ 
valued  potential  field  that  is  positioned  at 
the  robot’s  location  during  the  traverse  to 
prevent  waypoints  being  selected  too  close 
to  it.”  (Tews  and  others,  2004) 


Figure  27.  “The  shadow  regions  are  cast 
behind  objects  from  the  observer’s 
position.  Safety  zones  are  regions  inside 
shadows  for  the  robot  to  manoeuvre 
without  being  exposed  by  accidentally 
overrunning  a  frontier”  (Tews  and  others, 
2004). 


Under  the  same  label,  another  group  of  researchers  from  the  University  of 
Southern  California  studied  stealthy  traverses  using  artificial  intelligence  devices.  Robot 
is  tasked  with  moving  from  one  location  to  another  while  remaining  hidden  from  an 
observer  at  a  known  location.  Since  the  researchers  assume  that  robots  have  no  priori 
model  of  the  environment,  a  stealthy  transverse  behavior  makes  opportunistic  use  of 
terrain  features  to  hide  from  the  observer  (figures  28  and  29).  To  perform  that 
experiment,  the  researchers  developed  and  employed  a  behavior-based  algorithm  which 
combines  a  number  of  well  known  techniques,  such  as  potential  fields,  ray-tracing  and 
connected  components  (Birgersson  and  others,  2003). 

According  to  their  research,  the  robot  used  a  pair  of  two  2-dimensional  internal 
maps  to  represent  the  world:  one  map  describes  the  location  of  obstacles,  while  the  other 
describes  the  world's  observability  (i.e.,  which  areas  can  or  cannot  be  seen  by  the 
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observer).  These  maps  induced  an  artificial  potential  field  that  was  used  to  guide  the 
machines  towards  the  goal,  away  from  obstacles,  and  towards  unobservable  regions. 
These  experiments  show  a  clear  improvement  in  the  stealthiness  of  the  robot,  when 
evaluated  in  both  real  and  simulated  experiments,  comparing  it  against  a  regular  goal- 
seeking/obstacle-avoidance  behavior.  Although  the  algorithm  does  still  have  some 
weaknesses,  due  to  resulting  primarily  from  local  minima  in  the  potential  fields,  the 
experiments  confirmed  that  the  algorithm  does  indeed  produce  stealthy  behavior. 


Figure  28.  This  picture  above  shows  a 
real  robot  movement  which  applied  the 
sneaky  behavior.  The  point  number  1 
represents  the  start  position;  an  observer 
is  placed  on  the  point  2;  and  the  point  3 
denotes  the  goal  spot. 


Figure  29.  Above,  another  layout  pres¬ 
ents  a  path  following  by  a  stealth  be- 
haveior-based  robot.  As  the  figure  26, 
points  1,  2  and  3  designate  the  start  po¬ 
sition,  observer  location  and  target  point, 
respectively. 


According  to  Ravela  and  others,  the  motivation  for  using  terrain  features: 

(...)  lies  in  the  need  for  robust  stealth  navigation,  i.e.  to  be  able  to  react  to 
uncertainty  in  the  environment  arising  from  inaccurate  terrain  maps,  inaccurate 
sensory  modalities  such  as  odometry  or  to  react  to  a  dynamic  environment  where 
the  target  is  moving.  (Ravela,  1994:1093) 

A  civilian/military  application  of  low-visibility  path  can  be  found  on  the  study 
carried  out  by  three  researchers,  Teng,  DeMenthon,  and  Davis,  all  of  them  from  the 
Department  of  Computer  Science  of  the  University  of  Maryland  Institute  for  Advanced 
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Computer  Studies.  In  their  writings  they  described  a  method  for  solving  visibility-based 
path  planning  problems  over  natural  terrain  using  massively  parallel  hypercube  machines 
(Teng  and  others,  1993).  These  problems  arose  in  the  development  of  both  autonomous 
and  teleoperated  systems  for  vehicle  navigation  in  a  battlefield  scenario,  in  which  a 
typical  example  is  to  find  a  path  that  is  hidded  from  moving  adversaries.  According  to 
their  paper,  that  kind  of  problem,  i.e.  path  planning  with  moving  obstacles,  can  be 
generalized  as  a  time-varying  constrained  path  planning  problem,  and  it’s  proven  to  be 
computationally  hard  (Reif  and  Sharir,  1985).  Their  summarized  problem  is  presented  in 
the  passage  below  (see  figure  30). 

Given  a  digital  terrain  map  with  elevation  data  at  each  grid  cell,  and  information 
about  a  friendly  agent  and  adversary  agents  moving  on  the  ground,  we  want  to 
find  a  path  for  the  friendly  agent  from  its  current  position  to  a  final  goal  such  that 
the  movement  of  the  friendly  agent  is  hidden  from  the  adversaries  by  taking 
advantage  of  the  terrain.  We  refer  to  this  property  as  safety.  (Teng  and  others, 
1993) 

It  is  important  to  emphasize  that  visibility  constraints  are  potentially  harder  than 
obstacle  constraints  since  their  sizes  and  distributions  can  change  in  time  with  little 
coherence,  nearly  randomly.  Also,  the  analysis  of  visibility  requires  intensive 
computation,  showing  that  their  problem  is  potentially  harder  than  path  planning  under 
moving  obstacle  constraints.  However,  by  transforming  these  problems  to  a  discretized 
formalism,  many  basic  computations  can  be  arranged  in  a  regular  pattern,  and  thus  could 
be  done  in  parallel  efficiently.  Therefore,  the  researchers  decided  to  combine, 
simultaneously,  the  three  categories  of  path  planning  (path  planning  with  moving 
obstacles,  for  terrain  navigation,  and  visibility  analysis  on  terrain),  with  digital 
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approximation  techniques  on  a  powerful  parallel  machine  which  could  solve  realistic 
problems  in  an  acceptable  amount  of  time. 
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Figure  30.  “Finding  a  path  through  the  Safety-Reachable  (RS)  corridor;  thick  vertical  line 
segments  show  the  RS  regions  and  the  dashed  line  is  an  example  of  a  valid  path,  in  which 
the  angle  p  represents  the  speed  constraint”  (Teng  and  others,  1993). 

Although  that  kind  of  problem  has  been  proved  to  be  hard,  it  is  nevertheless  a 
very  important  problem.  Researchers  have  developed  algorithms  to  plan  the  best  path 
under  different  assumptions  and  approximations.  For  example,  an  algorithm  proposed  by 
researchers  Kamal  Kant  and  Steven  Zucker  used  a  path-velocity  decomposition  approach 
that  first  finds  a  minimum  length  path  among  the  static  obstacles  and  then  determines  the 
speed  along  the  path  to  avoid  collision  with  the  moving  obstacles  (Kant  and  Zucker, 
1986).  Theoretically,  this  same  approach  can  be  employed  on  the  current  jungle  warfare 
supply  problem,  in  which  the  “path  length”  is  replace  by  “risk  level”  or  “travel 
survivability  potential”. 

When  discussing  detection  actions,  both  in  military  as  well  in  civilian  prospect,  it 
becomes  important  to  refer  to  the  search  theory  foundations.  This  knowledge  provides  the 
most  effective  quest  techniques  when  developing  and  sustaining  effectual  hide-out 
behavior.  In  consequence,  the  tactical  analysis  conducted  in  selecting  stealthy  routes  also 
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involves  comprehension  of  the  main  search  theory  concepts  (Morse  and  Kimball, 
1998:86). 

As  stated  in  many  Search,  Detection  and  Damage  Assessment’s  literature  sources, 
the  efficiency  of  search  is  measured  in  terms  of  the  ability  of  the  observer  to  identify  the 
aim  (Przemieniecki,  2000:261;  Washburn  and  Kress,  2009:133;  Zehna,  2005:125,  271). 
This  will  be  altered  by  the  speed  of  the  searcher  as  well  as  the  target  visibility.  Rapid 
search  or  poor  visibility  reduces  the  likelihood  of  the  target  being  identified.  A  slow 
search  will  reduce  the  area  that  can  be  covered  in  the  time  available. 

The  critical  aspects  of  a  search  problem  are  concerned  with  the  nature  of  the 
target  in  space  and  time,  and  the  method  and  efficiency  of  search.  It  is  convenient  to 
consider  the  target  in  dimensional  terms  as  well  to  be  described  in  terms  of  knowledge  of 
its  movement  patterns  (Zehna,  2005:258).  Also,  the  target  is  assumed  to  have  a  particular 
probability  of  being  in  a  specific  location,  which  is  critical  to  the  eventual  search  methods 
adopted.  For  example,  during  war  campaigns,  ground  supply  convoys  are  more  likely  to 
be  found  along  the  existing  roads  in  theatre  of  operations  instead  of  any  location  far  away 
from  military  service  camps.  From  the  same  prospect,  the  detection  function  also  takes 
into  account  the  observer  attributes,  such  as  its  velocity,  altitude,  available  frequency 
search  time,  as  well  as  the  search  models  (exhaustive,  random,  and  inverse  cube  law) 
chosen  by  camp  commanders  (Przemieniecki,  2000:269). 

Therefore  many  concerns  and  assumptions  arise  about  the  target,  such  as  number, 
position,  movement,  as  well  the  effectiveness  of  the  searcher  using  different  quest 
strategies.  In  general,  the  goal  is  to  increase  the  probability  of  “success”  by  choosing  the 
best  plan  within  the  limits  of  time  or  resources.  In  Table  1,  it’s  possible  to  identify  the 
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major  search  strategy  attributes  which  need  to  be  weighed  for  military  planners.  In  this 
specific  study,  combat  logisticians  conceive  their  war  supply  plan  of  action  before 
performing  on  the  field  of  battle  (Stone,  1975). 


Table  1.  Summary  of  Present  Capability  to  Find  Optimal  Search  Plans  (Stone,  1975) 


Search  Effort 

Stationary 

Moving 

Discrete  Time 

Continuous  Time 

Short  Range 
Infinitely  Divisible 

Uniformly  solvable 
for  concave  detection 
functions 

Solvable  with 
exponential 
detection  function 

Necessary  and 
sufficient 
conditions 

Discrete  Effort 

Solvable  in  principle 

Necessary  but  not 
sufficient 
conditions 

Little  or  no 
progress 

Constraints  on 
Search 

Solvable  only  in 
special  situations 

Some  progress 

Little  or  no 
progress 

Long  Range 

Some  approximate 
solutions 

Little  or  no 
progress 

Little  or  no 
progress 

Comparing  both  the  content  presented  in  chapter  I  which  highlighted  that  jungle 
warfare  is  founded  on  resistance  strategy  “trading  space  for  time”  (Epstein  1985:14).  The 
research  findings  above  infer  that  there  is  not  an  optimal  search  plan  for  the  long  range 
search  effort.  That  conclusion  benefits  the  defender  since  the  aggressor  must  spend  their 
human  and  materiel  resources,  which  are  limited  by  nature  to  locate  their  opponent 
without  any  expectation  of  an  end  result.  Therefore,  from  the  jungle  combat  logistics 
standpoint,  timely  supply  isn’t  a  cut-off  constraint,  and  a  reasonable  time-window  will 
always  exist. 

Under  the  detectability  problems,  the  Department  of  Industrial  Engineering’s 
researchers  from  the  Seoul  National  University  applied  the  theory  of  search  concepts  to 
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improve  the  single-searcher  path-constrained  discrete-time  Markovian-target  search 

problem  solution.  They  proposed  an  algorithm  that  uses  an  approximate  non-detection 

probability  of  a  search  path,  referred  to  as  the  depth-l  approximation,  which  is  computed 

from  the  conditional  probability  reflecting  the  search  history  over  the  time  windows  of  a 

predetermined  length,  /.  The  problem  then  is  formulated  as  a  shortest  path  problem  on  an 

acyclic  layered  network  whose  number  of  layers  is  of  the  order  of  search  duration  T.  The 

objective  is  to  find  a  searcher’s  path  that  minimizes  the  non-detection  probability,  i.e.  the 

probability  that  the  target  survives  the  search. 

A  single  target  moves  among  a  finite  set  of  cells  S  =  {1;  2;  ...;  C}  over  the 
discrete  time  periods  t  =  {1;  2;  .  .  .  ;  T).  At  the  beginning  of  each  time  period,  it 
moves  from  a  cell,  say  i,  to  a  cell,  say  j,  of  its  neighbors,  N(i),  with  transition 
probability  py.  We  adopt  the  hexagon  cells  (see  figures  29  and  30)  instead  of 
square  cells  which  appear  to  enhance  the  model’s  applicability.  In  this  case,  given 
a  cell,  the  neighbors  of  the  cell  are  defined  to  be  the  6  adjacent  cells  and  the  cell 
itself.  Hence  the  number  of  neighbor  of  a  cell  v  =  7  except  the  case  when  the  cell 
is  near  to  the  search  area  boundary.  A  single  searcher  examines  one  cell  during 
each  time  period.  If  the  target  happens  to  be  in  the  same  cell,  say  i,  the  searcher 
detects  it  with  probability  qt=  (1  —  for  some  at>  0.  Our  goal  is  to  find  a 
path  of  the  searcher,  represented  by  the  cells  it  visits  over  the  time  periods,  U;  ij; . 

.  .  ;  if  that  minimizes  the  non-detection  probability.  (...)  When  a  search  path  is 
fixed,  the  target’s  non-detection  probability  can  be  computed  by  conditioning  on 
the  target’s  location  over  the  search  periods.  The  difficulty  of  the  problem  lies  in 
that  the  number  of  the  possible  search  paths  grows  exponentially  as  T  or  C 
increases,  while  the  enumeration  is  unavoidable  in  an  exact  method  due  to  the 
intractability  of  the  problem.  To  cope  with  the  difficulty,  we  propose  a  heuristic 
minimizing  an  approximate,  instead  of  exact,  non-detection  probability  of  a  path. 
(Hong  and  others,  2009:352) 

Though  the  main  dissimilarity  between  the  previous  surveys  are  associated  to  a 
single  searcher  and  its  behavior,  omnidirectional  sensing  stationary  conduct  and  target 
seeker  attitude,  pursue  the  non-sighting  probability  minimization  while  the  (single)  target 
keeps  moving  (figures  3 1  and  32). 
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Figure  31.  A  search  area  with  400  hex-  Figure  32.  An  example  of  heuristic  op- 

agons  cells  (Hong  and  others,  2009:353).  timal  path  (Hong  and  others,  2009:363). 


Synthesis 

As  previously  mentioned,  the  number  of  academic  papers,  empirical  analysis  as 
well  as  virtual  experiments,  using  modeling  and  simulations  approaches,  regarding 
stealthy  conduct  are  small  compared  to  its  technological  research.  However,  among  these 
extracts,  important  and  fundamental  insights  were  gathered  that  built  the  theoretical 
perspective  of  optimal  supply  technique  in  jungle  combat  conditions. 

The  first  conclusion  from  this  review  confirms  the  logical  and  rational  importance 
of  environmental  features  on  the  shadow  making  process,  as  well  as  its  practical 
applicability  in  both  military  and  civilian  sectors.  Further,  the  nonmilitary  research  results 
demonstrate  a  powerful  area  of  study,  in  which  unobtrusive  and  sneaky  movements 
enable  systems  achieving  a  high  performance  level  when  the  mutual  sense  of  presence 
inhibits  total  effectiveness. 
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The  second  inference  is  related  to  the  usefulness  of  setting  up  an  omnidirectional 
sensing  observer  assumption.  This  point  of  view  suggests  a  conservative  posture,  and  will 
guarantee  that  the  final  research  outcomes  will  provide  the  worst  possible  operational 
situation  for  provision  action  plans  on  the  battlefield.  Consequently,  camp  commanders 
as  well  as  their  logistic  planners  identify  their  maximum  loss  in  terms  of  materiel, 
personal  and  equipment  which  improves  both  their  logistics  and  tactical  decisions  and 
their  mission  feasibility.  In  other  words,  based  on  “calculated  risk”  they  match  support 
plan  with  combat  actions  (Eccles,  1981:80). 

The  third  upshot  refers  to  the  evidenced  complexity  and  circuitousness  involved 
in  developing  an  efficient  search  strategy,  which  essentially  means  to  attain  the  intended 
objective  without  wasting  resources  (time,  effort  and  expense).  The  vast  number  of 
options,  in  addition  to  the  occasional  complex  relationship  among  mathematical  variables 
can  produce  a  diversity  of  outcomes,  suggests  a  prudent  and  careful  research  posture  in 
order  to  avoid  setting  out  for  tortuous  and  endless  paths. 
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III.  Methodology 


Introduction 

A  practical  and  reasonable  way  to  carry  out  a  resupply  mission  with  combat  units 
and  provides  low-visibility  to  aerial  observers  while  reducing  the  risk  of  enemy  detection 
during  the  travel  is  to  take  advantage  of  the  surroundings,  especially  the  shadows 
provided  by  large  and  bulky  tropical  riverine  flora.  Choosing  routes  that  offer  large  dark 
areas  upon  the  river  surface,  so  the  boat  can  travel  the  watercourse  network  safely 
improves  the  chance  of  survivability. 

Taking  advantage  of  the  surprise  effect,  of  flexibility  and  adaptability,  peculiar 
attributes  of  irregular  warfare  operations,  are  needed  for  jungle  warriors.  It  becomes 
imperative  that  they  also  are  shared  by  jungle  combat  logisticians  (Gregory,  1987).  In 
order  to  meet  the  demands  imposed  by  forces  fighting  in  a  characteristically  nonlinear 
dynamic  battlefield,  the  logistics  domain  must  go  through  a  commensurable 
transformation;  maybe  a  large  “paradigm  shift”  in  its  attitude  at  the  tactical  level  of  war 
(Asprey,  2002;  Pinkston,  1996).  For  example,  instead  of  transporting  a  large  quantity  of 
provisions,  they  should  plan  shorter  but  regular  loadings,  coherent  with  the  logical  nature 
of  a  high  mobile  style  of  war,  which  eliminates  large  static  stockpiling. 

Similarly,  transportation  cost  or  delivery  time  precision  (punctuality)  are  not  the 
major  metrics  to  measure  the  success  of  a  military  logistics  operation  in  a  jungle  scenario. 
The  basis  of  successful  guerrilla  combat  is  the  resistance  strategy  and  offensive  action 
combined  with  surprise.  The  jungle  combatants  manage  the  time  element  (DA,  2009:10; 
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Trinquier,  2006:70).  Rather,  the  ability  to  load  the  entire  cargo  in  response  to  operational 
needs  is  much  more  important. 

Regarding  the  stealth  efficiency  of  each  path,  the  darkness  magnitude,  depends 
directly  on  the  source  of  light  as  well  the  shadow-casting  object.  In  this  study,  the  former 
is  the  sun,  while  the  latter  is  the  botanic  system  bordering  the  watercourses.  Considering 
that  the  intensity  of  sun  radiation  is  a  function  of  the  solar  direction  relative  to  the  local 
plane  of  the  earth’s  surface  at  that  instant,  and  also  taking  into  account  that  variables  such 
as  solar  azimuth  and  solar  altitude  angles  change  continuously  throughout  the  day,  solar 
radiation  needs  to  be  continuously  updated  while  the  boat  is  moving  on  the  network.  To 
make  simpler  computations,  the  solar  declination  angle  may  be  assumed  constant  and  is 
calculated  only  once  per  day. 

Given  a  triad  of  information  date/time/location  as  well  a  structured  network 
(source,  destination,  path  attributes,  and  existing  connections)  the  best  movement  of  an 
ship  inside  a  temporal  window  may  be  computed  through  three  basic  and  fundamental 
steps: 

First,  set  up  a  solar  position  algorithm  capable  to  provide  shadow  lengths  based 
on  actual  move  position  over  time.  The  position  of  any  place  upon  the  surface  of  the 
Earth  may  be  at  once  determined  when  its  latitude  and  longitude  are  known  (Lawson, 
2008:6).  The  former  gives  its  distance,  north  or  south,  from  the  equator;  the  latter,  its 
distance  east  or  west  of  a  certain  point.  In  England,  the  point  agreed  upon  is  Greenwich. 
In  common  usage,  “latitude”  refers  to  geodetic  or  geographic  latitude  and  is  the  angle 
between  the  equatorial  plane  and  the  direction  of  gravity  at  any  point  to  the  reference 
ellipsoid,  which  approximates  the  shape  of  earth  to  account  for  flattening  of  the  poles  and 
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bulging  of  the  equator.  This  value  usually  differs  from  the  geocentric  latitude,  which,  as 
the  name  implies,  is  the  angle,  at  the  center  of  the  earth,  between  the  plane  of  the  equator 
and  a  line  drawn  from  the  observer  to  that  center.  The  line  evidently  does  not  coincide 
with  the  direction  of  gravity,  since  the  earth  is  not  spherical  (Duffett- Smith,  1988). 
Geocentric  latitude  is  employed  in  astronomical  calculations,  since  it  relates  to  the  moon 
and  eclipses,  in  which  it  becomes  necessary  to  reduce  observations  to  the  center  of  the 
earth  (Young,  1902).  To  simplify  the  calculations,  the  Ptolemaic  view  of  the  sun’s  motion 
around  the  earth  is  used.  In  other  words  the  site  is  fixed  and  is  the  centre  of  origin. 

Second,  determine  the  risk  measure  of  effectiveness  between  each  node,  through 
the  travel  time.  This  information  is  obtained  from  the  study  of  cinematic  behavior  of  the 
sun  through  the  sky,  using  the  astronomical  algorithm  defined  in  the  first  step,  as  well  as 
the  mathematical  elemental  procedures  comparing  both  the  local  path  azimuth  and  the 
local  solar  azimuth,  at  each  particular  time.  Using  trigonometric  identities  with  solar  and 
earthly  attributes,  it’s  possible  to  figure  out  the  available  shadow,  on  each  path  at  any 
period  of  time.  Therefore,  once  the  boat  leaves  its  dock,  which  is  assumed  fix  on  the 
depot  location,  the  total  travel  risk  level  will  depend  exclusively  on  the  time  departure,  as 
well  as  the  path  it  had  taken.  Assuming  that  after  accessing  a  route  the  move  can’t  return 
all  possible  outcomes  need  to  be  previously  tested  in  order  to  define  the  potential  optimal 
route,  if  that  is  possible  and  really  exists.  Route  identification  is  performed  employing  an 
algorithm  that  is  capable  of  minimizing  the  formerly  defined  success  mission  attributes 
called  exposure  to  enemy  fire  level,  or  simply,  risk  level. 

Finally,  set  up  a  graph  algorithm  capable  to  work  on  network  flow  problems, 
particularly  one  distance-time  parameters  minimization  in  route  optimization  studies. 
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Applying  the  elemental  principle  of  isomorphism  -  preservation  of  structural  logic  -  it’s 
possible  to  appropriately  shift  the  above-mentioned  variables  and  solve  the  travel  risk 
assessment  quandary  proposed  in  this  research. 

Manual  Solar  Position  Calculation  Model 

The  solar  position  in  the  sky  relative  to  a  location  on  the  surface  of  the  Earth  can 
be  specified  by  two  angles:  the  solar  altitude  and  the  solar  azimuth  (Iqbal,  1983)  (figure 
33).  The  former  is  the  angle  between  the  sun’s  position  and  the  horizontal  plane  of  the 
earth’s  surface.  At  sunset/sunrise  the  altitude  is  0  and  is  90  degrees  when  the  sun  is  at  the 
little  zenith.  The  latter  specifies  the  angle  between  the  line  from  the  observer  to  the  sun 
projected  on  the  ground  and  the  line  from  the  observer  due  south.  A  positive  azimuth 
angle  generally  indicates  the  sun  is  east  of  south,  and  a  negative  azimuth  angle  generally 
indicates  the  sun  is  west  of  south  (Calow,  1999:137). 


Figure  33.  Major  attributes  used  to  calculate  the  risk  on  each  path  (river). 

The  astronomical  scientific  literature  contains  many  fast  algorithms  for  solar 
position  calculation.  They  require  a  small  computation  effort,  although  their  uncertainties 
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are  greater  than  0.01°  in  solar  zenith  and  azimuth  angle  calculations;  further,  some  of 
them  are  only  valid  for  a  specific  number  of  years.  For  example,  Michalsky’s 
calculations  are  limited  to  the  period  from  1950  to  2050  with  uncertainty  greater  than 
±0.01°(Michalsky,  1988),  and  the  calculations  of  Blanco-Muriel  and  others  are  restricted 
to  the  period  from  1999  to  2015  with  maximal  error  of  0.008°  (Blanco-Muriel  and  others, 
2001).  Some  algorithms  like  the  Spencer  formula  and  the  Pitman  and  Vant-Hull 
algorithm,  cannot  reduce  the  error  more  than  0.25°  and  0.02°  respectively  (Reda  and 
Andreas,  2003),  but  also  present  the  ecliptic  sun  coordinates  based  exclusively  on  its 
declination  and  equation  of  time  (Blanco-Muriel  and  others,  2001),  which  require 
additional  reckoning  to  obtain  both  the  solar  azimuth  and  altitude  angles. 

At  this  time,  the  most  precise  solar  position  calculator  is  National  Renewable 
Energy  Laboratory  Solar  Position  Algorithms  (NREL's  -  SPA),  proposed  by  Ibrahim 
Reda  and  Afshin  Andreas  in  2003.  It  calculates  the  solar  zenith  and  azimuth  angle  with 
uncertainties  smaller  than  0.0003°,  in  the  range  from  0°  and  90°,  on  a  long  time  period 
(2,000B.C.-6,000A.C).  Nevertheless,  it  involves  a  large  amount  of  calculation  (Grena, 
2008),  requiring  the  use  of  tables  to  calculate  heliocentric  latitude,  longitude  and  radius 
vector  on  each  Lagrangian  point,  and  also  “breaks”  some  standards  to  accommodate  for 
solar  radiation  applications.  For  instance,  the  algorithm  estimates  the  azimuth  angle 
measuring  eastward  from  north  instead  of  westward  from  south;  also  it  computes  the 
observer’s  geographical  longitude  considering  negative  west,  or  positive  east  from 
Greenwich,  contrary  to  current  astronomical  conventions  (Reda  and  Andreas,  2003:1). 

Recently,  Roberto  Grena,  a  researcher  from  the  Centro  Ricerche  Casaccia  in 
Rome,  suggested  an  algorithm  that  allows  computing  the  sun  position  very  accurately, 
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with  a  maximal  error  of  0.0027°.  In  his  work,  the  researcher  considered  the  main  effects 
can  affect  the  solar  position  like  the  moon  perturbation,  nutation,  difference  between 
topocentric  and  geocentric  coordinates,  refraction  disturbance  effect,  and  also  introduced 
some  empirical  corrections  in  the  heliocentric  longitude  calculation  to  sum  up  all  the 
other  small  perturbations  too  complex  to  be  considered  one  by  one  (Grena,  2008).  His 
method  is  quite  adapted  to  the  period  2003-2023,  and  strongly  reducing  the  amount  of 
calculations  needed,  especially  the  number  of  trigonometric  functions,  whether  compared 
to  the  previous  mathematical  procedures. 

Another  alternative  option  as  a  step-by-step  problem-solving  procedure  to  track 
accurately  the  solar  position  with  small  computational  effort  is  to  employ  a  Carruthers- 
Szokolay  composite  algorithm.  The  azimuth  and  altitude  of  the  sun  were  calculated  using 
formulae  first  proposed  by  Spencer  (1965),  then  refined  by  Szokolay  (1996).  The  value 
for  solar  declination  was  determined  using  formulae  proposed  by  David  Carruthers 
(1990).  All  angles  are  in  radians  which  is  a  common  requirement  of  most 
implementations  of  trigonometric  functions  in  spreadsheets. 

Although  the  time  of  calculation  depends  on  the  machine  and  on  the  programming 
language  used,  an  approximate  comparison  can  be  found  by  simply  comparing  the 
number  of  times  each  algorithm  calls  a  trigonometric  function,  both  direct  and  inverse, 
since  such  functions  usually  have  a  computational  cost  much  bigger  than  that  of  all  the 
other  operations  (Grena,  2008). 

A  synopsis  of  solar  calculation  methods  are  presented  in  Table  2.  With  respect  to 
accuracy  and  timely  answer  attributes,  the  Reda  and  Andreas  algorithm  was  chosen  for 
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both  its  acceptable  computational  tractability  in  the  VBA  Excel  environment,  and  its 
more  than  satisfactory  precision  in  providing  solar  location  in  the  sky. 


Table  2.  Comparison  of  the  Computational  Cost  Required  by  the  Algorithms 

(Grena,  2008) 


Algorithm 

Additions 

and 

subtractions 

Products 

and 

divisions 

Direct 

trigonometric 

functions 

Inverse 

trigonometric 

functions 

Blanco-Muriel  and  others 

35 

36 

16 

4 

Michalsky 

28 

25 

19 

4 

Grena 

52 

54 

21 

4 

Carruthers  and  others 

26 

35 

16 

2 

Reda  and  Andreas 

>1000 

>1300 

>300 

7 

Manual  Risk  and  Vulnerability  Calculation  Model 

In  non-technical  contexts,  the  word  “risk”  refers  to  situations  in  which  it’s 
possible  but  not  certain  that  some  undesirable  event  will  occur.  In  technical  contexts,  on 
the  other  hand,  this  term  presents  several  more  differentiated  meanings,  as  it  can  be  sawn 
below  (Hansson,  2002:1): 

a)  An  unwanted  event  which  may  or  may  not  occur; 

b)  The  cause  of  an  unwanted  event  which  may  or  may  not  occur; 

c)  The  probability  of  an  unwanted  event  which  may  or  may  not  occur; 

d)  The  statistical  expectation  value  of  an  unwanted  event  which  may  or  may  not  occur; 

e)  The  fact  that  a  decision  is  made  under  conditions  of  known  probabilities. 

Although  most  of  the  above-mentioned  meanings  of  “risk”  have  been  found  in  the 
philosophical  literature,  the  more  informal  definition,  namely  as  “the  probability  of  an 
undesirable  event  may  or  may  not  occur”,  was  chosen  to  be  employed  in  this  research, 
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since  it’s  consistent  with  most  of  the  current  English  lexis.  The  Cambridge  International 
Dictionary  of  English,  for  example,  defines  risk  as  “the  possibility  of  something  bad 
happening”  (Procter,  2001),  while  the  Compact  Oxford  English  Dictionary  opted  for  “the 
possibility  that  something  unpleasant  will  happen”  (Soanes,  2008).  Following  the  same 
sense,  the  United  States  Department  of  Defense  defines  risk  as  “the  probability  and 
severity  of  loss  linked  to  hazards”  (Camm  and  Greefield,  2005:1 1). 

Considering  that  the  probabilistic  approach,  i.e.  existence  of  random  outcomes, 
translates  the  essence  of  risk  idea,  its  logical-mathematical  formalization  presents  itself  as 
the  better  phenomenon  description,  well-expressed  in  the  following  statement. 

“Let  A  be  an  event  belonging  to  the  field  of  event  A  of  an  experiment.  If  event  A 
occurred  nA  times  while  we  repeated  the  experiments  n  times,  then  nA  is  called  the 
frequency,  and  ^  =  hA  is  called  the  relative  frequency  of  the  event  A.  The  relative 

frequency  satisfies  certain  properties  which  can  be  used  to  built  up  an  axiomatic 
definition  of  the  notion  of  the  probability  P(A )  of  event  A  in  the  field  of  events  A.  A  real 
function  P  defined  on  the  field  of  events  is  called  a  probability,  if  it  satisfies  the 
following  properties”  (Bronshtein  and  others,  2007:748-749). 

a)  For  every  event  .4  e  4  we  have: 

0  <  PC4)  <  1  and  0  <  hA  <  1  (1) 

b)  For  the  impossible  event  O  and  the  certain  event  I,  we  have: 

P(0)  =  0,  P<j)  =  1,  and  h0  =  0,  h,  =  i  (2) 

c)  If  the  events  .4;  e  A  (/  =  1,  2,  3,  ...)  are  finite  or  countably  many  mutually  exclusive 
events  (.4,.4;;  =  0  for  c  =  fc),  then: 
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P(Aj  +A2  +  A3  +  -)  =  P(A1)  +  iJU2)-l-P(A3),..Jand 


kdi+da+dj...  —  ^Al  +  ^Al  +  ftflB  +  ■■■ 


(3) 


Therefore,  based  on  the  above-mentioned  theory  of  probability  axioms,  likewise 
the  theoretical  construction  derived  from  literature  review,  “risk  level”  calculation  is 
found  on  the  relationship  between  the  boat  width  (“simple  outcome”)  and  the 
perpendicularly  projected  shadow  length  to  the  path  azimuth  (“field  of  events”),  which 
represents  the  total  linear  shadow  capable  to  cover  the  move  during  travel  time  on  each 
path  (figure  34).  Formally: 


Risk  Level  = 


Boat  width 

Perpendicularly 'projected  shadow  lenght  to  the  path  azimuth 


(4) 


Azimuth,  in  the  military  and  land-navigation  contexts,  can  be  defined  as  “the 
horizontal  angle  of  the  observer's  bearing  in  surveying,  measured  clockwise  from  a 
referent  direction,  as  from  the  north,  or  from  a  referent  celestial  body,  usually  Polaris” 
(Harcourt,  2006).  It  can  be  computed  using  the  following  equation: 


Aztj  —  Ac£J.S 


/sln(gJ-)  —  sin(0j)  X  cos(5i;)> 


sin(5jj)  X  cos(^) 

If  sin(p,  —  Pj)  <  0,  then  Path  Azimuth  —  AzLj 
If  sin(p,  —  p.)  >  0,  then  Path  Azimuth  =  2rr  —  Aztj 
Where: 


(5) 


a)  AztJ :  Path  Azimuth. 

b)  d-i :  latitude  of  start  point, 

c)  8j :  latitude  of  end  point. 
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d)  pi  dongitude  of  start  point, 

e)  p,-: longitude-  of  endpoint 

f) 5i;:  distance  between  start  and  end  points. 


Figure  34.  The  major  attributes  to  risk  level  calculation  on  each  path. 


Combining  solar  altitude,  solar  azimuth  and  path  (river)  azimuth  equations  it’s 
possible  to  define  the  risk  level  calculation. 


risk  level  =  Tilf 


9 

sm(3i?-)  X  toti 
tan 


Where: 

a)  (p:  boat  width. 

b)  5ij-:  difference  between  solar  azimuth  (0)  and  path  azimuth  Oj), 

c)  £0^-:  Riverine  forest  height  between  points  i  and  j. 

d)  Solar  Altitude  between  points  i  and  j,  at  time  t. 


(6) 


It’s  important  to  emphasize  that  StJ-  value  depend  on  the  quarter-circle  in  which 
each  azimuth  is  found,  since  the  risk  is  established  in  terms  of  the  perpendicular  projected 
shadow. 
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a)  Both  Solar  and  Path  Azimuths  are  in  first  quadrant: 


Figure  35:  If  >  rj  then  =  <p  —  7]  (7) 


b)  Path  Azimuth  is  in  the  first  quadrant  and  Solar  Azimuth  is  in  the  second  one: 


Figure  37:  If  (0  —  7)  <  ?j  then 
8tj=  <t>-V 


(9) 


Figure  38:  If  (4>  —  7)  >  V  then 

=  n  —  <p  +  7]  (10) 
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c)  Path  Azimuth  is  in  the  first  quadrant  and  Solar  Azimuth  is  in  the  third  one: 


Figure  39:  If  (0  —n)>r}  then 

$ij  =  4>~n-rj  (11) 


Figure  40:  If  (0  —  n)  <  rj  then 

dij  =  4>  -  (rj  -  n)  (12) 


d)  Path  Azimuth  is  in  the  first  quadrant  and  Solar  Azimuth  is  in  the  fourth  one: 


3tt 

Figure  41:  If  ~ )<  V  then 

5ij=  4>-n~  V 


(13) 


Figure  42:  If  (0 - — )  >  v  then 

Stj  —  2  k  —  <p  +7]  (14) 
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e)  Solar  Azimuth  is  in  the  first  quadrant  and  Path  Azimuth  is  in  the  second  one: 


Figure  43:  If  (?J  —  j)<  <p  then  Figure  44:  If  (Vj  —  >  <p>  then 

SiS=n-rj  +  <l>  (15)  5i}-  =  rj  ~  <t>  (16) 

f)  Solar  Azimuth  is  in  the  first  quadrant  and  Path  Azimuth  is  in  the  third  one: 


Figure  45:  If  (17  —  n)  >  <p  then 

Su=V-4>  (17) 


Figure  46:  If  (tj  —  tt)  <<p  then 

3i}-  =  <P~  0?-O  (18) 
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g)  Both  Solar  Azimuth  and  Path  Azimuth  are  in  the  second  quadrant: 


Figure  47:  If  17  >  <p  then  Figure  48:  If  f]  <  <p  then 

=?j-0  (19)  =  <p-?}  (20) 


h)  Solar  Azimuth  is  in  second  quadrant  and  Path  Azimuth  is  in  the  third  one: 
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i)  Both  Solar  and  Path  Azimuths  are  in  the  third  quadrant: 


Figure  51:  If  is  greater  the  0  then 


Si2-  =  ri  -  0 


n 'W 


Figure  52:  If  rj  is  less  than  0  then 
8ij=  <P~V 


(24) 


j)  Path  Azimuths  is  in  the  third  quadrant  and  Solar  Azimuth  is  in  the  fourth  one: 


Figure  53:  If  (<p  —  — )  <  (rj  —  n)  then 
5ij  =  4>~V  (25) 


Figure  54:  If  {<$  — — )  >  (77  —  tt)  then 
Sij=  n  +  i]-  4>  (26) 
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k)  Both  Path  and  Solar  Azimuths  are  in  the  fourth  quadrant: 


Figure  55:  If  rj>  <p  then  Figure  56:  If  rj  <  <p  then 

Bit=  V-4>  (27)  8if=  <p-v  (28) 

Beyond  the  above  considerations,  some  assumptions  related  to  shadow  making 
are  required: 

a)  The  enemy  searcher  has  an  unlimited  omnidirectional  sensing,  implying  no  random 
search  model.  As  a  result,  this  research  intents  to  study  and  analyze  the  worst  possible 
tactical  and  operational  situation  to  combat  logistics  planners.  The  choice  for  that 
approach  is  due  to  the  critical  features  of  a  search  problem  concerned  with  the  nature 
of  the  target  in  space  and  time  and  the  method  and  efficiency  of  the  search; 

b)  Point  light  sources  do  not  exist  in  practice,  and  hard  shadows  give  a  rather  unrealistic 
feeling  to  images.  Even  the  sun,  probably  the  most  common  shadow-creating  light 
source  in  our  daily  life  has  a  significant  angular  extent  and  does  not  create  hard 
shadows.  In  the  more  realistic  case  of  a  light  source  with  finite  extent,  a  point  on  the 
receiver  can  have  a  partial  view  of  the  light,  i.e.  only  a  fraction  of  the  light  source  is 
visible  from  that  point; 
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c)  The  path  surroundings  consist  laterally  of  the  light  obstacles  (i.e.  riverine  vegetation) 
that  can  occlude  the  move  from  the  observer's  sensors; 

d)  The  distance  between  the  water  surface  and  riverine  forest  ground  is  minor  and 
insignificant  in  the  final  calculation  outcome; 

e)  Shades  created  by  sunlight  are  hard  shadows,  that  is  to  say  they  are  produced  by  a 
point  light  source,  rather  than  soft  ones.  Therefore,  the  shadows  across  the  river  are 
made  only  to  shade,  and  appear  entirely  black.  As  a  result,  the  move  will  travel 
completely  or  partially  hidden,  depending  on  what  percentage  of  the  boat  is  in  the 
shadow. 

Survivability  Alternative  Approaches 

The  risk  between  each  network  node  can  be  easily  computed  using  formula  (6). 
However,  the  total  travel  risk  doesn’t  necessarily  mean  the  sum  of  the  partial  ones.  Since 
risk  was  defined  as  a  state  of  uncertainty  for  which  the  probability  distribution  is  known, 
computing  the  total  travel  risk  depends  on  how  much  each  individual  move  affects  the 
total  mobility  performance  in  terms  of  exposure  level.  Therefore,  the  approach  in  which 
these  probabilistic  outcomes  are  counted  will  define  the  travel  risk  through  the  total  route. 

In  this  research  are  suggested  four  different  mathematical-based  perspectives 
whose  applicability  is  dependent  on  the  camp  commander’s  point  of  view,  as  well  as  his 
logistics  staff,  in  terms  of  mission  criticality. 

The  first  alternative,  called  the  fastest  travel,  computes  the  total  travel  risk  taken 
into  account  the  shortest  path  between  depot  and  deployed  combat  unit,  regardless  of 
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individual  risk  paths.  The  goal  is  to  assess  the  total  exposure  level  considering  the  move 
needs  to  reach  the  destination  as  soon  as  possible. 

A  second  and  more  conservative  approach  is  finding  a  route  with  zero-exposure  to 
enemy  fire.  This  point  of  view  is  elementary  since  it  is  based  on  the  sum  of  individual 
risks.  However,  it  neglects  the  exposure  time  on  each  path,  as  well  as  different  route 
combinations.  It’s  a  relative  measure  instead  of  an  absolute  one.  Consequently,  it  is 
impossible  to  compare  two  or  more  options  in  which  the  number  of  paths  transversed  on 
each  alternative  route  is  different. 

For  example,  if  in  a  hypothetical  movement,  from  point  A  to  B,  an  option  1  has 
individual  path  risks  of  (0,  0.5,  0.5,  0,  0.5)  and  option  2  has  risk  of  (0.5,  0.5,  0.5),  the  two 
options  have  the  same  total  risk  (1.5).  However,  the  exposure  to  enemy  fire  is  at  different 
times,  which,  by  inference,  suggests  that  option  2  can  be  more  dangerous  than  option  1. 
Also,  these  sums  can  result  in  values  greater  than  1,  which  annuls  its  probabilistic 
meaning.  Mathematically  it  can  be  defined  as  the  following: 


T otal  risk  —  A  — 


(29) 


Where: 

a)  Ait:  risk  on  arc  i  at  time  t; 

b)  n :  total  number  of  arcs  on  path; 

The  third  perspective  is  found  on  the  rudimentary  principles  of  serial  reliability  of 
components  in  military  operations  and  systems  (Przemieniecki,  2000:231-232;  Zehna, 
2005:305-309).  In  this  approach,  the  reliability  of  a  military  operation  is  measured  in 
terms  of  its  probability  of  success.  Thus,  the  probability  of  success  of  a  series  of 
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components  comprising  the  overall  mission  is  the  probability  of  success  of  independent 
events.  Since  all  components  in  that  system  operate  in  order  for  the  whole  operation  to  be 
successful.  The  problem  with  this  method  is  related  to  difficulty  in  measuring  the  effect 
of  each  leg  on  the  overall  mission  fulfillment.  For  example,  if  the  movement  is 
completely  vulnerable,  i.e.  unhidden,  in  a  unique  course,  the  military  operation  reliability 
will  be  zero.  However,  this  tool  becomes  important  when  it  is  necessary  to  identify  the 
options  in  which  there  is  a  great  possibility  that  the  boat  will  be  sunk.  The  route 
survivability  is  computed  using  the  following  formula: 


n 


Survivability  =  fi  =  (1  —  Ait) 


(30) 


Where: 

a)  jiit:  risk  on  arc  i  at  time  t; 

b)  n:  total  number  of  arcs  on  path; 

The  fourth  option  is  a  time-dependent  survivability  measure,  which  emphasizes 
the  importance  of  the  period  of  time  that  the  move  is  completely  or  partially  disclosure. 
Under  this  perspective,  it  is  possible  to  compare  different  scenarios,  regardless  of  start 
times,  total  number  of  paths  taken,  and  total  travel  time.  The  weighted-survivability  is 
calculated  employing  the  equation  below. 


weighted  —  survivability  =  = 


(31) 


Where: 


a)  git-  survivability  on  arc  i  at  time  t; 
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b)  n:  total  number  of  arcs  on  path; 

c)  ri;.:  travel  time  between  nodes  i  and  j; 

d)  T :  total  travel  time. 

It  is  important  to  highlighted  that  there  is  no  an  absolute  “best”  or  “superior” 
option  to  calculate  either  total  travel  risk  or  total  travel  vulnerability  since  all  models  are 
found  on  probabilistic  analysis  (Leedy  and  Ormrod,  2009:259-260).  The  alternative 
methods,  as  descriptive  statistics,  are  measures  of  the  moves  behavioral  tendency.  These 
methods  are  based  on  problem  assumptions  and  are  subject  to  interpretation.  A 
prespecified  focus  or  goal  must  to  be  defined  in  advance  in  order  to  choose  the  most 
appropriate  tool  for  analyzing  the  problem. 

Finally,  beyond  the  above  considerations,  some  assumptions  related  to  risk  as 
well  as  survivability  evaluation  are  required: 

a)  Risk  between  any  two  nodes  is  calculated  based  on  the  second  node’s  data,  in  terms  of 
its  geographical  coordinates,  and  solar  position  data  (sun  azimuth  and  sun  altitude); 

b)  Since  risk  is  dependent  of  the  trinomial  move/location/time,  it  is  in  constant  change. 
However,  to  simplify  the  calculations  risk  is  considered  constant  on  each  path. 

Graph  Algorithms 

A  graph  G  =  (V,  E)  is  defined  by  a  set  of  vertices  V,  and  a  set  of  edges  E 
consisting  of  ordered  or  unordered  pairs  of  vertices  from  V.  In  a  network,  the  vertices 
may  represent  the  nodes  or  junctions.  Certain  pairs  of  which  are  connected  by  edges, 
which  can  be  wires,  roads,  rivers,  pipes,  and  so  on  (Papadimitriou  and  Steiglitz,  1998:20). 
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The  basic  operation  in  most  graph  algorithms  is  completely  and  systematically 
traversing  the  graph,  visiting  every  vertex  and  every  edge  exactly  once  in  some  well- 
defined  order.  There  are  two  primary  kinds  of  traversal  algorithms:  Breadth-first  Search 
(BFS)  and  Depth-first  Search  (DFS)  (Castillo  and  others,  1996:160;  Gross  and  Yellen, 
1998:130).  For  certain  problems,  it  makes  absolutely  no  difference  which  one  you  use, 
but  in  other  cases  the  distinction  is  crucial.  Both  graph  traversal  procedures  share  one 
fundamental  idea,  namely,  that  it  is  necessary  to  mark  the  vertices  seen  before  so  they 
aren’t  explored  again.  BFS  and  DFS  differ  only  in  the  order  in  which  they  explore 
vertices. 

The  BFS  algorithm  starts  at  a  given  vertex,  which  is  at  level  0.  In  the  first  stage, 
all  vertices  are  visited  at  level  1 .  In  the  second  stage,  all  vertices  of  the  second  level  are 
visited.  These  are  new  vertices,  which  are  adjacent  to  level  1  vertices  (Cormen  and 
others,  2001:531).  The  BFS  traversal  terminates  when  every  vertex  has  been  visited.  It’s 
used  to  solve  the  following  problems: 

a)  Testing  whether  a  graph  is  connected. 

b)  Computing  a  spanning  forest  of  graph. 

c)  Computing,  for  every  vertex  in  graph,  a  path  with  the  minimum  number  of  edges 
between  start  vertex  and  current  vertex  or  reporting  that  no  such  path  exists. 

d)  Computing  a  cycle  in  graph  or  reporting  that  no  such  cycle  exists. 

The  DFS  algorithm  starts  at  a  specific  vertex-source  (S)  in  G,  which  becomes  the 
current  vertex.  The  algorithm  traverses  the  graph  by  any  edge  (Xj,Xj)  incident  to  the 
current  vertex  u.  If  the  edge  (X„  Xj)  leads  to  an  already  visited  vertex  Xj,  then  we 
backtrack  to  current  vertex  Xt.  Backtracking  ensures  correctness  by  enumerating  all 
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possibilities.  It  ensures  efficiency  by  never  visiting  a  state  more  than  once.  If,  on  the  other 
hand,  edge  (Xh  Xj)  leads  to  an  unvisited  vertex  Xj,  then  we  go  to  Xj  and  Xj  becomes  our 
current  vertex.  We  proceed  in  this  manner  until  we  reach  a  “deadend”.  At  this  point  we 
start  back  tracking.  The  process  terminates  when  backtracking  leads  back  to  the  start 
vertex  (Cormen  and  others,  2001:541-542).  Edges  that  lead  to  new  vertex  are  called 
discovery  or  tree  edges  and  edges  that  lead  to  already  visited  edges  are  called  “back 
edges”.  DFS  algorithm  is  used  to  solve  the  following  problems. 

a)  Testing  whether  a  graph  is  connected. 

b)  Computing  a  spanning  forest  of  graph. 

c)  Computing  a  path  between  two  vertices  of  graph  or  equivalently  reporting  that  no  such 
path  exists. 

d)  Computing  a  cycle  in  graph  or  equivalently  reporting  that  no  such  cycle  exists. 

There  are  two  points  to  remember  about  using  BFS  to  find  a  shortest  path  from  Xj 
to  Xj.  First,  the  shortest  path  tree  is  only  useful  if  BFS  was  performed  with  x  as  the  root  of 
the  search.  Second,  BFS  only  gives  the  shortest  path  if  the  graph  is  unweighted. 
Therefore,  either  simple  DFS  or  BFS  is  able  to  solve  the  purposed  problem.  However,  the 
Dijkstra's  single  source  shortest  path  algorithm,  a  well-known  algorithmic  program  used 
in  routing  problems,  which  uses  the  BFS  principle,  is  a  good  option  to  compute  the  total 
travel  risk. 

Dijkstra 's  Algorithm 

The  Dijkstra’s  Algorithm  is  one  of  the  most  celebrated  algorithms  in  the 
Computer  Science  and  very  popular  in  the  Operations  Research  (Avron  and  others, 
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2008:589;  Du  and  others,  2009:187).  One  of  the  main  reasons  for  its  popularity  is  its  use 
in  generating  (exact)  optimal  solutions  to  a  large  class  of  shortest  path  problems  (Bubak 
and  others,  2004:965;  LaValle,  2006:47).  Also,  the  main  point  is  that  the  shortest  path 
problem  is  extremely  important  theoretically  as  well  as  practically,  because  any 
combinatorial  optimization  problem  is  formulated  as  a  shortest  path  problem 
(Papadimitriou  and  Steiglitz  1998:75).  This  class  of  problems  is  extremely  large  and 
includes  numerous  practical  problems  that  have  nothing  to  do  with  actual  (“genuine”) 
shortest  path  problems. 

The  typical  descriptions  of  this  algorithm  start  by  postulating  that  the  shortest 
paths  be  enumerated  in  the  order  of  increasing  distances  from  the  source  (Sack  and 
Urrutia,  2001:636).  Basically,  it  works  by  solving  a  subproblem  k,  which  computes  the 
shortest  path  from  the  source  v  to  vertices  among  the  k  closest  vertices  to  the  source 
(Zang,  2002:17).  However  for  the  Dijkstra’s  algorithm  works  it  should  be  directed- 
weighted  graph  and  the  edges  should  be  non-negative  (Cormen  and  others,  2001:582; 
Manber,  1989:204-206).  If  the  edges  are  negative  then  the  actual  shortest  path  cannot  be 
obtained.  At  the  kth  round,  there  will  be  a  set  called  frontier  of  k  vertices  that  will  consist 
of  the  vertices  closest  to  the  source  and  the  vertices  that  lie  outside  frontier  are  computed 
and  put  into  new  frontier.  The  shortest  distance  obtained  is  stored  in  w.SP.  It  holds  the 
estimate  of  the  distance  from  v  to  w.  Dijkstra’s  algorithm  finds  the  next  closest  vertex  by 
maintaining  the  new  frontier  vertices  in  a  priority-min  queue  (Cormen  and  others, 
2001:599). 

The  algorithm  works  by  keeping  the  shortest  distance  of  vertex  w  from  the  source 
v  in  an  array,  w.SP.  The  shortest  distance  of  the  source  to  itself  is  zero.  The  w.SP  for  all 
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other  vertices  is  set  to  infinity  to  indicate  that  those  vertices  are  not  yet  processed.  After 
the  algorithm  finishes  the  processing  of  the  vertices,  w.SP  will  have  the  shortest  distance 
of  vertex  v  to  w.  Two  sets  are  maintained,  frontier  and  new  frontier,  which  helps  in  the 
processing  of  the  algorithm.  Frontier  has  k  vertices  which  are  closest  to  the  source,  will 
have  already  computed  shortest  distances  to  these  vertices,  for  paths  restricted  up  to  k 
vertices.  The  vertices  that  reside  outside  of  frontier  are  put  in  a  set  called  new  frontier. 
The  following  lines  present  the  Dijkstra’s  algorithm’s  pseudocode  (Manber,  1989:206) 

Input:  G  =  (V,  E)  (a  weighted  directed  graph),  and  v  (the  source  vertex) 

Output:  for  each  vertex  w,  w.SP  is  the  length  of  the  shortest  path  from  v  to  w. 

{all  lengths  are  assumed  to  be  nonnegative}. 

begin 

for  all  vertices  w  do 
w.mark  :=  false  ; 

v. SP  :=  co ; 

v.SP  :=  0; 

while  there  exist  an  unmarked  vertex  do 

let  w  be  an  unmarked  vertex  such  w.SC  is  minimal ; 

w. mark  ■=  true ; 

for  all  edges  (w,v)  such  that  z  is  unmarked  do 
if  w.  SP  +  lenght(w,z)  <  z.  SP  then 
z.SP  :■=  w. SP  +  length  (w, z) 

end 

In  this  research,  as  previously  mentioned,  the  path  weights  are  based  on  travel  risk 
instead  of  travel  cost.  Dijkstra’s  algorithm  is  still  appropriate  since  the  graph  is  direct- 
weighted  and  edges  are  always  non-negative. 
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IV.  Results  and  Analysis 


Introduction 

In  order  to  determine  the  absolute  performance  of  the  designed  model 
configuration  under  the  assumptions  previously  set  up  both  verification  and  validation 
processes  were  performed  taking  into  account  actual  non-classified  information  from  the 
Amazonian  region  in  which  the  scenario  problem  was  projected  (Asner  and  others,  2002; 
McClain  and  others,  2001;  Spencer,  D.,  2004). 

Since  these  procedures  are  concerned  with  accuracy  transformation  and 
representation,  are  essential  prerequisites  to  the  credible  and  reliable  use  of  this  model  as 
well  as  its  results  in  military  planning  processes  (Abu-Taieh  and  others,  2009:58;  Banks 
and  Carson,  1996;  Kheir,  1996:220;  Law,  2007:243-244).  Therefore,  based  on  current 
simulation  literature,  some  methods  were  selected  to  verify  and  validate  the  present 
model. 

For  verification,  the  data  analysis  process  was  employed.  The  model  validation, 
on  the  other  hand,  was  obtained  using  cause-effect  graphing  method,  extreme  conditions 
test,  and  predictive  validation  dynamic  methods  (Abu-Taieh  and  others,  2009:60-65). 

Problem  Model  Representation 

One  of  the  most  difficult  problems  facing  a  simulation  analysis  is  trying  to 
determine  whether  a  model  is  an  accurate  representation  of  the  system  being  studied 
(Kent  and  Williams,  1990:232;  Law,  2007:243;  Lenhard  and  others,  2006:150).  Based  on 
this  assertion,  real  data  about  tropical  and  subtropical  moist  broadleaf  forest  environment, 
particularly  from  the  Brazilian  Amazon  basin,  was  collected  before  starting  any  kind  of 
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model  assessment  (figure  57).  This  data  was  translated  into  a  theoretical  network,  which 
constitutes  graphic  representation  of  a  hypothetical  theater  of  operations  in  which 
military  supply  missions  will  be  performed  (figures  58  and  59). 


Figure  57.  Theoretical  network  plot. 


Figure  58.  Amazon  region  satellite  picture. 
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Figure  59.  Theoretical  network  based  on  real  watercourse  system. 


The  theoretical  network  is  composed  of  64  distinct  nodes  with  arcs  connections 
the  nodes  as  shown  in  the  figure  57.  The  nodes  are  identified  by  its  geographical 
coordinates  (Appendix  A).  These  points  represent  all  possible  locations  for  both  a  small 
deployed  combat  force  and  a  military  depot  station  responsible  for  supporting  the  combat 
force. 

Verification  Process 

The  verification  process  refers  to  a  testing  procedure  that  determines  whether  a 
system  is  consistent  with  its  designations  as  well  as  both  problem  assumptions  and 
circumstantial  issues  related  to  the  question  have  been  correctly  translated  into  a 
computer  program  (Law,  2007:243).  The  typical  questions  to  be  answered  during  this 
process  include  (Sokolowiski  and  Banks,  2009:126): 

a)  Does  the  program  code  of  the  executable  model  correctly  implement  the  conceptual 
model? 

b)  Does  the  conceptual  model  satisfy  the  intended  uses  of  the  model? 
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c)  Does  the  executable  model  produce  results  when  it’s  needed  and  in  the  required 


format? 

To  answer  these  questions  the  following  experiment  was  performed: 
a)  The  logical-mathematical  consistency  between  the  conceptual  model  and  coding 
system  were  evaluated  comparing  the  set  of  functions  and  variables  with  its 
counterparts  as  well  as  the  theoretical  constructs  that  represents  the  phenomena  in 
study,  like  physical  laws,  hypothesis  and  related  theories. 

Data  Analysis 

This  verification  tool  compares  data  definitions  and  operations  in  the  conceptual 
model  to  those  in  the  executable  one.  It  includes  data  dependency  analysis,  which 
determines  which  variable  depends  on  the  other  one,  and  data  flow  analysis.  This  is 
useful  for  undefined  and  unused  data  and  is  also  helpful  in  identifying  inconsistencies  in 
data  structure  and  improper  linkages  among  data  items  (Ledin,  2001:209-210). 

Though  all  computation  are  performed  in  the  VBA  environment,  it’s  possible  to 
attest  model  congruence  in  terms  of  parameter  dependences,  by  bringing  functions  and 
variables  employed  in  the  coding  process  from  VBA  to  an  Excel  Spreadsheet  and 
examining  their  connections.  For  example,  as  mentioned  on  chapter  III,  risk  level 
depends  directly  on  the  pair  sun-move  position,  since  their  relative  location  generates  (or 
not)  shadow  areas.  That  internal  relationship  can  be  identified  comparing  the  content  of 
the  red  rectangles  (figures  60  and  61),  with  the  real-world  source  (figures  62  and  63),  as 
well  as  the  changes  due  to  alteration  in  one  argument,  particularly  in  the  timeframe  (from 
8:00am  to  7:00am). 
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Figure  60.  Excel  Spreadsheet  screenshot  showing  the  relationship  among  the  solar 
position  parameters  and  timeframe  (on  March  21th  2010,  at  8:00am). 
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Figure  61.  Excel  Spreadsheet  screenshot  showing  the  relationship  among  the  solar 
position  parameters  and  timeframe  (on  March  21th  2010,  at  7:00am). 


It  is  possible  to  note  on  the  above  pictures  that  before  the  sunrise  time  (over 
7.2183  hours,  i.e.  7h  13min)  solar  altitude  has  a  negative  value,  which  can  be  confirmed 
on  the  NOAA  Solar  Calculator  website  (figures  62  and  63). 
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Figure  62.  NOAA  Solar  Calculation  screenshot  showing  the  relationship  among  the  solar 
position  parameters  and  timeframe  (on  March  21th  2010,  at  8:00am). 
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Figure  63.  NOAA  Solar  Calculation  screenshot  showing  the  relationship  among  the  solar 
position  parameters  and  timeframe  (on  March  21th  2010,  at  7:00am). 
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Since  path  azimuth  is  essential  data  to  risk  computation,  its  estimation  is 
evaluated  in  the  VBA  environment  (figure  64)  compared  against  the  values  provided  by 
distance  and  azimuth  computation,  found  on  the  Federal  Communication  Commission 
(FCC)  website  (figure  65).  In  this  condition,  the  path  azimuth  will  be  360°  -  280.368°, 
which  equals  79.63°. 


|  (General) 
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Figure  64.  Excel  VBA  screenshot  showing  the  path  azimuth  calculation, 
between  nodes  1  and  2,  computed  by  program  computer. 
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Figure  65.  FCC  Distance  and  Azimuth  Calculation  showing  the  path 
azimuth  between  nodes  1  and  2. 
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Using  data  analysis  techniques  it  was  possible  to  verify  the  computer  program 
against  conceptual  model  comparing  functions  and  variables  that  represents  the 
phenomena  in  study.  This  assesses  as  a  result,  its  logical-mathematical  consistency. 

Validation  Process 

According  to  current  simulation  literature,  validation  process  refers  to  a  testing 
procedure  that  determines  whether  a  system  satisfies  the  user  requirements  by  comparing 
simulation  results  derived  from  the  model  against  the  real-world  system.  Also,  it’s 
understood  as  a  “process  of  determining  the  degree  to  which  the  model  is  an  accurate 
representation  of  the  real  world”  (DoD,  1996).  “Degree”  is  a  subjective  term,  so  a  model 
is  considered  valid  when  it’s  a  reasonable  representation  of  the  real  system. 

In  general,  the  validation  of  models  of  physical  phenomena  is  straightforward 
since  the  laws  of  nature  are  usually  well  known  and  mathematically  precise.  In  this 
specific  case,  the  model  is  validated  by  comparing  the  results  of  running  the  model  with 
observations  from  the  real  world  (Zacharias  and  others,  2008:343).  However,  in  order  to 
standardize  its  process,  some  typical  and  standard  questions  need  to  be  answered,  which 
include  (Sokolowiski  and  Banks,  2009:126): 

a)  Is  the  conceptual  model  a  correct  representation  of  the  simulation? 

b)  How  close  are  the  results  produced  by  the  executable  model  to  the  behavior  of  the  real 
world? 

c)  Under  what  range  of  inputs  are  the  model’s  results  credible  and  useful? 

To  achieve  these  objectives  the  following  tests  were  performed: 
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a)  The  model  was  evaluated  on  the  extremes  scenarios,  that  is,  under  the  non-existence  of 
riverine  forest  as  well  as  in  the  presence  of  extremely  high  trees  surrounding  the 
routes; 

b)  The  boat  speed  was  progressively  changed  from  a  small  value  to  a  large  one; 

c)  The  boat  width  was  progressively  changed  from  a  small  value  to  a  large  one. 

Extreme  Condition  Tests 

This  first  procedure  consists  in  selecting  utmost  numerical  inputs,  in  which  the 
model  outcome  is  easily  predicable,  considering  its  senselessness  and  improbable 
combination  of  levels  in  parts  of  the  model.  This  is  used  to  check  whether  the  model 
responds  accordantly,  in  other  words,  presenting  also  extremes  results  (Abu-Taieh  and 
others,  2009:63).  The  two  more  significant  instances  in  this  study  are  related  to  the  total 
absence  of  vegetation  along  the  river,  which  doesn’t  produce  shadow  at  all,  and 
incredibly  high  trees,  the  “perfect  shadow  makers”. 

The  figures  66,  67  and  68  present  the  input  and  outcome  screenshots,  showing 
what  information  was  placed  as  well  as  the  respective  results.  In  these  pictures,  it  is 
evident  the  complete  lack  of  trees  doesn’t  produce  shadow  areas  and,  as  result,  the 
relative  risk  level  decreases  over  time,  depending  on  the  move  start  time.  At  the  same 
time,  the  probability  of  mission  success,  i.e.  survivability,  is  zero.  The  four  subsequent 
figures  (69  to  72),  graphically  display  the  model  behavior,  exactly  as  it  was  expected. 
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Figure  66.  Input  Spreadsheet’s  Excel  screenshot,  considering 
30-minute  interval  departure. 


Figure  67.  Output  Spreadsheet’s  Excel  screenshot,  considering 
30-minute  interval  departure. 
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Figure  68.  Output  Spreadsheet’s  Excel  screenshot,  considering 
15 -minute  interval  departure. 
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Figure  69.  Graphic  presenting  the  boat 
capability  of  surviving  the  supply  mis¬ 
sion  under  the  absence  of  tree  through¬ 
out  the  complete  route,  with  30-minute 
interval  departure.  Regardless  start  time, 
if  one  path  has  risk  1,  the  probability  of 
mission  success  will  be  zero. 


Weighted  Survivability  Analysis 
30-minute  Interval  Departure 


Figure  70.  Graphic  presenting  the  boat 
capability  of  surviving  the  supply  mis¬ 
sion  under  tree  absence  throughout  the 
complete  route,  emphasizing  the  impor¬ 
tance  of  travel  time,  with  30-minute  in¬ 
terval  departure. 
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Risk  Analysis 
30-minute  Interval  Departure 


Figure  71.  Graphic  presenting  the  relative 
travel  risk  throughout  the  complete  route, 
with  30-minute  interval  departure. 
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Figure  72.  Graphic  presenting  the  boat 
capability  of  surviving  the  supply  mis¬ 
sion  under  tree  absence  throughout  the 
complete  route,  emphasizing  the  impor¬ 
tance  of  travel  time,  with  15 -minute  in¬ 
terval  departure. 


Figures  73,  74  and  75  present  the  input  and  outcome  screenshots  showing  what 
data  was  placed  as  well  as  the  result  model.  The  existence  of  higher  riverine  arboreal 
biome  produces  larger  shadow  areas  and,  as  a  result,  the  risk  level  is  null  on  all  paths. 
The  probability  of  mission  success,  is  100%.  The  four  subsequent  figures  (76  to  79), 
graphically  display  the  model  behavior,  exactly  as  it  was  expected. 
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Figure  73.  Input  Spreadsheet’s  Excel  screenshot, 
considering  30-minute  interval  departure. 


Figure  74.  Output  Spreadsheet’s  Excel  screenshot, 
considering  30-minute  interval  departure. 
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Figure  75.  Output  Spreadsheet’s  Excel  screenshot, 
considering  15 -minute  interval  departure. 
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Figure  76.  Graphic  presenting  the  boat 
capability  of  surviving  the  supply  mis¬ 
sion  under  higher  tree  path  condition 
throughout  the  complete  route,  with  30- 
minute  interval  departure. 
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Figure  77.  Graphic  presenting  the  boat 
capability  of  surviving  the  supply  mis¬ 
sion  under  higher  tree  path  condition 
throughout  the  complete  route,  empha¬ 
sizing  the  importance  of  travel  time, 
with  30-minute  interval  departure. 
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Figure  78.  Graphic  presenting  the  relative 
travel  risk  under  higher  tree  path  condi¬ 
tion  throughout  the  complete  route,  with 
30-minute  interval  departure. 
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Figure  79.  Graphic  presenting  the  boat 
capability  of  surviving  the  supply  mission 
under  higher  tree  path  condition  through¬ 
out  the  complete  route,  emphasizing  the 
importance  of  travel  time,  with  15 -minute 
interval  departure. 


Cause-Effect  Graphing 

The  second  verification  method  compares  the  cause-and-effect  relationships  in  the 
real  world  to  those  in  the  conceptual  model  (Dasso  and  Funes,  2007:172;  Haug  and 
others,  2001 :47).  Causes  are  events  or  conditions  that  may  occur  in  the  real  world  and  the 
effects  are  the  consequences  or  state  changes  that  result  from  the  causes.  In  this  study,  the 
most  significant  parameters  that  affect  the  enemy  exposure  level  (either  increase  or 
decrease)  are  boat  speed  and  width,  as  well  as  the  actual  riverine  forest  configuration  on 
the  network. 

Changing  boat  speed  gradually  (figures  80  to  84)  from  small  velocities  to  large 
ones,  it  is  possible  to  identify  smaller  variation  on  both  the  travel  risk  and  move 
survivability  (figures  85  to  100).  This  fact  is  due  to,  outside  the  first  hours  of  the  morning 
as  well  as  the  lasts  ones  of  the  afternoon  a  boat’s  need  to  transverse  the  predefined 
network  completely  from  node  1  to  node  38  regardless  of  its  speed.  Since  the  move 
becomes  faster  it  attains  high-level  risk  locations  earlier  than  if  it  was  slower.  Therefore, 
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considering  the  larger  distances  between  each  node,  the  unique  way  to  perform  a  zero- 


risk  travel  is  to  move  extremely  fast  over  the  river  which,  as  a  result,  will  require 
powerful  and  noisy  engine.  This  results  in  unfeasible  conduct  for  a  stealthy  movement. 
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Figure  80.  Input  Spreadsheet’s  Excel  screenshot,  with  rate 
1.86  miles  per  hour. 
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Figure  81.  Output  Spreadsheet’s  Excel  screenshot,  with  rate 
1 .86  miles  per  hour. 
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Figure  82.  Output  Spreadsheet’s  Excel  screenshot,  with  rate 
3.1 1  miles  per  hour. 
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Figure  83.  Output  Spreadsheet’s  Excel  screenshot,  with  rate 
6.21  miles  per  hour. 
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Figure  84.  Output  Spreadsheet’s  Excel  screenshot,  with  rate 
9.32  miles  per  hour. 


Survivability  Analysis 
30-minute  Interval  Departure 


Figure  85.  Graphic  presenting  the  boat  ca¬ 
pability  of  surviving  the  supply  mission, 
given  rate  1.86  miles  per  hour,  with  30- 
minute  interval  departure. 


Weighted  Survivability  Analysis 
30-minute  Interval  Departure 


Figure  86.  Graphic  presenting  the  boat  ca¬ 
pability  of  surviving  the  supply  mission, 
given  rate  1.86  miles  per  hour,  empha¬ 
sizing  the  importance  of  travel  time,  with 
30-minute  interval  departure. 
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Risk  Analysis 
30-minute  Interval  Departure 


Figure  87.  Graphic  presenting  the  relative 
travel  risk,  given  rate  1.86  miles  per  hour, 
with  30-minute  interval  departure. 


Survivability  Analysis 
30-minute  Interval  Departure 


Figure  89.  Graphic  presenting  the  boat  ca¬ 
pability  of  surviving  the  supply  mission, 
given  rate  3.11  miles  per  hour,  with  30- 
minute  interval  departure. 


Weighted  Survivability  Analysis 
15-minute  Interval  Departure 


Figure  88.  Graphic  presenting  the  boat  ca¬ 
pability  of  surviving  the  supply  mission, 
given  rate  1.86  miles  per  hour,  empha¬ 
sizing  the  importance  of  travel  time,  with 
15-minute  interval  departure. 
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Figure  90.  Graphic  presenting  the  boat  ca¬ 
pability  of  surviving  the  supply  mission, 
given  rate  3.11  miles  per  hour,  empha¬ 
sizing  the  importance  of  travel  time,  with 
30-minute  interval  departure. 
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Risk  Analysis 
30-minute  Interval  Departure 


Figure  91.  Graphic  presenting  the  relative 
travel  risk  throughout  the  complete  route, 
given  rate  3.11  miles  per  hour,  with  30- 
minute  interval  departure. 


Weighted  Survivability  Analysis 
15-minute  Interval  Departure 


Figure  92.  Graphic  presenting  the  boat  ca¬ 
pability  of  surviving  the  supply  mission, 
given  rate  3.11  miles  per  hour,  empha¬ 
sizing  the  importance  of  travel  time,  with 
15 -minute  interval  departure. 
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Figure  93.  Graphic  presenting  the  boat  ca¬ 
pability  of  surviving  the  supply  mission, 
given  rate  6.21  miles  per  hour,  with  30- 
minute  interval  departure. 
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Figure  94.  Graphic  presenting  the  boat  ca¬ 
pability  of  surviving  the  supply  mission, 
given  rate  6.21  miles  per  hour,  empha¬ 
sizing  the  importance  of  travel  time,  with 
30-minute  interval  departure. 
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Risk  Analysis 
30-minute  Interval  Departure 


Figure  95.  Graphic  presenting  the  relative 
travel  risk  throughout  the  complete  route, 
given  rate  6.21  miles  per  hour,  with  30- 
minute  interval  departure. 
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Figure  97.  Graphic  presenting  the  boat  ca¬ 
pability  of  surviving  the  supply  mission, 
given  rate  9.32  miles  per  hour,  with  30- 
minute  interval  departure. 


Weighted  Survivability  Analysis 
15- minute  Interval  Departure 


Figure  96.  Graphic  presenting  the  boat  ca¬ 
pability  of  surviving  the  supply  mission, 
given  rate  6.21  miles  per  hour,  empha¬ 
sizing  the  importance  of  travel  time,  with 
15 -minute  interval  departure. 
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Figure  98.  Graphic  presenting  the  boat  ca¬ 
pability  of  surviving  the  supply  mission, 
given  rate  9.32  miles  per  hour,  empha¬ 
sizing  the  importance  of  travel  time,  with 
30-minute  interval  departure. 
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Risk  Analysis 
30-minute  Interval  Departure 
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Figure  99.  Graphic  presenting  the  relative 
travel  risk  throughout  the  complete  route, 
given  rate  9.32  miles  per  hour,  with  30- 
minute  interval  departure. 
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Figure  100.  Graphic  presenting  the  boat 
capability  of  surviving  the  supply  mission, 
given  rate  9.32  miles  per  hour,  empha¬ 
sizing  the  importance  of  travel  time,  with 
15 -minute  interval  departure. 


Regardless  of  boat  speed,  any  alteration  on  boat  width  (figure  101)  causes 
significant  variation  on  both  the  travel  risk  and  travel  survivability  (figures  102  to  105). 
Analyzing  the  boat’s  behavior  with  distinct  boat  lengths  highlights  the  increase  in  risk 
because  of  the  enlarged  exposure  area  figures  106  to  121). 
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Figure  101.  Input  Spreadsheet’s  Excel  screenshot, 
considering  a  boat  width  of  1.64  feet. 
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Figure  102.  Output  Spreadsheet’s  Excel  screenshot, 
considering  a  boat  width  of  1.64  feet. 
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Figure  103.  Output  Spreadsheet’s  Excel  screenshot, 
considering  a  boat  width  of  3.28  feet. 
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Figure  104.  Output  Spreadsheet’s  Excel  screenshot, 
considering  a  boat  width  of  4.92  feet. 
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Figure  105.  Output  Spreadsheet’s  Excel  screenshot, 
considering  a  boat  width  of  6.56  feet. 
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Survivability  Analysis 
30-minute  Interval  Departure 


Figure  106.  Graphic  presenting  the  boat 
capability  of  surviving  the  supply  mission, 
given  1.64  feet  of  width,  with  30-minute 
interval  departure. 
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Figure  108.  Graphic  presenting  the  rela¬ 
tive  travel  risk  throughout  the  complete 
route,  given  1.64  feet  of  width,  with  30- 
minute  interval  departure. 


Weighted  Survivability  Analysis 
30-minute  Interval  Departure 


Figure  107.  Graphic  presenting  the  boat 
capability  of  surviving  the  supply  mission, 
given  1.64  feet  of  width,  emphasizing  the 
importance  of  travel  time,  with  30-minute 
interval  departure. 
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Figure  109.  Graphic  presenting  the  boat 
capability  of  surviving  the  supply  mission, 
given  1.64  feet  of  width,  emphasizing  the 
importance  of  travel  time,  with  15 -minute 
interval  departure. 
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Survivability  Analysis 
30-minute  Interval  Departure 


Figure  110.  Graphic  presenting  the  boat 
capability  of  surviving  the  supply  mission, 
given  3.28  feet  of  width,  with  30-minute 
interval  departure. 
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30-minute  Interval  Departure 


Figure  111.  Graphic  presenting  the  boat 
capability  of  surviving  the  supply  mission, 
given  3.28  feet  of  width,  emphasizing  the 
importance  of  travel  time,  with  30-minute 
interval  departure. 
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Figure  112.  Graphic  presenting  the  relative 
travel  risk  throughout  the  complete  route, 
given  3.28  feet  of  width,  with  30-minute 
interval  departure. 
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Figure  113.  Graphic  presenting  the  boat 
capability  of  surviving  the  supply  mission, 
given  3.28  feet  of  width,  emphasizing  the 
importance  of  travel  time,  with  15 -minute 
interval  departure. 
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Survivability  Analysis 
30-minute  Interval  Departure 


Figure  114.  Graphic  presenting  the  boat 
capability  of  surviving  the  supply  mission, 
given  4.92  feet  of  width,  with  30-minute 
interval  departure. 
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Figure  116.  Graphic  presenting  the  relative 
travel  risk  throughout  the  complete  route, 
given  4.92  feet  of  width,  with  30-minute 
interval  departure. 


Weighted  Survivability  Analysis 
30-minute  Interval  Departure 


Figure  115.  Graphic  presenting  the  boat 
capability  of  surviving  the  supply  mission, 
given  4.92  feet  of  width,  emphasizing  the 
importance  of  travel  time,  with  30-minute 
interval  departure. 
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Figure  117.  Graphic  presenting  the  boat 
capability  of  surviving  the  supply  mission, 
given  4.92  feet  of  width,  emphasizing  the 
importance  of  travel  time,  with  15 -minute 
interval  departure. 
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Survivability  Analysis 
30-minute  Interval  Departure 


Figure  118.  Graphic  presenting  the  boat 
capability  of  surviving  the  supply  mis¬ 
sion,  given  6.56  feet  of  width,  with  30- 
minute  interval  departure. 
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Figure  120.  Graphic  presenting  the  rela¬ 
tive  travel  risk  throughout  the  complete 
route,  given  6.56  feet  of  width,  with  30- 
minute  interval  departure. 
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Figure  119.  Graphic  presenting  the  boat 
capability  of  surviving  the  supply  mission, 
given  6.56  feet  of  width,  emphasizing  the 
importance  of  travel  time,  with  30-minute 
interval  departure. 
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Figure  121.  Graphic  presenting  the  boat  ca¬ 
pability  of  surviving  the  supply  mission, 
given  6.56  feet  of  width,  emphasizing  the 
importance  of  travel  time,  with  15 -minute 
interval  departure. 
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In  this  study  the  magnitude  of  the  influence  of  the  variable  “boat  width”  on 
stealthy  boat  ability  became  evident.  Therefore,  the  more  dimensional  analysis, 
conducted  against  risk/survivability  assessment  can  help  military  planners  designing 
stealthier  watercraft  models  that  could  be  employed  in  dense  tropical  forest  streams  in 
order  to  avoid  or  reduce  aerial  detection  by  enemy  searchers.  Also,  analyzing  the  Excel 
output  screenshot,  it  is  possible  to  verify  that  faster  travel,  starting  at  earliest  travel  time, 
rarely  produces  optimal  stealthy  routes,  except  for  the  extreme  conditions  of  either 
excessively  high  trees  or  powerful  boat  engines. 

Predictive  Validation 

This  tool  is  a  dynamic  verification  procedure  used  to  compare  specific  outcomes 
in  the  real  world  to  corresponding  outomes  in  the  model.  It’s  considered  “a  primary 
method  for  providing  evidence  that  a  simulation  is  an  accurate  representation  of  an  actual 
system”  (Ledin,  2001:217).  Given  some  data,  the  model  is  executed  with  the  same  inputs 
and  its  results  are  compared  with  the  historical  or  experimental  data. 

In  this  experiment,  the  model  was  ran,  under  ordinary  conditions  (figure  91),  in 
order  to  assesses  its  result.  After  that,  parameters  directly  related  to  the  provided  outcome 
will  be  replaced  in  order  to  verify  how  the  model  reacts  to  the  new  condition. 

Once  the  model  ran,  the  best  routes  chosen,  i.e.  these  that  provide  null  risk,  were 
the  following  (figure  92): 

a)  1-5-10-18-23-32-33-37-38:  starting  between  8:00h  and  8:45h; 

b)  1-5-10-12-20-25-34-36-48-5 1-57-59-55-47-39-29-38:  starting  between 9:00h  and 9: 15h; 

c)  1-5-10-12-20-25-34-36-48-48-51-49-43-38:  starting  between  9:30h  and  10:30h; 
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d)  1-5-10-18-23-32-38:  starting  between  15:30hand  16:30h; 

e)  1-5-10-12-20-24-33-37-38:  starting  at  16:45h; 

f)  1-2-7-13-16-22-28-29-38:  starting  at  17:00h; 

g)  1-5-10-18-23-32-38:  starting  between  17:15hand  18:00h. 

Analysing  the  results  presented  above,  it’s  evident  that  curretly  defined  route  A 
(starting  from  1-5-10-...)  provides  the  most  stealthy  routes,  compared  againt  all  possible 
solutions  computed  by  Dijkstra's  algorithm.  If  variables  connected  to  these  routes  are 
radically  altered  (figure  94)  the  model  behaves  exactly  as  it  was  expected:  another  route, 
called  B,  arises  from  the  new  scenario  (figure  95),  in  which  the  last  option  doesn’t 
provide  at  all  the  best  alternative  in  terms  of  a  low  exposure  level  to  enemy  observation 
(figure  92). 


A1 

A 

B 

C 

D 

E 

F 

G 

H 

I  J 

K 

L 

M 

N 

O 

P 

Q 

R 

S 

T 

U 

V 

w 

1 

R 

2 

time-zone 

-3 

aylight  Saving 

0 

Sunrise 

7.218 

hours 

3 

Day 

21 

Hour 

8 

Solamoon 

13.272 

hours 

u 

4 

Month 

3 

Minute 

0 

Earliest  Start 

3 

Sunset 

191126 

hours 

5 

Year 

2010 

Second 

0 

Latest  Start 

18 

6 

Leap  Year 

0 

Riverine  forest  height 

7 

Boat  with= 

0.001242  miles 

2 

6.56167979 

feet 

0.0062  miles 

10 

8 

Boat  speeds 

3.106855  miles/hour 

5 

Km/h 

32.808  feet 

9 

Ratio  of  the  Earth= 

3,963.17  miles 

393.7  inch 

10 

11 

geographic  latitude 

geographic  longitude 

12 

Nodes 

D 

M 

S 

D 

M 

S 

ifl 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

13 

1 

Point  01 

-8 

-33.39 

-62 

-18 

-2.93 

1 

0.0062 

0.0062 

0.0062 

0.0062 

0.0062 

0.0062 

0.0062 

0.0062 

0.0062 

0.0062 

0.00 

: 

14 

2 

Point  02 

-8 

-56.77 

-62 

-20 

-10.73 

2 

0.0062 

0.0062 

0.0062 

0.0062 

0.0062 

0.0062 

0.0062 

0.0062 

0.0062 

0.0062 

0.00 

: 

15 

3 

Point  03 

-9 

-20.87 

-62 

-18 

-14.87 

3 

0.0062 

0.0062 

0.0062 

0.0062 

0.0062 

0.0062 

0.0062 

0.0062 

0.0062 

0.0062 

0.00 

: 

16 

4 

Point  04 

-10 

-10.27 

-62 

-16 

-27.98 

4 

0.0062 

0.0062 

0.0062 

0.0062 

0.0062 

0.0062 

0.0062 

0.0062 

0.0062 

0.0062 

0.00 

: 

17 

5 

Point  05 

-9 

-38.05 

-62 

-16 

-19.46 

5 

0.0062 

0.0062 

0.0062 

0.0062 

0.0062 

0.0062 

0.0062 

0.0062 

0.0062 

0.0062 

0.00 

: 

18 

6 

Point  06 

-8 

-30.17 

-62 

-15 

-30.58 

6 

0.0062 

0.0062 

0.0062 

0.0062 

0.0062 

0.0062 

0.0062 

0.0062 

0.0062 

0.0062 

0.00 

: 

19 

7 

Point  07 

-11 

-50.28 

-62 

-16 

-18.18 

7 

0.0062 

0.0062 

0.0062 

0.0062 

0.0062 

0.0062 

0.0062 

0.0062 

0.0062 

0.0062 

0.00 

: 

20 

8 

Point  08 

-9 

-50.01 

-62 

-14 

-17.58 

8 

0.0062 

0.0062 

0.0062 

0.0062 

0.0062 

0.0062 

0.0062 

0.0062 

0.0062 

0.0062 

0.00 

: 

21 

9 

Point  09 

-8 

-47.25 

-62 

-14 

-2.39 

9 

0.0062 

0.0062 

0.0062 

0.0062 

0.0062 

0.0062 

0.0062 

0.0062 

0.0062 

0.0062 

0.00 

: 

22 

10 

Point  10 

-10 

-57.86 

-62 

-14 

-29.54 

10 

0.0062 

0.0062 

0.0062 

0.0062 

0.0062 

0.0062 

0  0062 

0.0062 

0.0062 

0.0062 

0.00 

: 

23 

11 

Point  11 

-12 

-12.48 

-62 

-14 

-38.15 

11 

0.0062 

0.0062 

0.0062 

0.0062 

0.0062 

0.0062 

0.0062 

0.0062 

0.0062 

0.0062 

0.00 

: 

24 

12 

Point  12 

-10 

-17.27 

-62 

-12 

-42.80 

12 

0.0062 

0.0062 

0.0062 

0.0062 

0.0062 

0.0062 

0.0062 

0.0062 

0.0062 

0.0062 

0.00 

; 

25 

13 

Point  13 

-12 

-48.08 

-62 

-14 

-53.47 

13 

0.0062 

0.0062 

0.0062 

0.0062 

0.0062 

0.0062 

0.0062 

0.0062 

0.0062 

0.0062 

0.00 

: 

26 

14 

Point  14 

-13 

-35.56 

-62 

-15 

-2.00 

14 

0.0062 

0.0062 

0.0062 

0.0062 

0.0062 

0.0062 

0.0062 

0.0062 

0.0062 

0.0062 

0.00 

: 

27 

15 

Point  15 

-14 

-19.68 

-62 

-14 

-45.10 

15 

0.0062 

0.0062 

0.0062 

0.0062 

0.0062 

0.0062 

0.0062 

0  0062 

0.0062 

0.0062 

0.00 

: 

2S 

16 

Point  16 

-14 

-19.83 

-62 

-12 

-53.10 

16 

0.0062 

0.0062 

0.0062 

0.0062 

0.0062 

0.0062 

0.0062 

0.0062 

0.0062 

0.0062 

0.00 

; 

29 

17 

Point  17 

-13 

-57.84 

-62 

-12 

-19.15 

17 

0.0062 

0.0062 

0.0062 

0.0062 

0.0062 

0.0062 

0.0062 

0.0062 

0.0062 

0.0062 

0.00 

: 

30 

18 

Point  18 

-11 

-47.40 

-62 

-11 

-31.28 

18 

0.0062 

0.0062 

0.0062 

0.0062 

0.0062 

0.0062 

0.0062 

0.0062 

0.0062 

0.0062 

0.00 

: 

31 

19 

Point  19 

-15 

-58.30 

-62 

-12 

-58  26 

19 

0.0062 

0.0062 

0.0062 

0.0062 

0.0062 

0.0062 

0.0062 

0.0062 

0.0062 

0.0062 

0.00 

: 

32 

20 

Point  20 

-9 

-41.88 

-62 

-10 

-8.37 

20 

0.0062 

0.0062 

0.0062 

0.0062 

0.0062 

0.0062 

0.0062 

0.0062 

0.0062 

0.0062 

0.00 

: 

33 

21 

Point  21 

-18 

-9.03 

-62 

-11 

-52.29 

21 

0.0062 

0.0062 

0.0062 

0.0062 

0.0062 

0.0062 

0.0062 

0.0062 

0.0062 

0.0062 

0.00 

: 

34 

22 

Point  22 

-15 

-29.66 

-62 

-10 

-32.27 

22 

0.0062 

0.0062 

0.0062 

0.0062 

0.0062 

0.0062 

0.0062 

0.0062 

0.0062 

0.0062 

0.00 

: 

35 

23 

Point  23 

-14  -42.42 

-62 

-7 

-35.71 

23 

0.0062 

0.0062 

0.0062 

0.0062 

0.0062 

0.0062 

0.0062 

0.0062 

0.0062 

0.0062 

0.00 

: 

♦ 

36 

24 

Point  24 

-11 

-37.50 

-62 

-6 

-56.56 

0.0062 

0.0062 

0.0062 

0.0062 

input ,  Output  °«j _  _ I'M  mi 

Ready  S  (Ml,  100%  0  P  0 


Figure  122.  Input  Spreadsheet’s  Excel  screenshot, 
considering  a  normal  operational  scenario. 
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2 


Route  A 


Figure  123.  Theoretical  network  plot  showing  the  two  different  scenarios. 


Figure  124.  Output  Spreadsheet’s  Excel  screenshot, 
considering  a  normal  operational  scenario. 
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Figure  125.  Input  Spreadsheet’s  Excel  screenshot, 
considering  a  new  scenario  configuration. 
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Figure  126.  Output  Spreadsheet’s  Excel  screenshot, 
considering  a  new  scenario  configuration. 
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Synthesis 

The  investigation  performed  on  the  computer  program  against  the  theoretical 
model  proved  conclusively  its  feasibility  and  performance  in  assessing  move  risk  and 
survivability  under  stealthy  movement  in  jungle  environments. 

The  effects  of  path  surroundings,  such  as  riverine  forest  features,  as  well  as 
inherent  characteristics  of  boat,  essentially  its  speed  and  width,  influences  stealth 
behavior  movement.  Both  of  them  are  travel  start-time  dependent.  The  solar  noon  move 
experiments  greatly  increase  travel  risk,  since  the  sun  is  very  far  away  from  the  earth. 
This  is  an  essential  element  to  be  considered  in  the  transverse  planning  process,  as  the 
main  goal  is  to  avoid  enemy  detection. 

Another  important  consideration  is  referred  to  how  many  different  risk  assessment 
approaches  are  able  to  provide  concise  but  significant  insights  to  military  planners  about 
the  safety  routes  and  time-windows  framework.  At  this  time  the  military  objective  is  the 
primary  concern  when  choosing  the  risk  evaluation  method.  Further,  using  properly  the 
chosen  tool  is  important.  For  example,  selecting  distinct  time-range  to  start-time  provides 
the  user  flexibility  in  identifying  more  convenient  period  to  cross  the  network. 

To  sum  up,  the  model  reasonably  achieved  the  purposed  objective,  providing  a 
method  to  accurately  calculate  both  routes  and  period  of  times  that  minimize  detection 
probability  taken  into  account  the  assumptions  previously  made. 
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V.  Discussion 


Introduction 

This  research  provided  a  computational  tool  to  calculate  optimal  routes  and  time 
window  framework  to  supply  a  deployed  jungle  combat  unit  under  restriction  of  an 
omnidirectional  sensing  searcher.  This  final  chapter  serves  to  sum  up  the  entirety  of  this 
study  as  well  as  answer  the  research  questions. 

Insights 

This  research  revealed  valuable  insights  for  future  studies  concerning  combat 
logistics  research  topics,  especially  those  concerned  with  travel  risk  assessment  under 
enemy  detection  threat.  The  main  insights  are  found  in  the  following  paragraphs. 

First,  it  was  evident  the  benefits  of  the  concealment  schemes  employment  for 
avoiding  detection  and  movement  tracking,  in  terms  of  position  and  speed  taking.  This 
technique,  whose  efficiency  was  historically  proven,  is  an  inherently  military  activity  in 
which  the  major  goal  is  to  annul  or  reduce  materiel  and  personnel  loss  due  to  enemy 
target  acquisition.  However,  it  also  has  been  efficiently  employed  on  civilian  sectors, 
which  have  pursued  performance  improvements  in  artificial  intelligence  devices. 
Particularly  in  jungle  warfare,  because  of  the  peculiar  environmental  conditions,  hiding 
movements  are  relatively  simple  since  movement  is  able  to  take  advantage  of  shadows 
provided  by  botanical  features.  However,  since  this  method  conjugates  both  internal  and 
external  movement  attributes  such  as  boat  dimension  and  forest  features,  respectively, 
optimal  solutions  arise  from  a  balance  between  these  two  perspectives,  which  is  quite 
difficult  to  achieve. 
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Second,  changes  in  risk  quantifier  parameters  don’t  affect  proportionally  and 
equally  the  boat  vulnerability  through  the  route.  Increasing  in  boat  speed,  keeping  its 
feasibility,  for  instance,  demonstrates  low  impact  in  reducing  boat  vulnerability,  since  the 
distances  between  nodes  are  big  and  use  of  a  powerful  engine  is  impracticable  because  of 
disturbing  stealthy  behavior.  On  the  other  hand,  a  small  variation  on  boat  width  produces 
large  consequences  for  risk  travel,  being  an  important  factor  to  be  considered  during  the 
stealthy  supply  watercraft  design  process.  Also,  though  riverine  vegetation  has  a  weighty 
influence  on  offering  distinct  shadow  areas,  there  is  no  control  over  their  facet, 
suggesting  it  has  to  be  seen  as  a  guideline  to  develop  route  strategies,  travel  machines 
models  or  depot  location  arrangements. 

Third,  the  possibility  of  providing  different  analytic  approaches  to  the  same 
operational  scenario  contributes  to  a  military  planner,  ability  to  attain  simultaneously 
tactical  and  logistics  objectives,  since  the  outcomes  are  compatibly  displayed  with  their 
goals.  If  a  camp  commander  emphasizes  the  value  of  a  zero-detection  enemy  travel,  in 
transporting  ammunitions,  drugs  or  wounded  soldiers,  for  example,  the  relative  risk  and 
survivability  analysis  are  the  best  tools  to  be  employed.  However,  if  he  has  more  freedom 
of  action,  both  in  space  and  time  dimensions,  to  conduct  his  operations,  as  well  as  he 
pursues  a  better  tactic-logistic  integration,  the  weighted-survivability  analysis  is  the 
outcome  that  can  better  help  him  during  the  trade-off  process  between  available  choices: 
faster  delivery  but  more  exposure  of  supplies,  boats  and  its  operator,  or  sacrificing  the 
troops  readiness  against  a  guaranteed  resupply. 
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Assumptions  and  Limitations 

Some  fundamental  assumptions  and  limitations  were  presented  in  this  study  since 
its  objective  was  to  give  quantitative  visibility  to  solve  a  realistic,  specific  and  well- 
defined  military  problem:  supply  combat  units  in  jungle  environment  conditions. 

The  first  logical  constraint  is  related  to  set  up  a  combat  situation  framework  and 
its  level  of  detail.  These  most  typically  occur  when  there  is  not  sufficient  information 
about  the  battle  dynamics  to  simulate,  or  when  a  model  is  to  be  used  to  calculate  a 
measure  of  effectiveness  of  analysis  which  may  not  be  directly  related  to  specific  combat 
processes.  As  a  result,  in  this  research,  the  major  path  watercourse  characteristics  as  well 
as  height  of  riverine  vegetation  was  considered  constant,  both  in  time  and  space 
dimensions,  and  whose  values  are  simply  estimates  of  the  existent  real-world  entities. 
Also,  since  there  is  no  information  about  the  enemy  behavior,  its  tactics,  operational 
objectives,  and  available  means  to  achieve  them  as  well  as  their  real  capabilities,  the 
worst  scenario  was  modeled,  in  which  enemy  forces  have  total  visibility  of  battlefield. 

Although  the  model  accurately  calculates  both  safety  routes  and  better  start-time 
to  perform  a  stealthy  navigation  to  the  deployed  fight  unit  locations,  these  outcomes 
referenced  the  existence  of  solar  radiation.  Since  the  algorithm  used  in  this  research  was  a 
sun  position  calculator,  the  use  of  different  environment  light  sources,  such  as  luminosity 
of  the  moon,  won’t  provide  the  same  performance,  for  logical  reasons.  Therefore,  this 
methodology  is  valid  for  travel  performed  during  daylight. 

As  already  mentioned,  the  dynamics  of  the  watercourse  network  wasn’t  modeled. 
In  the  real  world,  the  amount  of  variables  that  affect  the  aerial  exposure  level  on  the 
jungle  environment  is  beyond  the  tree  height  and  move  features.  The  environmental 
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effect  of  fog  at  sunrise  and  sunset  due  to  high  humidity  tropical  rain  forests,  as  well  as  the 
excessive  condensation  on  dry  seasons,  was  omitted,  though  both  of  them  are  low- 
visibility  makers.  Additionally,  the  river  attributes  (flow  speed  and  depth)  were  also 
completely  neglected  in  addition  to  the  main  riverine  forest  aspect  concerning  its  fallen 
posture  over  streams,  which  reduce  boat  vulnerability  due  to  large  shadow  areas 
originated  near  to  its  borders. 

Another  relevant  aspect  not  modeled  in  this  research  is  the  return  trip  of  the  boat 
to  the  depot.  The  model  ignores  this  action  as  well  as  its  implications  for  future  missions. 
Also,  since  the  searcher  has  an  omnidirectional  sensing  of  observation,  move  evasive 
actions  or  deliberate  route  changes  against  the  previous  planned  ones,  after  aerial  enemy 
detection,  were  completely  omitted. 

Future  Research 

While  this  research  serves  to  address  some  important  questions  related  to  stealthy 
supply  navigation  under  jungle  warfare  constraints,  it  by  no  means  answers  all  of  them. 
Thus,  suggestions  for  future  analysis  are  the  following: 

a)  Keeping  both  a  single-source  and  destination  (combat  unit),  assess  the  performance  of 
the  Bellman-Ford  algorithm  in  dealing  with  stealthy  navigation  problem; 

b)  Increase  the  number  of  units  deployed  on  battlefield  and,  as  a  result,  evaluate  the 
optimal  start-time,  boat  dimensions; 

c)  Also  increase  the  number  of  different  moves  capable  to  support  several  locations, 
leaving  the  depot  either  simultaneously  or  at  distinct  times; 
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d)  Perform  assessment  under  conditions  of  changing  in  time  dependent  environment 
features  such  as  flow  speed,  on  both  direction  and  its  intensity,  and  depth  river, 
including  seasonal  alterations,  which  can  affect  boat  transverse  capability; 

e)  Given  a  network,  as  well  as  the  operations  theater  boundaries,  evaluate  the  optimal 
fixed  depot  location  able  to  minimize  the  total  travel  risk  between  an  origin  and  the 
deployed  fight  forces; 

f)  Include  the  moon  luminosity  influence  on  exposure  level  calculations  during  network 
night  crossing; 

g)  Conduct  a  random  enemy  aerial  search  strategy  on  the  vulnerability  evaluation,  since 
the  present  study  focused  on  omnidirectional  sensing  observer,  which  is  a  pessimistic 
as  well  as  unrealistic  (theoretical)  construct; 

h)  Develop  other  combat  logistics  effectiveness  measures  for  assessing  travel  network 
vulnerability  that  include  fuel  consumption  rate,  defense  posture,  water  searcher  in 
addition  to  aerial  ones. 
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Appendix  A:  Geographic  Coordinates  of  the  Theoretical  Network  Locations 


Geographic  Latitude 

Geographic  Lonj 

*itude 

Nodes 

Degree 

Minute 

Second 

Degree 

Minute 

Second 

1 

Node  01 

-1 

-8 

-33.39 

-62 

-18 

-2.93 

2 

Node  02 

-1 

-8 

-56.77 

-62 

-20 

-10.73 

3 

Node  03 

-1 

-9 

-20.87 

-62 

-18 

-14.87 

4 

Node  04 

-1 

-10 

-10.27 

-62 

-16 

-27.98 

5 

Node  05 

-1 

-9 

-38.05 

-62 

-16 

-19.46 

6 

Node  06 

-1 

-8 

-30.17 

-62 

-15 

-30.58 

7 

Node  07 

-1 

-11 

-50.28 

-62 

-16 

-18.18 

8 

Node  08 

-1 

-9 

-50.01 

-62 

-14 

-17.58 

9 

Node  09 

-1 

-8 

-47.25 

-62 

-14 

-2.39 

10 

Node  10 

-1 

-10 

-57.86 

-62 

-14 

-29.54 

11 

Node  11 

-1 

-12 

-12.48 

-62 

-14 

-38.15 

12 

Node  12 

-1 

-10 

-17.27 

-62 

-12 

-42.80 

13 

Node  13 

-1 

-12 

-48.08 

-62 

-14 

-53.47 

14 

Node  14 

-1 

-13 

-35.56 

-62 

-15 

-2.00 

15 

Node  15 

-1 

-14 

-19.68 

-62 

-14 

-45.10 

16 

Node  16 

-1 

-14 

-19.83 

-62 

-12 

-53.10 

17 

Node  17 

-1 

-13 

-57.84 

-62 

-12 

-19.15 

18 

Node  18 

-1 

-11 

-47.40 

-62 

-11 

-31.28 

19 

Node  19 

-1 

-15 

-58.30 

-62 

-12 

-58.26 

20 

Node  20 

-1 

-9 

-41.88 

-62 

-10 

-8.37 

21 

Node  21 

-1 

-18 

-9.03 

-62 

-11 

-52.29 

22 

Node  22 

-1 

-15 

-29.66 

-62 

-10 

-32.27 

23 

Node  23 

-1 

-14 

-42.42 

-62 

-7 

-35.71 

24 

Node  24 

-1 

-11 

-37.50 

-62 

-6 

-56.56 

25 

Node  25 

-1 

-10 

-56.82 

-62 

-6 

-25.95 

26 

Node  26 

-1 

-18 

-15.82 

-62 

-10 

-14.37 

27 

Node  27 

-1 

-17 

-14.76 

-62 

-4 

-55.59 

28 

Node  28 

-1 

-16 

-46.10 

-62 

-8 

-59.06 

29 

Node  29 

-1 

-16 

-54.72 

-62 

-6 

-33.52 

30 

Node  30 

-1 

-19 

-29.02 

-62 

-7 

-22.95 

31 

Node  31 

-1 

-18 

-31.40 

-62 

-6 

-45.54 

32 

Node  32 

-1 

-15 

-24.97 

-62 

-5 

-54.04 

33 

Node  33 

-1 

-14 

-41.00 

-62 

-4 

-51.02 

34 

Node  34 

-1 

-12 

-0.55 

-62 

-5 

-18.47 

35 

Node  35 

-1 

-12 

-11.75 

-62 

-4 

-32.33 

36 

Node  36 

-1 

-14 

-19.11 

-62 

-2 

-42.22 

37 

Node  37 

-1 

-15 

-25.17 

-62 

-3 

-43.39 
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38 

Node  38 

-1 

-16 

28.01 

-62 

-2 

-59.34 

39 

Node  39 

-1 

-18 

-43.43 

-62 

-4 

-56.95 

40 

Node  40 

-1 

-19 

-55.78 

-62 

-4 

-58.65 

41 

Node  41 

-1 

-20 

-35.41 

-62 

-5 

-32.12 

42 

Node  42 

-1 

-20 

-47.91 

-62 

-5 

-56.56 

43 

Node  43 

-1 

-17 

-7.34 

-62 

-1 

-15.74 

44 

Node  44 

-1 

-18 

-14.14 

-62 

-1 

-47.97 

45 

Node  45 

-1 

-22 

-21.48 

-62 

-2 

-43.31 

46 

Node  46 

-1 

-21 

-31.52 

-62 

-1 

-0.17 

47 

Node  47 

-1 

-20 

-59.15 

-62 

0 

-42.50 

48 

Node  48 

-1 

-15 

-56.90 

-61 

-58 

-28.26 

49 

Node  49 

-1 

-17 

-53.75 

-61 

-59 

-36.82 

50 

Node  50 

-1 

-18 

-31.37 

-62 

0 

-7.76 

51 

Node  51 

-1 

-18 

-10.76 

-61 

-55 

-56.94 

52 

Node  52 

-1 

-20 

-59.00 

-61 

-57 

-13.60 

53 

Node  53 

-1 

-22 

-50.83 

-61 

-58 

-56.76 

54 

Node  54 

-1 

-23 

-34.97 

-61 

-59 

-35.09 

55 

Node  55 

-1 

-22 

-15.64 

-61 

-56 

-23.47 

56 

Node  56 

-1 

-21 

-19.71 

-61 

-55 

-39.29 

57 

Node  57 

-1 

-20 

-11.71 

-61 

-54 

-39.91 

58 

Node  58 

-1 

-23 

-55.82 

-61 

-55 

-33.35 

59 

Node  59 

-1 

-22 

-48.18 

-61 

-53 

-44.36 

60 

Node  60 

-1 

-23 

-44.12 

-61 

-53 

-56.14 

61 

Node  61 

-1 

-24 

-31.23 

-61 

-54 

-4.97 

62 

Node  62 

-1 

-23 

-23.71 

-61 

-50 

-24.02 

63 

Node  63 

-1 

-24 

-7.87 

-61 

-50 

-21.08 

64 

Node  64 

-1 

-25 

-0.86 

-61 

-50 

-29.91 
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Appendix  B:  Model  Code 


Option  Explicit 

Function  radToDeg(angleRad) 

'Convert  Angles  in  Radians  to  ones  in  Degrees' 

radToDeg  =  (180#  *  angleRad/  Application.WorksheetFunction.pi()) 

End  Function 

Function  degToRad(angleDeg) 

'Convert  Angles  in  Degrees  to  ones  in  Radians' 

degToRad  =  (Application.WorksheetFunction.pi()  *  angleDeg/  180#) 

End  Function 


Function  calcJD(year,  month,  day) 

'Calculate  Julian  Day  from  Calendary  Day' 

Dim  A  As  Double  'Declarate  Auxiliar  Variable  "A"  to  Calculate  Julian  Date' 

Dim  B  As  Double  'Declarate  Auxiliar  Variable  "B"  to  Calculate  Julian  Date' 

Dim  JD  As  Double  'Declarate  the  Julian  Day' 

'Month:  January  =  1;  Day:  from  1  to  31' 

If  (month  <=  2)  Then 
year  =  year  -  1 
month  =  month  +12 
End  If 

A  =  Application.  WorksheetFunction.Floor(year/  100,  1) 

B  =  2  -  A  +  Application. WorksheetFunction.Floor(A  /  4,  1) 

JD  =  Application.WorksheetFunction.Floor(365.25  *  (year  +  4716),  1)  + 
Application.WorksheetFunction.Floor(30.6001  *  (month  +  1),  1)  +  day  +  B  -  1524.5 
calcJD  =  JD 

'Correct  Years  and  Months  values  to  the  Current  Calendary  pattern' 

If  month  =13  Then 
month  =  1 
year  =  year  +  1 
End  If 

If  month  =14  Then 
month  =  2 
year  =  year  +  1 
End  If 

End  Function 
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Function  calcTimeJulianCent(JD) 

'Calculate  the  Number  of  Julian  centuries  since  J2000.0' 

Dim  t  As  Double  'Declarate  the  number  of  Julian  centuries  since  J2000.0' 

t  =  (JD  -  2451545#)  /  36525# 
calcTimeJulianCent  =  t 

End  Function 

Function  calcJDFromJulianCent(t) 

'Convert  J2000.0  to  Julian  Date' 

Dim  JD  As  Double  'Declarate  the  Julian  Day' 

JD  =  t  *  36525#  +  2451545# 
calcJDFromJulianCent  =  JD 

End  Function 


Function  calcGeomMeanLongSun(t) 

'Calculate  the  Geometric  Mean  Longitude  of  the  Sun  in  Degrees' 

Dim  10  As  Double  'Declarate  the  Geometric  Mean  Longitude  of  the  Sun' 

10  =  280.46646  + 1  *  (36000.76983  +  0.0003032  *  t) 

Do 

If  (10  <=  360)  And  (10  >=  0)  Then  Exit  Do 
If  10  >  360  Then  10  =  10-360 
If  10  <0  Then  10  =  10  +  360 
Loop 

calcGeomMeanLongSun  =  10 
End  Function 


Function  calcGeomMeanAnomalySun(t) 

'Calculate  the  Geometric  Mean  Anomaly  of  the  Sun  in  Degrees' 

Dim  m  As  Double  'Declarate  Geometric  Mean  Anomaly  of  the  Sun' 

m  =  357.5291 1  + 1  *  (35999.05029  -  0.0001537  *  t) 
calcGeomMeanAnomalySun  =  m 

End  Function 


Function  calcEccentricityEarthOrbit(t) 

'Calculate  the  Eccentricity  of  the  Earth's  orbit' 

Dim  e  As  Double  'Declarate  the  Eccentricity  of  the  Earth's  Orbit' 
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e  =  0.016708634  - 1  *  (0.000042037  +  0.0000001267  *  t) 
calcEccentricityEarthOrbit  =  e 

End  Function 


Function  calcSunEqOfCenter(t) 

'Calculate  the  Equation  of  Center  of  the  Sun  in  Degrees' 


Dim  m  As  Double 
Dim  mrad  As  Double 
Dim  sinm  As  Double 
Dim  sin2m  As  Double 
Dim  sin3m  As  Double 
Dim  c  As  Double 


'Declarate  Geometric  Mean  Anomaly  of  the  Sun  in  Degrees' 
'Declarate  Geometric  Mean  Anomaly  of  the  Sun  in  Radians' 
'Declarate  the  Sine  of  the  "m"  Angle' 

’Declarate  the  Sin  of  the  "2m"  Angle' 

’Declarate  the  Sin  of  the  "3m"  Angle' 

'Declarate  the  Equation  of  Center  of  the  Sun  in  Degrees' 


m  =  calcGeomMeanAnomalySun(t) 

mrad  =  degToRad(m) 

sinm  =  Sin(mrad) 

sin2m  =  Sin(mrad  +  mrad) 

sin3m  =  Sin(mrad  +  mrad  +  mrad) 


c  =  sinm  *  (1.914602  - 1  *  (0.004817  +  0.000014  *  t))  +  sin2m  *  (0.019993  -  0.000101  *  t)  + 
sin3m  *  0.000289 

calcSunEqOfCenter  =  c 

End  Function 


Function  calcSunTrueLong(t) 

'Calculate  the  True  Longitude  of  the  Sun  in  Degrees' 

Dim  10  As  Double  'Declarate  the  Geometric  Mean  Longitude  of  the  Sun' 
Dim  c  As  Double  'Declarate  the  Equation  of  Center  of  the  Sun  in  Degrees' 

Dim  O  As  Double  ’Declarate  the  True  Longitude  of  the  Sun  in  Degrees' 

10  =  calcGeomMeanLongSun(t) 
c  =  calcSunEqOfCenter(t) 

O  =  10  +  c 

calcSunTrueLong  =  O 
End  Function 


Function  calc  SunTrue Anomaly (t) 

'Calculate  the  True  Anonamaly  of  the  Sun  in  Degrees' 


Dim  m  As  Double 
Dim  c  As  Double 
Dim  v  As  Double 


'Declarate  Geometric  Mean  Anomaly  of  the  Sun  in  Degrees' 
'Declarate  the  Equation  of  Center  of  the  Sun  in  Degrees' 
'Declarate  the  True  Anonamaly  of  the  Sun  in  Degrees' 
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m  =  calcGeomMeanAnomalySun(t) 
c  =  calcSunEqOfCenter(t) 


v  =  m  +  c 

calc  SunTrue  Anomaly  =  v 
End  Function 


Function  calcSunRadVector(t) 

'Calculate  the  Distance  to  the  Sun  in  Astronomical  Units  (UA)' 

Dim  v  As  Double  'Declarate  the  True  Anonamaly  of  the  Sun  in  Degrees' 

Dim  e  As  Double  'Declarate  the  Eccentricity  of  the  Earth's  orbit' 

Dim  R  As  Double  'Declarate  the  Distance  to  the  Sun  in  Astronomical  Units  (UA)' 

v  =  calcSunTrueAnomaly(t) 
e  =  calcEccentricityEarthOrbit(t) 

R  =  (1.000001018  *  (1  -  e  *  e))  /  (1  +  e  *  Cos(degToRad(v))) 
calcSunRadVector  =  R 

End  Function 


Function  calcSunApparentLong(t) 

'Calculate  the  Apparent  Longitude  of  the  Sun  in  Degrees' 

Dim  O  As  Double  'Declarate  the  True  Longitude  of  the  Sun  in  Degrees' 

Dim  omega  As  Double  'Declarate  the  Geocentric  Longitude  of  the  Sun  in  Degree' 
Dim  lambda  As  Double  'Declarate  the  Apparent  Longitude  of  the  Sun  in  Degrees' 

O  =  calcSunTrueLong(t) 

omega  =  125.04  -  1934.136  *  t 

lambda  =  O  -  0.00569  -  0.00478  *  Sin(degToRad(omega)) 
calc  Sun  ApparentLong  =  lambda 

End  Function 


Function  calcMeanObliquityOfEcliptic(t) 

'Calculate  the  Mean  Obliquity  of  the  Ecliptic' 

Dim  seconds  As  Double  'Declarate  seconds' 

Dim  eO  As  Double  'Declarate  the  Mean  Obliquity  of  the  Ecliptic' 

seconds  =  21.448  - 1  *  (46.815  + 1  *  (0.00059  - 1  *  (0.001813))) 
eO  =  23#  +  (26#  +  (seconds  /  60#))  /  60# 
calcMeanObliquityOfEcliptic  =  eO 

End  Function 
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Function  calcObliquityCorrection(t) 

'Calculate  the  Corrected  Obliquity  of  the  Ecliptic  in  Degrees' 

Dim  eO  As  Double  'Declarate  the  Mean  Obliquity  of  the  Ecliptic' 

Dim  omega  As  Double  'Declarate  the  Geocentric  Longitude  of  the  Sun  in  Degree’ 

Dim  e  As  Double  'Declarate  the  Eccentricity  of  the  Earth's  orbit' 

eO  =  calcMeanObliquityOfEcliptic(t) 

omega  =  125.04  -  1934.136  *  t 
e  =  eO  +  0.00256  *  Cos(degToRad(omega)) 
calcObliquityCorrection  =  e 

End  Function 


Function  calcSunRtAscension(t) 

'Calculate  the  Right  Ascension  of  the  Sun  in  Degrees' 


Dim  e  As  Double 
Dim  lambda  As  Double 
Dim  tananum  As  Double 
Dim  tanadenom  As  Double 
Dim  alpha  As  Double 


'Declarate  the  Eccentricity  of  the  Earth's  orbit' 

'Declarate  the  Apparent  Longitude  of  the  Sun  in  Degrees' 
'Declarate  the  Numerator  of  the  Auxiliar  Variable  "a"' 
'Declarate  the  Denominator  of  the  Auxiliar  Variable  "a"' 
'Declarate  the  Right  Ascension  of  the  Sun  in  Degrees' 


e  =  calcObliquityCorrection(t) 
lambda  =  calcSunApparentLong(t) 


tananum  =  (Cos(degToRad(e))  *  Sin(degToRad(lambda))) 
tanadenom  =  (Cos(degToRad(lambda))) 


’Original  NOAA  code  using  javascript  Math.Atan2(y,x)  convention:' 

’a)  var  alpha  =  radToDeg(Math.atan2(tananum,  tanadenom));' 

’b)  alpha  =  radToDeg(Application.WorksheetFunction.Atan2(tananum,  tanadenom)).’ 


’Translated  using  Excel  VBA  Application.WorksheetFunction.Atan2(x,y)  convention: 
alpha  =  radToDeg(Application.WorksheetFunction.Atan2(tanadenom,  tananum)) 


calcSunRtAscension  =  alpha 


End  Function 


Function  calcSunDeclination(t) 

’Calculate  the  Declination  of  the  Sun  in  Degrees' 


Dim  e  As  Double 
Dim  lambda  As  Double 
Dim  sint  As  Double 

Dim  theta  As  Double 


'Declarate  the  Eccentricity  of  the  Earth's  orbit' 

'Declarate  the  Apparent  Longitude  of  the  Sun  in  Degrees' 

'Declarate  the  Auxiliar  Variable  which  represents  the  produt  of  Earth's 
orbit  Eccentricity  and  the  Sun  Longitude' 

’Declarate  the  Declination  of  the  Sun  in  Degrees' 
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e  =  calcObliquityCorrection(t) 
lambda  =  calcSunApparentLong(t) 

sint  =  Sin(degToRad(e))  *  Sin(degToRad(lambda)) 

theta  =  radToDeg(Application.WorksheetFunction.Asin(sint)) 

calcSunDeclination  =  theta 

End  Function 


Function  calcEquationOfTime(t) 

'Calculate  the  Equation  of  Time  in  Degrees' 


Dim  epsilon  As  Double 
Dim  10  As  Double 
Dim  e  As  Double 
Dim  m  As  Double 
Dim  y  As  Double 

Dim  sin210  As  Double 

Dim  sinm  As  Double 
Dim  cos210  As  Double 

Dim  sin410  As  Double 

Dim  sin2m  As  Double 

Dim  Etime  As  Double 


'Declarate  the  Corrected  Obliquity  of  the  Ecliptic  in  Degrees 
'Declarate  the  Geometric  Mean  Longitude  of  the  Sun' 

'Declarate  the  Eccentricity  of  the  Earth's  Orbit' 

'Declarate  the  Geometric  Mean  Anomaly  of  the  Sun  in  Degrees' 
'Declarate  the  Auxiliar  Variable  whic  represents  the  Tangent  of 
Epsilon  Value' 

'Declarate  the  Sine  Squared  of  the  Geometric  Mean  Longitude  of  the 
Sun' 

'Declarate  the  Sine  of  the  Geometric  Mean  Anomaly  of  the  Sun' 
'Declarate  the  Cosine  Squared  of  the  Geometric  Mean  Longitude  of 
the  Sun' 

'Declarate  the  Cosine  Fouth  of  the  Geometric  Mean  Longitude  of  the 
Sun' 

'Declarate  the  Sin  Squared  of  the  Geometric  Mean  Anomaly  of  the 
Sun' 

'Declarate  the  Equation  of  Time  in  Degrees' 


epsilon  =  calcObliquityCorrection(t) 
10  =  calcGeomMeanLongSun(t) 
e  =  calcEccentricityEarthOrbit(t) 
m  =  calcGeomMeanAnomalySun(t) 


y  =  Tan(degToRad(epsilon)  /  2#) 
y  =  y  A  2 


sin210  =  Sin(2#  *  degToRad(lO)) 
sinm  =  Sin(degToRad(m)) 
cos210  =  Cos(2#  *  degToRad(lO)) 
sin410  =  Sin(4#  *  degToRad(lO)) 
sin2m  =  Sin(2#  *  degToRad(m)) 


Etime  =  y  *  sin210  -  2#  *  e  *  sinm  +  4#  *  e  *  y  *  sinm  *  cos210  -  0.5  *  y  *  y  *  sin410  -  1 .25  *  e 
*  e  *  sin2m 


calcEquationOfTime  =  radToDeg(Etime)  *  4# 
End  Function 
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Function  calcHourAngleDawn(lat,  solarDec,  solardepression) 
'Calculate  the  Local  Hour  Angle  of  the  Sun  at  Dawn  in  Radians' 


Dim  latRad  As  Double 
Dim  sdRad  As  Double 
Dim  HAarg  As  Double 

Dim  HA  As  Double 


'Declarate  the  Local  Latitude  in  Radians' 

'Declarate  the  Solar  Declination  Angle  in  Degrees' 

'Declarate  the  Auxiliar  Variable  to  Calculate  Local  Hour  Angle  of  t 
the  Sun  at  Dawn' 

'Declarate  the  Local  Hour  Angle  of  the  Sun  at  Dawn' 


latRad  =  degToRad(lat) 
sdRad  =  degToRad(solarDec) 


HAarg  =  (Cos(degToRad(90  +  solardepression))  /  (Cos(latRad)  *  Cos(sdRad))  -  Tan(latRad)  * 
Tan(sdRad)) 


HA  =  (Application.WorksheetFunction.Acos(Cos(degToRad(90  +  solardepression))  / 
(Cos(latRad)  *  Cos(sdRad))  -  Tan(latRad)  *  Tan(sdRad))) 


calcHourAngleDawn  =  HA 


End  Function 


Function  calcHourAngleSunrise(lat,  solarDec) 

'Calculate  the  Local  Hour  Angle  of  the  Sun  at  Sunrise  in  Radians' 


Dim  latRad  As  Double 
Dim  sdRad  As  Double 
Dim  HAarg  As  Double 

Dim  HA  As  Double 


'Declarate  the  Local  Latitude  in  Radians 
'Declarate  the  Solar  Declination  Angle  in  Radians' 

'Declarate  the  Auxiliar  Variable  to  Calculate  Local  Hour  Angle  of  the 
Sun  at  Sunrise' 

'Declarate  the  Local  Hour  Angle  of  the  Sun  at  Sunrise  in  Radians' 


latRad  =  degToRad(lat) 
sdRad  =  degToRad(solarDec) 


HAarg  =  (Cos(degToRad(90.833))  /  (Cos(latRad)  *  Cos(sdRad))  -  Tan(latRad)  *  Tan(sdRad)) 


HA  =  (Application.WorksheetFunction.Acos(Cos(degToRad(90.833))  /  (Cos(latRad) 
Cos(sdRad))  -  Tan(latRad)  *  Tan(sdRad))) 


calcHourAngleSunrise  =  HA 


End  Function 


Function  calcHour  Angle  Sunset(lat,  solarDec) 

'Calculate  the  Local  Hour  Angle  of  the  Sun  at  Sunset  in  Radians' 


Dim  latRad  As  Double 
Dim  sdRad  As  Double 
Dim  HAarg  As  Double 


'Declarate  the  Local  Latitude  in  Radians 
'Declarate  the  Solar  Declination  Angle  in  Radians' 

'Declarate  the  Auxiliar  Variable  to  Calculate  Local  Hour  Angle  of 
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the  Sun  at  Sunset' 

Dim  HA  As  Double  'Declarate  the  Local  Hour  Angle  of  the  Sun  at  Sunset  in  Radians' 

latRad  =  degToRad(lat) 
sdRad  =  degToRad(solarDec) 

HAarg  =  (Cos(degToRad(90.833))  /  (Cos(latRad)  *  Cos(sdRad))  -  Tan(latRad)  *  Tan(sdRad)) 

HA  =  (Application.WorksheetFunction.Acos(Cos(degToRad(90.833))  /  (Cos(latRad)  * 
Cos(sdRad))  -  Tan(latRad)  *  Tan(sdRad))) 

calcHourAngleSunset  =  -HA 

End  Function 


Function  calcHourAngleDusk(lat,  solarDec,  solardepression) 
'Calculate  the  Local  Hour  Angle  of  the  Sun  at  Dusk' 


Dim  latRad  As  Double 
Dim  sdRad  As  Double 
Dim  HAarg  As  Double 

Dim  HA  As  Double 


'Declarate  the  Local  Latitude  in  Radians 
'Declarate  the  Solar  Declination  Angle  in  Radians' 

'Declarate  the  Auxiliar  Variable  to  Calculate  Local  Hour  Angle  of  the 
Sun  at  Dusk’ 

'Declarate  the  Local  Hour  Angle  of  the  Sun  at  Dusk  in  Radians 


latRad  =  degToRad(lat) 
sdRad  =  degToRad(solarDec) 


HAarg  =  (Cos(degToRad(90  +  solardepression))  /  (Cos(latRad)  *  Cos(sdRad))  -  Tan(latRad)  * 
Tan(sdRad)) 


HA  =  (Application.WorksheetFunction.Acos(Cos(degToRad(90  +  solardepression))  / 
(Cos(latRad)  *  Cos(sdRad))  -  Tan(latRad)  *  Tan(sdRad))) 

calcHourAngleDusk  =  -HA 

End  Function 


Function  calcDawnUTC(JD,  latitude,  longitude,  solardepression) 

'Calculate  the  Universal  Coordinated  Time  of  Dawn  for  a  given  Day  and  Location' 


Dim  t  As  Double 
Dim  eqtime  As  Double 
Dim  solarDec  As  Double 
Dim  hourangle  As  Double 

Dim  delta  As  Double 

Dim  timeDiff  As  Double 
Dim  timeUTC  As  Double 


'Declarate  the  number  of  Julian  centuries  since  J2000.0' 
'Declarate  the  Equation  of  Time  in  Degrees' 

'Declarate  the  Solar  Declination  in  Degrees' 

'Declarate  the  Local  Hour  Angle  at  Sunrise,  given  a  Latitude  and 
Solar  Declination' 

'Declarate  the  difference  between  Local  Longitude  and  Local 
Hour  Angle  at  Sunrise,  in  Radians' 

'Declarate  the  4  times  Delta  Value' 

'Declarate  the  Universal  Coordinated  Time  of  Dawn  in  Minutes' 
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Dim  newt  As  Double 


'Declarate  the  new  Mt' 


t  =  calcTimeJulianCent(JD) 

’First  pass  to  approximate  Sunrise’ 
eqtime  =  calcEquationOfTime(t) 
solarDec  =  calcSunDeclination(t) 
hourangle  =  calcHourAngleSunrise(latitude,  solarDec) 

delta  =  longitude  -  radToDeg(hourangle) 

timeDiff  =  4  *  delta 

timeUTC  =  720  +  timeDiff  -  eqtime 

’Second  pass  includes  fractional  Julian  Day  in  Gamma  calculation’ 
newt  =  calcTimeJulianCent(calcJDFromJulianCent(t)  +  timeUTC  /  1440#) 
eqtime  =  calcEquationOfTime(newt) 
solarDec  =  calcSunDeclination(newt) 

hourangle  =  calcHourAngleDawn(latitude,  solarDec,  solardepression) 

delta  =  longitude  -  radToDeg(hourangle) 

timeDiff  =  4  *  delta 

timeUTC  =  720  +  timeDiff  -  eqtime 

calcDawnUTC  =  timeUTC 

End  Function 


Function  calcSunriseUTC(JD,  latitude,  longitude) 

’Calculate  the  Universal  Coordinated  Time  (UTC)  of  Sunrise  for  the  given  Day  and  at  the  given 
Location  on  Earth' 

Dim  t  As  Double 
Dim  eqtime  As  Double 
Dim  solarDec  As  Double 
Dim  hourangle  As  Double 

Dim  delta  As  Double 

Dim  timeDiff  As  Double 
Dim  timeUTC  As  Double 

Dim  newt  As  Double 

t  =  calcTimeJulianCent(JD) 

'First  pass  to  approximate  Sunrise’ 
eqtime  =  calcEquationOfTime(t) 
solarDec  =  calcSunDeclination(t) 
hourangle  =  calcHourAngleSunrise(latitude,  solarDec) 


'Declarate  the  number  of  Julian  centuries  since  J2000.0’ 
'Declarate  the  Equation  of  Time  in  Minutes' 

'Declarate  the  Solar  Declination  in  Degrees' 

'Declarate  the  Local  Hour  Angle  at  Sunrise,  given  a  Latitude  and 
Solar  Declination’ 

'Declarate  the  difference  between  Local  Longitude  and  Local 
Hour  Angle  at  Sunrise,  in  Radians' 

'Declarate  4  times  Delta  Value' 

'Declarate  the  Universal  Coordinated  Time  of  Sunrise  in 
Minutes' 

'Declarate  the  new  "t"' 
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delta  =  longitude  -  radToDeg(hourangle) 

timeDiff  =  4  *  delta 

timeUTC  =  720  +  timeDiff  -  eqtime 

'Second  pass  includes  fractional  Julian  Day  in  Gamma  calculation' 

newt  =  calcTimeJulianCent(calcJDFromJulianCent(t)  +  timeUTC  /  1440#) 

eqtime  =  calcEquationOfTime(newt) 

solarDec  =  calcSunDeclination(newt) 

hourangle  =  calcHourAngleSunrise(latitude,  solarDec) 

delta  =  longitude  -  radToDeg(hourangle) 

timeDiff  =  4  *  delta 

timeUTC  =  720  +  timeDiff  -  eqtime 

calcSunriseUTC  =  timeUTC 

End  Function 


Function  calcSolNoonUTC(t,  longitude) 

'Calculate  the  the  Universal  Coordinated  Time  of  Solar  Noon  Declination  for  the  given  Day  and 
at  the  given  Location  on  Earth' 


Dim  newt  As  Double 
Dim  eqtime  As  Double 
Dim  solarNoonDec  As  Double 
Dim  solNoonUTC  As  Double 


'Declarate  the  new  "t"’ 

'Declarate  the  Equation  of  Time  in  Minutes' 
'Declarate  the  Solar  Noon  Declination' 

'Declarate  the  Solar  Noon  Declination  at  Universal 
Coordinated  Time' 


newt  =  calcTimeJulianCent(calcJDFromJulianCent(t)  +  0.5  +  longitude  /  360#) 


eqtime  =  calcEquationOfTime(newt) 
solarNoonDec  =  calcSunDeclination(newt) 
solNoonUTC  =  720  +  (longitude  *  4)  -  eqtime 


calc  SolNoonUTC  =  solNoonUTC 


End  Function 


Function  calcSunsetUTC(JD,  latitude,  longitude) 

'Calculate  the  the  Universal  Coordinated  Time  of  Sunset  for  the  given  Day  and  at  the  given 
Location  on  Earth' 


Dim  t  As  Double 
Dim  eqtime  As  Double 
Dim  solarDec  As  Double 
Dim  hourangle  As  Double 
Dim  delta  As  Double 

Dim  timeDiff  As  Double 
Dim  timeUTC  As  Double 


'Declarate  the  number  of  Julian  centuries  since  J2000.0' 
'Declarate  the  Equation  of  Time  in  Minutes' 

'Declarate  the  Solar  Declination  Angle' 

'Declarate  the  Local  Hour  Angle  in  Degrees  at  Sunset' 

'Declarate  the  Difference  between  Local  Longitude  and  Local 
Hour  Angle  in  Radians' 

'Declarate  4  times  Delta  Value' 

'Declarate  the  Universal  Coordinated  Time  of  Sunset  in  Minutes' 
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Dim  newt  As  Double 


Declarate  the  new  Mt"? 


t  =  calcTimeJulianCent(JD) 

'Calculates  Sunrise  and  approximate  Length  of  Day' 
eqtime  =  calcEquationOfTime(t) 
solarDec  =  calcSunDeclination(t) 
hourangle  =  calcHourAngleSunset(latitude,  solarDec) 

delta  =  longitude  -  radToDeg(hourangle) 

timeDiff  =  4  *  delta 

timeUTC  =  720  +  timeDiff  -  eqtime 

'First  pass  used  to  include  fractional  day  in  Gamma  Calculation' 

newt  =  calcTimeJulianCent(calcJDFromJulianCent(t)  +  timeUTC  /  1440#) 

eqtime  =  calcEquationOfTime(newt) 

solarDec  =  calcSunDeclination(newt) 

hourangle  =  calcHourAngleSunset(latitude,  solarDec) 

delta  =  longitude  -  radToDeg(hourangle) 

timeDiff  =  4  *  delta 

timeUTC  =  720  +  timeDiff  -  eqtime 

calcSunsetUTC  =  timeUTC 

End  Function 


Function  calcDuskUTC(JD,  latitude,  longitude,  solardepression) 

'Calculate  the  Universal  Coordinated  Time  of  Dusk  for  the  given  Day  at  the  given  Location  on 
Earth' 

Dim  t  As  Double 
Dim  eqtime  As  Double 
Dim  solarDec  As  Double 
Dim  hourangle  As  Double 
Dim  delta  As  Double 

Dim  timeDiff  As  Double 
Dim  timeUTC  As  Double 
Dim  newt  As  Double 

t  =  calcTimeJulianCent(JD) 

'First  calculates  Sunrise  and  approximate  Length  of  Day' 
eqtime  =  calcEquationOfTime(t) 
solarDec  =  calcSunDeclination(t) 
hourangle  =  calcHourAngleSunset(latitude,  solarDec) 

delta  =  longitude  -  radToDeg(hourangle) 


'Declarate  the  number  of  Julian  centuries  since  J2000.0' 
'Declarate  the  Equation  of  Time  in  Minutes' 

'Declarate  the  Solar  Declination  at  Sunset  in  Degrees' 

'Declarate  the  Local  Hour  Angle  of  Sunset  in  Degrees' 

'Declarate  the  Difference  between  the  Local  Longitude  and  Local 
Hour  Angle  in  Radians' 

'Declarate  4  times  Delta  Value' 

'Declarate  the  Universal  Coordinated  Time  of  Dusk  in  Minutes' 
'Declarate  the  new  ”t'" 
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timeDiff  =  4  *  delta 

timeUTC  =  720  +  timeDiff  -  eqtime 

'First  pass  used  to  include  fractional  day  in  Gamma  calculation' 

newt  =  calcTimeJulianCent(calcJDFromJulianCent(t)  +  timeUTC  /  1440#) 

eqtime  =  calcEquationOfTime(newt) 

solarDec  =  calcSunDeclination(newt) 

hourangle  =  calcHourAngleDusk(latitude,  solarDec,  solardepression) 

delta  =  longitude  -  radToDeg(hourangle) 

timeDiff  =  4  *  delta 

timeUTC  =  720  +  timeDiff  -  eqtime 

calcDuskUTC  =  timeUTC 

End  Function 


Function  dawn(lat,  Ion,  year,  month,  day,  timezone,  dlstime,  solardepression) 
'Calculate  the  Time  of  Dawn,  in  days,  for  a  given  Date  and  Location' 


Dim  longitude  As  Double 
Dim  latitude  As  Double 
Dim  JD  As  Double 
Dim  riseTimeGMT  As  Double 
Dim  riseTimeLST  As  Double 


'Declarate  the  Local  Longitude  in  Degrees' 

'Declarate  the  Local  Latitude  in  Degrees' 

'Declarate  the  Julian  Day' 

'Declarate  the  Sunrise  for  this  date' 

'Declarate  the  Sunrise  for  this  date  adjusting  for  the  Time 
Zone  and  Daylight  Savings,  in  Minutes' 


'Change  the  Sign  Convention  for  Longitude  from  Negative  to  Positive  in  Western  Hemisphere' 
longitude  =  Ion  *  -1 
latitude  =  lat 

If  (latitude  >  89.8)  Then  latitude  =  89.8 
If  (latitude  <  -89.8)  Then  latitude  =  -89.8 


JD  =  calcJD(year,  month,  day) 


riseTimeGMT  =  calcDawnUTC(JD,  latitude,  longitude,  solardepression) 


'Adjusting  for  Time  Zone  and  Daylight  Savings  Time  in  Minutes 
riseTimeLST  =  riseTimeGMT  +  (60  *  timezone)  +  (dlstime  *  60) 


'Convert  the  Dawn  from  Minutes  to  Days' 
dawn  =  riseTimeLST  /  1440 


End  Function 


Function  sunrise(lat,  Ion,  year,  month,  day,  timezone,  dlstime) 

'Calculate  the  Sunrise,  in  days,  for  a  given  Date  and  Location' 

Dim  longitude  As  Double  'Declarate  the  Local  Longitude  in  Degrees' 
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Dim  latitude  As  Double  'Declarate  the  Local  Latitude  in  Degrees' 

Dim  JD  As  Double  'Declarate  the  Julian  Day' 

Dim  riseTimeGMT  As  Double  'Declarate  the  Sunrise  time  for  this  date' 

Dim  riseTimeLST  As  Double  'Declarate  the  Sunrise  on  this  date  adjusting  for  the  Time 

Zone  and  Daylight  Savings,  in  Minutes' 

'Change  the  Sign  Convention  for  Longitude  from  Negative  to  Positive  in  Western  Hemisphere' 
longitude  =  Ion  *  -1 
latitude  =  lat 

If  (latitude  >  89.8)  Then  latitude  =  89.8 
If  (latitude  <  -89.8)  Then  latitude  =  -89.8 

JD  =  calcJD(year,  month,  day) 

riseTimeGMT  =  calcSunriseUTC(JD,  latitude,  longitude) 

'Adjusting  for  Time  Zone  and  Daylight  Savings  Time  in  Minutes 
riseTimeLST  =  riseTimeGMT  +  (60  *  timezone)  +  (dlstime  *  60) 

'Convert  the  Sunrise  from  Minutes  to  Days' 
sunrise  =  riseTimeLST  /  1440 

End  Function 


Function  solarnoon(lat,  Ion,  year,  month,  day,  timezone,  dlstime) 

'Calculate  the  Universal  Coordinated  Time  of  Solar  Noon,  in  days,  for  the  given  Day  at  the  given 
Location  on  Earth' 


Dim  longitude  As  Double 
Dim  latitude  As  Double 
Dim  JD  As  Double 
Dim  t  As  Double 
Dim  newt  As  Double 
Dim  eqtime  As  Double 
Dim  solarNoonDec  As  Double 
Dim  solNoonUTC  As  Double 


'Declarate  the  Local  Longitude  in  Degrees' 

’Declarate  the  Local  Latitude  in  Degrees' 

'Declarate  the  Julian  Day’ 

'Declarate  the  number  of  Julian  centuries  since  J2000.0’ 
'Declarate  the  new  "t"’ 

'Declarate  the  Equation  of  Time' 

'Declarate  the  Solar  Noon  Declination  Angle' 

'Declarate  the  Solar  Noon  on  Universal  Coordinated  Time' 


’Change  the  Sign  Convention  for  Longitude  from  Negative  to  Positive  in  Western  Hemisphere' 
longitude  =  Ion  *  -1 
latitude  =  lat 

If  (latitude  >  89.8)  Then  latitude  =  89.8 
If  (latitude  <  -89.8)  Then  latitude  =  -89.8 

JD  =  calcJD(year,  month,  day) 
t  =  calcTimeJulianCent(JD) 

newt  =  calcTimeJulianCent(calcJDFromJulianCent(t)  +  0.5  +  longitude  /  360#) 
eqtime  =  calcEquationOfTime(newt) 
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solarNoonDec  =  calcSunDeclination(newt) 
solNoonUTC  =  720  +  (longitude  *  4)  -  eqtime 

’Adjusting  for  Time  Zone  and  Daylight  Savings  Time  in  Minutes 
solarnoon  =  solNoonUTC  +  (60  *  timezone)  +  (dlstime  *  60) 

'Convert  the  Sunrise  from  Minutes  to  Days' 
solarnoon  =  solarnoon  /  1440 

End  Function 


Function  sunset(lat,  Ion,  year,  month,  day,  timezone,  dlstime) 
'Calculate  the  time  of  Sunset,  in  days,  for  the  given  Date  and  Location. 


Dim  longitude  As  Double 
Dim  latitude  As  Double 
Dim  JD  As  Double 
Dim  setTimeGMT  As  Double 
Dim  setTimeLST  As  Double 


'Declarate  the  Local  Longitude  in  Degrees' 

'Declarate  the  Local  Latitude  in  Degrees' 

'Declarate  the  Julian  Day' 

'Declarate  the  Sunset  time  of  this  Date’ 

'Declarate  the  the  Sunset  time  for  the  given  Date  and 
Location. 


'Change  the  Sign  Convention  for  Longitude  from  Negative  to  Positive  in  Western  Hemisphere' 
longitude  =  Ion  *  -1 
latitude  =  lat 

If  (latitude  >  89.8)  Then  latitude  =  89.8 
If  (latitude  <  -89.8)  Then  latitude  =  -89.8 


JD  =  calcJD(year,  month,  day) 

setTimeGMT  =  calcSunsetUTC(JD,  latitude,  longitude) 

Adjusting  for  Time  Zone  and  Daylight  Savings  Time  in  Minutes 
setTimeLST  =  setTimeGMT  +  (60  *  timezone)  +  (dlstime  *  60) 


'Convert  the  Sunrise  from  Minutes  to  Days' 
sunset  =  setTimeLST  /  1440 


End  Function 


Function  dusk(lat,  Ion,  year,  month,  day,  timezone,  dlstime,  solardepression) 
'Calculate  the  Time  of  Dusk,  in  days,  for  a  given  Date  and  Location’ 


Dim  longitude  As  Double 
Dim  latitude  As  Double 
Dim  JD  As  Double 
Dim  setTimeGMT  As  Double 
Dim  setTimeLST  As  Double 


'Declarate  the  Local  Longitude  in  Degrees' 

'Declarate  the  Local  Latitude  in  Degrees' 

'Declarate  the  Julian  Day' 

'Declarate  the  Sunset  for  this  day’ 

'Declarate  the  the  Time  of  Dusk,  in  days,  for  a  given  Date 
and  Location' 
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'Change  the  Sign  Convention  for  Longitude  from  Negative  to  Positive  in  Western  Hemisphere' 
longitude  =  Ion  *  -1 
latitude  =  lat 

If  (latitude  >  89.8)  Then  latitude  =  89.8 
If  (latitude  <  -89.8)  Then  latitude  =  -89.8 

JD  =  calcJD(year,  month,  day) 

setTimeGMT  =  calcDuskUTC(JD,  latitude,  longitude,  solardepression) 

'Adjusting  for  Time  Zone  and  Daylight  Savings  Time  in  Minutes 
setTimeLST  =  setTimeGMT  +  (60  *  timezone)  +  (dlstime  *  60) 

'Convert  the  Sunrise  from  Minutes  to  Days' 
dusk  =  setTimeLST  /  1440 

End  Function 


Function  solarazimuth(lat,  Ion,  year,  month,  day,  hours,  minutes,  seconds,  timezone,  dlstime) 
'Calculate  the  Solar  Azimuth  Angle  (Degrees  from  North)  for  a  given  Date  and  Location. 


Dim  longitude  As  Double 
Dim  latitude  As  Double 
Dim  zone  As  Double 
Dim  daySavings  As  Double 
Dim  hh  As  Double 
Dim  mm  As  Double 
Dim  ss  As  Double 
Dim  timenow  As  Double 
Dim  JD  As  Double 
Dim  t  As  Double 
Dim  R  As  Double 
Dim  alpha  As  Double 
Dim  theta  As  Double 
Dim  Etime  As  Double 
Dim  eqtime  As  Double 
Dim  solarDec  As  Double 
Dim  earthRadVec  As  Double 
Dim  solarTimeFix  As  Double 

Dim  trueSolarTime  As  Double 
Dim  hourangle  As  Double 

Dim  harad  As  Double 
Dim  csz  As  Double 
Dim  zenith  As  Double 
Dim  azDenom  As  Double 
Dim  azRad  As  Double 
Dim  azimuth  As  Double 


'Declarate  the  Local  Longitude  in  Degrees' 

'Declarate  the  Local  Latitude  in  Degrees' 

'Declarate  the  Time  Zone  in  Hours' 

'Declarate  the  Daylight  Savings  (0  =  no  or  1  =  yes)' 

'Declarate  the  Local  Hour' 

'Declarate  the  Local  Minute' 

'Declarate  the  Local  Second’ 

'Declarate  the  Time  in  Hours  (GMT)' 

'Declarate  the  Julian  Day' 

'Declarate  the  number  of  Julian  centuries  since  J2000.0’ 
'Declarate  the  Distance  to  the  Sun  in  Astronomical  Units  (UA) 
'Declarate  the  Right  Ascension  of  the  Sun  in  Degrees' 
'Declarate  the  Declination  of  the  Sun  in  Degrees' 

'Declarate  the  Equation  of  Time  in  Radians' 

'Declarate  the  Equation  of  Time  in  Degrees' 

'Declarate  the  Solar  Declination  Angle' 

'Declarate  the  Earth's  Radius  Vector' 

'Declarate  the  Corrected  Solar  Time,  in  Minutes,  using 
Equation  of  Time' 

'Declarate  the  True  Solar  Time  in  Minutes' 

'Declarate  the  Local  Hour  Angle  at  Sunrise,  given  a  Latitude 
and  Solar  Declination,  in  Degrees' 

'Declarate  the  Local  Hour  Angle  in  Radians' 

'Declarate  the  Solar  Zenith  in  Radians' 

'Declarate  the  Solar  Zenith  in  Degrees' 

'Declarate  the  Denominator’ 

'Declarate  the  Solar  Azimuth  in  Radians' 

'Declarate  the  Solar  Azimuth  Angle' 
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Dim  exoatmElevation  As  Double 


'Declarate  the  Differente  between  90  degrees  and  Zenith' 
'Declarate  the  Auxiliar  Variable  "stepl'" 

'Declarate  the  Auxiliar  Variable  Mstep2"? 

'Declarate  the  Auxiliar  Variable  ”step3"' 

'Declarate  the  Refraction  Correction' 

'Declarate  the  Tangent  of  the  Differente  between  90 
degrees  and  Zenith' 

'Declarate  the  Solar  Zenith  after  the  Refraction 
Correction' 


Dim  stepl  As  Double 
Dim  step2  As  Double 
Dim  step3  As  Double 


Dim  refractionCorrection  As  Double 
Dim  te  As  Double 


Dim  solarzen  As  Double 


'Change  the  Sign  Convention  for  Longitude  from  Negative  to  Positive  in  Western  Hemisphere' 
longitude  =  Ion  *  -1 
latitude  =  lat 

If  (latitude  >  89.8)  Then  latitude  =  89.8 
If  (latitude  <  -89.8)  Then  latitude  =  -89.8 

'Change  Time  Zone  to  Positive  hours  in  Western  Hemisphere' 
zone  =  timezone  *  -1 
daySavings  =  dlstime  *  60 
hh  =  hours  -  (daySavings  /  60) 
mm  =  minutes 
ss  =  seconds 

timenow  =  hh  +  mm  /  60  +  ss  /  3600  +  zone 

JD  =  calcJD(year,  month,  day) 
t  =  calcTimeJulianCent(JD  +  timenow  /  24#) 

R  =  calcSunRadVector(t) 
alpha  =  calcSunRtAscension(t) 
theta  =  calcSunDeclination(t) 

Etime  =  calcEquationOfTime(t) 

eqtime  =  Etime 
solarDec  =  theta 
earthRadVec  =  R 

solarTimeFix  =  eqtime  -  4#  *  longitude  +  60#  *  zone 
trueSolarTime  =  hh  *  60#  +  mm  +  ss  /  60#  +  solarTimeFix 

Do  While  (trueSolarTime  >  1440) 

trueSolarTime  =  trueSolarTime  -  1440 

Loop 

hourangle  =  trueSolarTime  /  4#  -  180# 

If  (hourangle  <  -180)  Then  hourangle  =  hourangle  +  360# 

harad  =  degToRad(hourangle) 
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csz  =  Sin(degToRad(latitude))  *  Sin(degToRad(solarDec))  +  Cos(degToRad(latitude))  * 
Cos(degToRad(solarDec))  *  Cos(harad) 

If  (csz  >1#)  Then 
csz  =  1# 

Elself  (csz  <  -1#)  Then 
csz  = -1# 

End  If 

zenith  =  radToDeg(Application.WorksheetFunction.Acos(csz)) 
azDenom  =  (Cos(degToRad(latitude))  *  Sin(degToRad(zenith))) 

If  (Abs(azDenom)  >  0.001)  Then 

azRad  =  ((Sin(degToRad(latitude))  *  Cos(degToRad(zenith)))  -  Sin(degToRad(solarDec)))  / 
azDenom 

If  (Abs(azRad)  >1#)  Then 
If  (azRad  <  0)  Then 
azRad  =  -1# 

Else 

azRad  =1# 

End  If 
End  If 

azimuth  =  180#  -  radToDeg( Application. WorksheetFunction.Acos(azRad)) 

If  (hourangle  >  0#)  Then 
azimuth  =  -azimuth 
End  If 
Else 

If  (latitude  >  0#)  Then 
azimuth  =  1 80# 

Else 

azimuth  =  0# 

End  If 
End  If 

If  (azimuth  <  0#)  Then 
azimuth  =  azimuth  +  360# 

End  If 

exoatmElevation  =  90#  -  zenith 

If  (exoatmElevation  >  85#)  Then 
refractionCorrection  =  0# 

Else 

te  =  Tan(degToRad(exoatmElevation)) 

If  (exoatmElevation  >  5#)  Then 

refractionCorrection  =  58.1  /  te  -  0.07  /  (te  *  te  *  te)  +  0.000086  /  (te  *  te  *  te  *  te  *  te) 
Elself  (exoatmElevation  >  -0.575)  Then 
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stepl  =  (-12.79  +  exoatmElevation  *  0.711) 
step2  =  (103.4  +  exoatmElevation  *  (stepl)) 
step3  =  (-518.2  +  exoatmElevation  *  (step2)) 
refractionCorrection  =  1735#  +  exoatmElevation  *  (step3) 
Else:  refractionCorrection  =  -20.774  /  te 
End  If 

refractionCorrection  =  refractionCorrection  /  3600# 

End  If 

solarzen  =  zenith  -  refractionCorrection 
solarazimuth  =  azimuth 
End  Function 


Function  solarelevation(lat,  Ion,  year,  month,  day,  hours,  minutes,  seconds,  timezone,  dlstime) 
'Calculate  the  Solar  Elevation  or  Altitude  (Degrees  from  North)  for  a  given  Date,  Time,  and 
Location' 


Dim  longitude  As  Double 
Dim  latitude  As  Double 
Dim  zone  As  Double 
Dim  daySavings  As  Double 
Dim  hh  As  Double 
Dim  mm  As  Double 
Dim  ss  As  Double 
Dim  timenow  As  Double 
Dim  JD  As  Double 
Dim  t  As  Double 
Dim  R  As  Double 

Dim  alpha  As  Double 
Dim  theta  As  Double 
Dim  Etime  As  Double 
Dim  eqtime  As  Double 
Dim  solarDec  As  Double 
Dim  earthRadVec  As  Double 
Dim  solarTimeFix  As  Double 

Dim  trueSolarTime  As  Double 
Dim  hourangle  As  Double 

Dim  harad  As  Double 

Dim  csz  As  Double 

Dim  zenith  As  Double 

Dim  azDenom  As  Double 

Dim  azRad  As  Double 

Dim  azimuth  As  Double 

Dim  exoatmElevation  As  Double 


'Declarate  the  Local  Longitude  in  Degrees' 

'Declarate  the  Local  Latitude  in  Degrees' 

'Declarate  the  Time  Zone  in  Hours' 

'Declarate  the  Daylight  Savings  (0  =  no  or  1  =  yes)' 
'Declarate  the  Local  Hour’ 

'Declarate  the  Local  Minute' 

'Declarate  the  Local  Second' 

'Declarate  the  Time  in  Hours  (GMT)' 

'Declarate  the  Julian  Day’ 

'Declarate  the  number  of  Julian  centuries  since  J2000.0' 
'Declarate  the  Distance  to  the  Sun  in  Astronomical  Units 
(UA) 

'Declarate  the  Right  Ascension  of  the  Sun  in  Degrees' 
'Declarate  the  Declination  of  the  Sun  in  Degrees' 

'Declarate  the  Equation  of  Time  in  Radians' 

'Declarate  the  Equation  of  Time  in  Degrees' 

'Declarate  the  Solar  Declination  Angle' 

'Declarate  the  Earth's  Radius  Vector' 

'Declarate  the  Corrected  Solar  Time,  in  Minutes,  using 
Equation  of  Time' 

'Declarate  the  True  Solar  Time  in  Minutes' 

'Declarate  the  Local  Hour  Angle  at  Sunrise,  given  a  Latitude 
and  Solar  Declination,  in  Degrees' 

'Declarate  the  Local  Hour  Angle  in  Radians' 

'Declarate  the  Solar  Zenith  in  Radians' 

'Declarate  the  Solar  Zenith  in  Degrees' 

'Declarate  the  Denominator 
'Declarate  the  Solar  Azimuth  in  Radians' 

'Declarate  the  Solar  Azimuth  Angle' 

'Declarate  the  Differente  between  90  degrees  and  Zenith' 
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Dim  stepl  As  Double 
Dim  step2  As  Double 
Dim  step3  As  Double 


'Declarate  the  Auxiliar  Variable  "stepl"1 
'Declarate  the  Auxiliar  Variable  "step2"f 
'Declarate  the  Auxiliar  Variable  "step3'" 

'Declarate  the  Refraction  Correction' 

'Declarate  the  Tangent  of  the  Different  between  90 
degrees  and  Zenith' 

'Declarate  the  Solar  Zenith  after  the  Refraction 
Correction' 


Dim  refractionCorrection  As  Double 
Dim  te  As  Double 


Dim  solarzen  As  Double 


'Change  the  Sign  Convention  for  Longitude  from  Negative  to  Positive  in  Western  Hemisphere' 
longitude  =  Ion  *  -1 
latitude  =  lat 

If  (latitude  >  89.8)  Then  latitude  =  89.8 
If  (latitude  <  -89.8)  Then  latitude  =  -89.8 

'Change  the  Time  Zone  to  Positive  hours  in  Western  Hemisphere 

zone  =  timezone  *  -1 

daySavings  =  dlstime  *  60 

hh  =  hours  -  (daySavings  /  60) 

mm  =  minutes 

ss  =  seconds 

timenow  =  hh  +  mm  /  60  +  ss  /  3600  +  zone 

JD  =  calcJD(year,  month,  day) 
t  =  calcTimeJulianCent(JD  +  timenow  /  24#) 

R  =  calcSunRadVector(t) 
alpha  =  calcSunRtAscension(t) 
theta  =  calcSunDeclination(t) 

Etime  =  calcEquationOfTime(t) 

eqtime  =  Etime 
solarDec  =  theta 
earthRadVec  =  R 

solarTimeFix  =  eqtime  -  4#  *  longitude  +  60#  *  zone 
trueSolarTime  =  hh  *  60#  +  mm  +  ss  /  60#  +  solarTimeFix 

Do  While  (trueSolarTime  >  1440) 
trueSolarTime  =  trueSolarTime  -  1440 
Loop 

hourangle  =  trueSolarTime  /  4#  -  180# 

If  (hourangle  <  -180)  Then  hourangle  =  hourangle  +  360# 

harad  =  degToRad(hourangle) 

csz  =  Sin(degToRad(latitude))  *  Sin(degToRad(solarDec))  +  Cos(degToRad(latitude))  * 
Cos(degToRad(solarDec))  *  Cos(harad) 
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If  (csz  >  1#)  Then 
csz  =  1# 

Elself  (csz  <  -1#)  Then 
csz  =  -1# 

End  If 

zenith  =  radToDeg(Application.WorksheetFunction.Acos(csz)) 

azDenom  =  (Cos(degToRad(latitude))  *  Sin(degToRad(zenith))) 

If  (Abs(azDenom)  >  0.001)  Then 

azRad  =  ((Sin(degToRad(latitude))  *  _ 

Cos(degToRad(zenith))) 

Sin(degToRad(solarDec)))  /  azDenom 
If  (Abs(azRad)  >1#)  Then 
If  (azRad  <  0)  Then 
azRad  =  -1# 

Else 

azRad  =  1# 

End  If 
End  If 

azimuth  =  180#  -  radToDeg(Application.WorksheetFunction.Acos(azRad)) 

If  (hourangle  >  0#)  Then 
azimuth  =  -azimuth 
End  If 
Else 

If  (latitude  >  0#)  Then 
azimuth  =  1 80# 

Else 

azimuth  =  0# 

End  If 
End  If 

If  (azimuth  <  0#)  Then 
azimuth  =  azimuth  +  360# 

End  If 

exoatmElevation  =  90#  -  zenith 

If  (exoatmElevation  >  85#)  Then 
refractionCorrection  =  0# 

Else 

te  =  Tan(degToRad(exoatmElevation)) 

If  (exoatmElevation  >  5#)  Then 

refractionCorrection  =  58.1  /  te  -  0.07  /  (te  *  te  *  te)  +  0.000086  /  (te  *  te  *  te  *  te  *  te) 
Elself  (exoatmElevation  >  -0.575)  Then 

stepl  =  (-12.79  +  exoatmElevation  *  0.711) 
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step2  =  (103.4  +  exoatmElevation  *  (stepl)) 
step3  =  (-518.2  +  exoatmElevation  *  (step2)) 
refractionCorrection  =  1735#  +  exoatmElevation  *  (step3) 
Else 

refractionCorrection  =  -20.774  /  te 
End  If 

refractionCorrection  =  refractionCorrection  /  3600# 

End  If 

solarzen  =  zenith  -  refractionCorrection 
solarelevation  =  90#  -  solarzen 
End  Function 


Sub  solarposition(lat,  Ion,  year,  month,  day,  hours,  minutes,  seconds,  timezone,  dlstime, 
solarazimuth,  solarelevation) 

'Calculate  the  Solar  Position  given  a  Date,  Time  and  Location' 


Dim 

Dim 

Dim 

Dim 

Dim 

Dim 

Dim 

Dim 

Dim 

Dim 

Dim 

Dim 

Dim 

Dim 

Dim 

Dim 

Dim 

Dim 

Dim 

Dim 

Dim 

Dim 

Dim 

Dim 

Dim 

Dim 

Dim 


longitude  As  Double 
latitude  As  Double 
zone  As  Double 
daySavings  As  Double 
hh  As  Double 
mm  As  Double 
ss  As  Double 
timenow  As  Double 
JD  As  Double 
t  As  Double 
R  As  Double 


'Declarate  the  Local  Longitude  in  Degrees' 

'Declarate  the  Local  Latitude  in  Degrees' 

'Declarate  the  Time  Zone  in  Hours' 

'Declarate  the  Daylight  Savings  (0  =  no  or  1  =  yes)' 
'Declarate  the  Local  Hour' 

'Declarate  the  Local  Minute' 

'Declarate  the  Local  Second' 

'Declarate  the  Time  in  Hours  (GMT)' 

'Declarate  the  Julian  Day' 

'Declarate  the  number  of  Julian  centuries  since  J2000.0' 
'Declarate  the  Distance  to  the  Sun  in  Astronomical  Units 


alpha  As  Double 
theta  As  Double 
Etime  As  Double 
eqtime  As  Double 
solarDec  As  Double 
earthRadVec  As  Double 
solarTimeFix  As  Double 

trueSolarTime  As  Double 
hourangle  As  Double 

harad  As  Double 
csz  As  Double 
zenith  As  Double 
azDenom  As  Double 
azRad  As  Double 
azimuth  As  Double 
exoatmElevation  As  Double 


(UA) 

'Declarate  the  Right  Ascension  of  the  Sun  in  Degrees' 
'Declarate  the  Declination  of  the  Sun  in  Degrees' 
'Declarate  the  Equation  of  Time  in  Radians' 

'Declarate  the  Equation  of  Time  in  Degrees' 

'Declarate  the  Solar  Declination  Angle' 

'Declarate  the  Earth's  Radius  Vector' 

'Declarate  the  Corrected  Solar  Time,  in  Minutes,  using 
Equation  of  Time' 

'Declarate  the  True  Solar  Time  in  Minutes' 

’Declarate  the  Local  Hour  Angle  at  Sunrise,  given  a 
Latitude  and  Solar  Declination,  in  Degrees' 

'Declarate  the  Local  Hour  Angle  in  Radians' 

'Declarate  the  Solar  Zenith  in  Radians' 

'Declarate  the  Solar  Zenith  in  Degrees' 

'Declarate  the  Denominator’ 

'Declarate  the  Solar  Azimuth  in  Radians' 

'Declarate  the  Solar  Azimuth  Angle' 

'Declarate  the  Differente  between  90  degrees  and  Zenith' 
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Dim  stepl  As  Double 
Dim  step2  As  Double 
Dim  step3  As  Double 


'Declarate  the  Auxiliar  Variable  "stepl'” 

'Declarate  the  Auxiliar  Variable  "step2"? 

'Declarate  the  Auxiliar  Variable  "step3"' 

'Declarate  the  Refraction  Correction' 

'Declarate  the  Tangent  of  the  Differente  between  90 
degrees  and  Zenith' 

'Declarate  the  Solar  Zenith  after  the  Refraction 
Correction' 


Dim  refractionCorrection  As  Double 
Dim  te  As  Double 


Dim  solarzen  As  Double 


'Change  the  Sign  Convention  for  Longitude  from  Negative  to  Positive  in  Western  Hemisphere' 
longitude  =  Ion  *  -1 
latitude  =  lat 

If  (latitude  >  89.8)  Then  latitude  =  89.8 
If  (latitude  <  -89.8)  Then  latitude  =  -89.8 

'Change  the  Time  Zone  to  Positive  hours  in  Western  Hemisphere 

zone  =  timezone  *  -1 

daySavings  =  dlstime  *  60 

hh  =  hours  -  (daySavings  /  60) 

mm  =  minutes 

ss  =  seconds 

timenow  =  hh  +  mm  /  60  +  ss  /  3600  +  zone 

JD  =  calcJD(year,  month,  day) 
t  =  calcTimeJulianCent(JD  +  timenow  /  24#) 

R  =  calcSunRadVector(t) 
alpha  =  calcSunRtAscension(t) 
theta  =  calcSunDeclination(t) 

Etime  =  calcEquationOfTime(t) 

eqtime  =  Etime 
solarDec  =  theta 
earthRadVec  =  R 

solarTimeFix  =  eqtime  -  4#  *  longitude  +  60#  *  zone 
trueSolarTime  =  hh  *  60#  +  mm  +  ss  /  60#  +  solarTimeFix 

Do  While  (trueSolarTime  >  1440) 
trueSolarTime  =  trueSolarTime  -  1440 
Loop 

hourangle  =  trueSolarTime  /  4#  -  180# 

If  (hourangle  <  -180)  Then  hourangle  =  hourangle  +  360# 

harad  =  degToRad(hourangle) 

csz  =  Sin(degToRad(latitude))  *  Sin(degToRad(solarDec))  +  Cos(degToRad(latitude))  * 
Cos(degToRad(solarDec))  *  Cos(harad) 
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If  (csz  >  1#)  Then 
csz  =  1# 

Elself  (csz  <  -1#)  Then 
csz  =  -1# 

End  If 

zenith  =  radToDeg(Application.WorksheetFunction.Acos(csz)) 

azDenom  =  (Cos(degToRad(latitude))  *  Sin(degToRad(zenith))) 

If  (Abs(azDenom)  >  0.001)  Then 

azRad  =  ((Sin(degToRad(latitude))  *  _ 

Cos(degToRad(zenith))) 

Sin(degToRad(solarDec)))  /  azDenom 
If  (Abs(azRad)  >1#)  Then 
If  (azRad  <  0)  Then 
azRad  =  -1# 

Else 

azRad  =  1# 

End  If 
End  If 

azimuth  =  180#  -  radToDeg(Application.WorksheetFunction.Acos(azRad)) 

If  (hourangle  >  0#)  Then 
azimuth  =  -azimuth 
End  If 
Else 

If  (latitude  >  0#)  Then 
azimuth  =  1 80# 

Else 

azimuth  =  0# 

End  If 
End  If 

If  (azimuth  <  0#)  Then 
azimuth  =  azimuth  +  360# 

End  If 

exoatmElevation  =  90#  -  zenith 

If  (exoatmElevation  >  85#)  Then 
refractionCorrection  =  0# 

Else 

te  =  Tan(degToRad(exoatmElevation)) 

If  (exoatmElevation  >  5#)  Then 

refractionCorrection  =  58.1  /  te  -  0.07  /  (te  *  te  *  te)  +  0.000086  /  (te  *  te  *  te  *  te  *  te) 
Elself  (exoatmElevation  >  -0.575)  Then 
stepl  =  (-12.79  +  exoatmElevation  *  0.71 1) 
step2  =  (103.4  +  exoatmElevation  *  (stepl)) 
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step3  =  (-518.2  +  exoatmElevation  *  (step2)) 
refractionCorrection  =  1735#  +  exoatmElevation  *  (step3) 
Else 

refractionCorrection  =  -20.774  /  te 
End  If 

refractionCorrection  =  refractionCorrection  /  3600# 

End  If 

solarzen  =  zenith  -  refractionCorrection 

solarazimuth  =  azimuth 
solarelevation  =  90#  -  solarzen 

End  Sub 


Function  travel_risk(start_point  As  Double,  end_point  As  Integer,  JD  As  Double)  As  Double 
'Calculate  the  total  risk  during  travel  on  the  Military  network' 


Dim  i  As  Long 

Dim  j  As  Integer 

Dim  hour  As  Double 

Dim  min  As  Double 

Dim  sec  As  Double 

Dim  day  As  Double 

Dim  month  As  Double 

Dim  year  As  Double 

Dim  LY  As  Integer 

Dim  TZ  As  Integer 

Dim  D  As  Integer 

Dim  avg  lat  As  Double 

Dim  avg_lon  As  Double 

Dim  shadow_proj  As  Double 

Dim  diff_AZij_solAZ  As  Double 

Dim  SolAz  As  Double 
Dim  AZij  As  Double 
Dim  Sol_alt  As  Double 
Dim  SolAz_d  As  Double 
Dim  Real_shadow  As  Double 
Dim  virtual  shadow  As  Double 
Dim  tree  As  Double 
Dim  boat_speed  As  Double 
Dim  boat  width  As  Double 
Dim  river_speed  As  Double 
Dim  ratio  earth  As  Double 
Dim  pi  As  Double 
Dim  e  As  Double 
Dim  point(65,  10)  As  Double 


'Declarate  the  Start  Point' 

'Declarate  the  End  Point' 

'Declarate  Hours' 

'Declarate  Minutes' 

'Declarate  Seconds' 

'Declatate  Day' 

'Declarate  Month’ 

'Declarate  Year' 

'Declarate  Leap  Year' 

'Declarate  Time -Zone' 

'Declarate  Daylight  Savings' 

'Declarate  Average  Latitude  in  Radians' 

'Declarate  Average  Longitude  in  Radians' 

'Declarate  Solar  Projection’ 

'Declarate  the  Difference  between  Path  Azimuth  and  Solar 
Azimuth' 

'Declarate  Solar  Azimuth  Angle’ 

'Declarate  Azimuth  between  Points  "i"  and  "j"' 

'Declarate  Solar  Altitude’ 

'Declarate  Solar  Azimuth  in  Degrees' 

'Declarate  Real  Shadow  Lenght' 

'Declarate  Virtual  Shadow  Lenght’ 

'Declarate  Riverine  Forest  Height' 

'Declarate  Boat  Speed’ 

'Declarate  Boat  Width' 

'Declarate  River  Speed’ 

'Declarate  Ratio  Earth' 

'Declarate  the  Mathematical  Constant  Pi' 

'Declarate  the  Mathematical  Constant  Euler's  number  e' 
’Point[i][l]  Latitute  degrees  at  point  ”i"' 

'Point[i][2]  Latitute  minutes  at  point  "i"’ 
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?Point[i][3]  Latitute  seconds  at  point  "i?" 
?Point[i][4]  Longitude  degrees  at  point  "i"’ 
?Point[i][5]  Longitude  minutes  at  point  "i"' 
?Point[i][6]  Longitude  seconds  at  point  "i"’ 
?Point[i][7]  Solar  Azimuth  at  point  "i"' 
?Point[i][8]  Solar  Altitude  at  point  "i?" 
?Point[i][9]  Geocentric  Latitude  of  point  "i"’ 
'Point[i][10]  Longitude  of  Point  "i"' 


Dim  lat_d  As  Double 
Dim  lat_m  As  Double 
Dim  lat_s  As  Double 
Dim  lat  As  Double 
Dim  lon_d  As  Double 
Dim  lon_m  As  Double 
Dim  lon_s  As  Double 
Dim  Ion  As  Double 
Dim  Rad_lon_i  As  Double 
Dim  Rad_lonj  As  Double 
Dim  Rad_lat_i  As  Double 
Dim  Rad_lat ;_j  As  Double 
Dim  DD  lat  i  As  Double 
Dim  DD  lon  i  As  Double 
Dim  DD  lat _j  As  Double 
Dim  DD  lon _j  As  Double 
Dim  LC  As  Double 
Dim  RL  As  Double 
Dim  COSS  As  Double 
Dim  ACOSS  As  Double 
Dim  DISTS  As  Double 
Dim  travel_time  As  Double 
e  =  Exp(l) 


'Declarate  Degrees  of  Latitude  of  Point  "i?M 
'Declarate  Minutes  of  Latitude  of  Point  "i1" 

'Declarate  Seconds  of  Latitude  of  Point  "i'M 
'Declarate  Latitude  of  Point  "i"  in  Radians' 

'Declarate  Degrees  of  Latitude  of  Point  "i"' 

'Declarate  Minutes  of  Latitude  of  Point  "i"' 

'Declarate  Seconds  of  Latitude  of  Point  "i"' 

'Declarate  Longitude  of  Point  "i”  in  Radians' 

'Declarate  Longitude  of  Point  ”i"  in  Radians' 

'Declarate  Longitude  of  Point  "j"  in  Radians' 

'Declarate  Latitude  of  Point  "i"  in  Radians' 

'Declarate  Latitude  of  Point  "j"  in  Radians' 

'Declarate  Latitude  of  Point  ”i"  in  Degrees' 

'Declarate  Longitude  of  Point  "i"  in  Degrees' 

'Declarate  Latitude  of  Point  "j"  in  Degrees' 

'Declarate  Longitude  of  Point  "j”  in  Degrees' 

'Declarate  Longitude  Correction' 

'Declarate  Reference  Longitude' 

'Declarate  Cosine  of  Distance  between  Points  "i"  and  "j"' 
'Declarate  ArcCosine  of  Distance  between  Points  "i"  and  "j"' 
'Declarate  Cosine  of  Distance  between  Points  "i"  and  "j"' 
'Declarate  travel  time  between  points  "i"  and  "j"' 


i  =  start_point  'Declarate  the  start  point' 
j  =  end_point  'Declarate  the  end  point' 


'Get  variables' 

TZ  =  Worksheets("Input").Cells(2,  6).Value 
LY  =  Worksheets("Input").Cells(6,  4).Value 
D  =  Worksheets("Input").Cells(2,  8).Value 
hour  =  Worksheets("Input").Cells(3,  6).Value 
min  =  Worksheets("Input").Cells(4,  6).Value 
sec  =  Worksheets("Input").Cells(5,  6).Value 
day  =  Worksheets("Input").Cells(3,  4).Value 
month  =  Worksheets("Input").Cells(4,  4).Value 
year  =  Worksheets("Input").Cells(5,  4).Value 
boat_width  =  Worksheets("Input").Cells(7,  4).Value 
boat_speed  =  Worksheets("Input").Cells(8,  4).Value 
ratio  earth  =  Worksheets("Input").Cells(9,  4).Value 


'Get  Time -Zone' 

'Get  Leap  Year' 

'Get  Daylight  Savings' 
'Get  Hour' 

'Get  Minute' 

'Get  Second' 

'Get  Day' 

'Get  Month' 

'Get  Year' 

'Get  Boat  width' 

'Get  Boat  Speed' 

'Get  the  Ratio  of  Earth' 
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pi  =  Application.WorksheetFunction.piQ 


'Get  the  Mathematical  constant  Pi' 


'Calculate  Geocentric  Latitude  and  Longitude  of  Start  Point  in  Radians' 

point(i,  1)  =  Worksheets("Input").Cells(12  +  i,  4).Value  'Get  Degrees  of  Latitude  of  Point  "i"' 
lat_d  =  point(i,  1) 

point(i,  2)  =  Worksheets("Input").Cells(12  +  i,  5).Value  'Get  Minutes  of  Latitude  of  Point  "i"' 
latm  =  point(i,  2) 

point(i,  3)  =  Worksheets("Input").Cells(12  +  i,  6).Value  'Get  Seconds  of  Latitude  of  Point  "i"’ 
lat  s  =  point(i,  3) 

point(i,  4)  =  Worksheets("Input").Cells(12  +  i,  7).Value  'Get  Degrees  of  Longitude  of  Point  "i'M 
lon  d  =  point(i,  4) 

point(i,  5)  =  Worksheets("Input").Cells(12  +  i,  8).Value  'Get  Minutes  of  Longitude  of  Point  "i"f 
lon_m  =  point(i,  5) 

point(i,  6)  =  Worksheets("Input").Cells(12  +  i,  9).Value  'Get  Seconds  of  Longitude  of  Point  "i"' 
Ion  s  =  point(i,  6) 

'Calculate  Latitude  of  Start  Point  in  Degrees' 

DD  lat  i  =  lat  d  +  lat  m  /  60  +  lat  s  /  3600 

'Calculate  Latitude  of  Start  Point  in  Radians' 

Radlati  =  WorksheetFunction.Radians(DDlati) 

'Calculate  Longitude  of  Start  Point  in  Degrees' 

DD  lon  i  =  lon  d  +  lon  m  /  60  +  Ion  s  /  3600 

'Calculate  Longitude  of  Start  Point  in  Radians' 

Radloni  =  WorksheetFunction.Radians(DDloni) 
point(i,  10)  =  Rad_lon_i 

'Calculate  Geocentric  Latitude  of  Start  Point  in  Radians' 
lat  =  Atn((l  -  (1  /  2983)  A  2)  *  Tan(Rad_lat_i)) 
point(i,  9)  =  lat 

'Calculate  Geocentric  Latitude  and  Longitude  of  End  Point  in  Radians' 

point(j,  1)  =  Worksheets("Input").Cells(12  +  j,  4).Value  Get  Degrees  of  Latitude  of  Point  "i"' 

lat_d  =  point(j,  1) 

point(j,  2)  =  Worksheets("Input").Cells(12  +  j,  5).Value  'Get  Minutes  of  Latitude  of  Point  "i"' 
lat  m  =  point(j,  2) 

point(j,  3)  =  Worksheets("Input").Cells(12  +  j,  6).Value  'Get  Seconds  of  Latitude  of  Point  "i"' 
lat_s  =  point(j,  3) 

point(j,  4)  =  Worksheets("Input").Cells(12  +  j,  7).Value  'Get  Degrees  of  Longitude  of  Point  "i"' 
lon_d  =  point(j,  4) 

point(j,  5)  =  Worksheets("Input").Cells(12  +  j,  8).Value  'Get  Minutes  of  Longitude  of  Point  "i"' 
lon  m  =  point(j,  5) 

point(j,  6)  =  Worksheets("Input").Cells(12  +  j,  9).Value  'Get  Seconds  of  Longitude  of  Point  "i"' 
lon_s  =  point(j,  6) 
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'Calculate  Latitude  of  End  Point  in  Degrees' 

DD  latj  =  lat  d  +  lat  m  /  60  +  lat  s  /  3600 

'Calculate  Latitude  of  End  Point  in  Radians' 

Rad_latJ  =  WorksheetFunction.Radians(DD_latJ) 

'Calculate  Longitute  of  End  Point  in  Degrees' 

DD_lon_j  =  lon_d  +  lon_m  /  60  +  lon_s  /  3600 

'Calculate  Longitude  of  End  Point  in  Radians' 

Radlonj  =  WorksheetFunction.Radians(DDlonJ) 
point(j,  10)  =  Rad_lonj 

'Calculate  Geocentric  Latitude  of  End  Point  in  Radians' 
lat  =  Atn((l  -  (1  /  2983)  A  2)  *  Tan(RadJatJ)) 
point(j,  9)  =  lat 

'Calculate  Cosine  of  Distance  between  Points  i  and  j' 

COSS  =  Sin(point(j,  9))  *  Sin(point(i,  9))  +  Cos(point(j,  9))  *  Cos(point(i,  9))  *  Cos(point(j, 

10)  -  point(i,  10)) 

'Calculate  Distance  between  Points  i  and  j' 

ACOSS  =  Application.WorksheetFunction.Acos(COSS) 

'Calculate  Distance  between  Points  i  and  j  in  miles' 

DISTS  =  ACOSS  *  ratioearth 

'Calculate  Path  Azimuth  (between  Points  i  and  j)' 

AZij  =  Application. WorksheetFunction.Asin(Sin(pi  /  2  -  point(j,  9))  *  Sin(point(j,  10)  -  point(i, 
10))  /  Sin(ACOSS)) 

If  AZij  <  0  Then  AZij  =  AZij  +  2  *  pi 


'Calculate  Path  Azimuth  in  Degrees' 
Dim  AZij_d  As  Double 
AZij_d  =  AZij  *  180 /pi 


'Get  information  about  hight  riverine  forest' 

tree  =  Worksheets("Input").Cells(12  +  i,  12  +  j). Value  'Get  Riverine  Forest  Height' 

'Calculate  Average  Latitude  between  Points  i  and  j' 
avg  lat  =  (Rad  lat  i  +  Rad  latj)  /  2 

'Calculate  Solar  Azimuth' 

SolAz  d  =  solarazimuth(DD_latj,  DDlonj,  year,  month,  day,  hour,  min,  sec,  TZ,  D) 
SolAz  =  SolAz_d  *  pi  /  180 

'Calculate  Solar  Altitude' 

Dim  Sol_alt_d  As  Double 

Sol  alt  d  =  solarelevation(DD_latj,  DD  lonj,  year,  month,  day,  hour,  min,  sec,  TZ,  D) 
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Sol_alt  =  Sol_alt_d  *  pi  / 180 

'Calculate  the  Difference  between  Solar  Azimuth  and  Path  Azimuth' 

'Both  Solar  and  Path  Azimuths  are  in  first  quadrant' 

If  AZij  <=  pi  /  2  And  SolAz  <=  pi  /  2  Then 
If  SolAz  >  AZij  Then 
diffAZijsolAZ  =  SolAz  -  AZij 
Else:  diff_AZij_solAZ  =  AZij  -  SolAz 
End  If 
End  If 

'Path  Azimuth  is  in  the  first  quadrant  and  Solar  Azimuth  is  in  the  second  one 
If  AZij  <=  pi  /  2  And  SolAz  >=  pi  /  2  And  SolAz  <=  pi  Then 
If  (SolAz  -  pi  /  2)  <  AZij  Then 
diff_AZij_solAZ  =  SolAz  -  AZij 
Else:  diff_AZij_solAZ  =  pi  -  SolAz  +  AZij 
End  If 
End  If 

'Path  Azimuth  is  in  the  first  quadrant  and  Solar  Azimuth  is  in  the  third  one' 
If  AZij  <=  pi  /  2  And  SolAz  >=  pi  And  SolAz  <=  pi  Then 
If  (SolAz  -  pi)  >  AZij  Then 
diff  AZij  solAZ  =  SolAz  -  pi  -  AZij 
Else:  diff  AZij  solAZ  =  AZij  -  (SolAz  -  pi) 

End  If 
End  If 

'Path  Azimuth  is  in  the  first  quadrant  and  Solar  Azimuth  is  in  the  fourth  one' 
If  AZij  <=  pi  /  2  And  SolAz  >=  3  *  pi  /  2  Then 
If  (SolAz  -  3  *  pi  /  2)  <  AZij  Then 
diff  AZij  solAZ  =  SolAz  -  AZij  -  pi 
Else:  diff_AZij_solAZ  =  2  *  pi  -  SolAz  +  AZij 
End  If 
End  If 

'Solar  Azimuth  is  in  the  first  quadrant  and  Path  Azimuth  is  in  the  second  one 
If  AZij  >=  pi  /  2  And  AZij  <=  pi  And  SolAz  <=  pi  /  2  Then 
If  (AZij  -  pi  /  2)  <  SolAz  Then 
diff_AZij_solAZ  =  pi  -  AZij  +  SolAz 
Else:  diff  AZij  solAZ  =  AZij  -  SolAz 
End  If 
End  If 

'Solar  Azimuth  is  in  the  first  quadrant  and  Path  Azimuth  is  in  the  third  one' 
If  AZij  >=  pi  And  AZij  <=  3  *  pi  /  2  And  SolAz  <=  pi  /  2  Then 
If  (AZij  -  pi)  <  SolAz  Then 
diff_AZij_solAZ  =  AZij  -  pi  -  SolAz 
Else:  diff_AZij_solAZ  =  SolAz  -  (AZij  -  pi) 
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End  If 
End  If 

'Both  Solar  Azimuth  and  Path  Azimuth  are  in  the  second  quadrant' 

If  AZij  >=  pi  /  2  And  AZij  <=  pi  And  SolAz  >=  pi  /  2  And  SolAz  <=  pi  Then 
If  AZij  <  SolAz  Then 
diffAZijsolAZ  =  SolAz  -  AZij 
Else:  diff  AZij  solAZ  =  AZij  -  SolAz 
End  If 
End  If 

'Solar  Azimuth  is  in  second  quadrant  and  Path  Azimuth  is  in  the  third  one' 

If  AZij  >=  pi  And  AZij  <=  3  *  pi  /  2  And  SolAz  >=  pi  /  2  And  SolAz  <=  pi  Then 
If  AZij  >  SolAz  Then 
diff  AZij  solAZ  =  AZij  -  SolAz 
Else:  diff_AZij_solAZ  =  SolAz  -  (AZij  -  pi) 

End  If 
End  If 

'Both  Solar  and  Path  Azimuths  are  in  the  third  quadrant' 

If  AZij  >=  pi  And  AZij  <=  3  *  pi  /  2  And  SolAz  >=  pi  And  SolAz  <=  3  *  pi  /  2  Then 
If  SolAz  >  AZij  Then 
diff_AZij_solAZ  =  SolAz  -  AZij 
Else:  diff  AZij  solAZ  =  AZij  -  SolAz 
End  If 
End  If 

'Path  Azimuths  is  in  the  third  quadrant  and  Solar  Azimuth  is  in  the  fourth  one' 

If  AZij  <=  3  *  pi  /  2  And  AZij  >=  pi  And  SolAz  >=  3  *  pi  /  2  Then 
If  SolAz  >  AZij  Then 
diff  AZij  solAZ  =  SolAz  -  AZij 
Else:  diff  AZij  solAZ  =  (2  *  pi  -  SolAz)  +  (AZij  -  pi) 

End  If 
End  If 

'Both  Path  and  Solar  Azimuths  are  in  the  fourth  quadrant' 

If  AZij  >=  3  *  pi  /  2  And  SolAz  >=  3  *  pi  /  2  Then 
If  SolAz  >  AZij  Then 
diff_AZij_solAZ  =  SolAz  -  AZij 
Else:  diff_AZij_solAZ  =  AZij  -  SolAz 
End  If 
End  If 

'Calculate  Risk  between  Points  i  and  j' 

'Calculate  Shadow  length' 

If  Solalt  <  pi  /  2  Or  Solalt  >  pi  /  2  Then 
shadow_proj  =  tree  /  Tan(Sol_alt) 

If  Sol_alt  =  pi  /  2  Then  shadow_proj  =  0 
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End  If 


'Calculate  Virtual  Shadow  Lenght' 

virtual_shadow  =  shadow_proj  *  Sin(diff_AZij_solAZ) 

'Calculate  Real  Shadow  Lenght' 

Real_shadow  =  Application.WorksheetFunction.Max(0,  virtual_shadow) 

If  Realshadow  /  boatwidth  >  1  Then 
travelrisk  =  0 

Elself  Real_shadow  /  boat_width  >  0  Then 
travel_risk  =  Real_shadow  /  boat_width 
Elself  Real  shadow  /  boat  width  =  0  Then 
travel_risk  =  1 
End  If 

'Current  time' 

Dim  curtime  As  Double 
Dim  DS  As  Double 

DS  =  Sheets(”input").Cells(2,  8) 

TZ  =  Sheets("input").Cells(2,  6) 

cur  time  =  (JD  -  Int(JD))  *  24  +  12  -  DS  +  TZ 
If  cur_time  >  24  Then  cur_time  =  cur_time  -  24 

If  cur_time  <  Sheets("input").Cells(2,  10)  Then 
travelrisk  =  0 

Elself  cur  time  >  Sheets("input").Cells(4,  10)  Then 
travelrisk  =  0 
End  If 

End  Function 


Sub  Dijkstra_time() 

Dim  Rad_lat_i  As  Double 
Dim  Rad  lat  j  As  Double 
Dim  Rad_lon_i  As  Double 
Dim  Rad  lonj  As  Double 
Dim  DD  lat  i  As  Double 
Dim  DD  lon  i  As  Double 
Dim  DD  lat  J  As  Double 
Dim  DD  lonj  As  Double 
Dim  COSS  As  Double 
Dim  ACOSS  As  Double 
Dim  DISTS  As  Double 
Dim  lat_d  As  Double 
Dim  lat  m  As  Double 


'Declarate  Latitude  of  Point  "i"  in  Radians' 

'Declarate  Latitude  of  Point  "j”  in  Radians' 

'Declarate  Longitude  of  Point  ”i”  in  Radians' 

'Declarate  Longitude  of  Point  ”j”  in  Radians' 

'Declarate  Latitude  of  Point  "i"  in  degrees' 

'Declarate  Longitude  of  Point  ”i"  in  Degrees' 

'Declarate  Latitude  of  Point  "j"  in  Degrees' 

'Declarate  Longitude  of  Point  "j"  in  Degrees' 

'Declarate  Cosine  of  Distance  between  Points  ”i"  and  "j"' 
'Declarate  ArcCosine  of  Distance  between  Points  "i"and  "j”' 
'Declarate  Cosine  of  Distance  between  Points  "i”  and  "j"' 
'Declarate  Degrees  of  Latitude  of  Point  ”i'" 

'Declarate  Minutes  of  Latitude  of  Point  ”i'" 
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Dim  lat_s  As  Double 
Dim  lat  As  Double 
Dim  lon_d  As  Double 
Dim  lon_m  As  Double 
Dim  lon_s  As  Double 
Dim  Ion  As  Double 
Dim  point()  As  Double 


Dim  boat_speed  As  Double 
Dim  boat_width  As  Double 
Dim  river_speed  As  Double 
Dim  ratio_earth  As  Double 
Dim  pi  As  Double 
Dim  hour  As  Double 
Dim  min  As  Double 
Dim  sec  As  Double 
Dim  day  As  Double 
Dim  month  As  Double 
Dim  year  As  Double 
Dim  JD  As  Double 
Dim  start  time  As  Double 
Dim  UT  As  Double 
Dim  DS  As  Integer 
Dim  t  As  Double 
Dim  numnodes  As  Long 
Dim  Q()  As  Double 
Dim  adj()  As  Integer 
Dim  adj_risk()  As  Double 
Dim  adj_time()  As  Double 
Dim  previous()  As  Long 
Dim  risk()  As  Double 
Dim  time()  As  Double 
Dim  old_time()  As  Double 
Dim  gap  As  Integer 
Dim  k  As  Integer 
Dim  m  As  Integer 
Dim  i  As  Integer 
Dim  j  As  Integer 
Dim  infty  As  Double 
Dim  count  As  Double 


'Declarate  Seconds  of  Latitude  of  Point  "i,M 
'Declarate  Latitude  of  Point  "i"  in  Radians' 
'Declarate  Degrees  of  Latitude  of  Point  "i"? 
'Declarate  Minutes  of  Latitude  of  Point  "i”' 
'Declarate  Seconds  of  Latitude  of  Point  "i"' 
'Declarate  Longitude  of  Point  "i"in  Radians' 
’Point[i][l]  Latitute  degrees  at  point  "i"’ 
?Point[i][2]  Latitute  minutes  at  point  "i"' 
'Point[i][3]  Latitute  seconds  at  point  ”i”' 
'Point[i][4]  Longitude  degrees  at  point  "i”' 
?Point[i][5]  Longitude  minutes  at  point  "i”’ 
fPoint[i][6]  Longitude  seconds  at  point  "i”' 
’Point[i][7]  Solar  Azimuth  at  point  "i'” 

?Point[i][8]  Solar  Altitude  at  point  i' 

?Point[i][9]  Geocentric  Latitude  of  point  i' 
'Point[i][10]  Longitude  of  point  i' 

'Declarate  the  Boat  Speed' 

'Declarate  the  Boat  Width' 

'Declarate  the  River  Speed' 

'Declarate  the  Ratio  Earth' 

'Declarate  the  Mathematical  Constant  Pi' 
'Declarate  Hours' 

'Declarate  Minutes' 

'Declarate  Seconds' 

'Declatate  Day' 

'Declarate  Month' 

'Declarate  Year' 

'Declarate  the  Julian  Date' 

'Declarate  Start  Time' 

'Declarate  Universal  Time' 

'Declarate  Daylight  Saving' 

'Declarate  Auxiliar  Variable  “t”' 

'Declarate  the  Total  Number  of  Nodes' 
'Declarate  the  Set  Q' 

'Declarate  the  Adjacent  Node' 

'Declarate  the  Risk  of  Adjacent  Node' 
'Declarate  the  Time  between  Adjacent  Nodes' 
'Declarate  the  Previous  Node' 

'Declarate  the  Risk' 

'Declarate  the  Time' 

'Declarate  the  Old  Time' 

'Declarate  the  Gap' 

'Declarate  the  Auxiliar  Variable  "k"' 
'Declarate  the  Auxiliar  Variable  "m"’ 
'Declarate  the  Node  "i"' 

'Declarate  the  Node  "j'" 

'Declarate  the  Infinity  value' 

'Declarate  the  Counter' 
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Dim  TZ  As  Double 
Dim  Source  As  Integer 
Dim  tempi  As  Double 
Dim  temp2  As  Double 
Dim  temptime  As  Double 
Dim  alt_risk  As  Double 
Dim  sink  As  Integer 
Dim  rownum  As  Integer 
Dim  firstnode  As  Integer 
Dim  secondnode  As  Integer 
Dim  temp  As  Integer 

'Count  Number  of  Nodes' 
numnodes  =13 
While  Sheets("Input").Cells(numnodes,  2)  <>  "" 
numnodes  =  numnodes  +  1 
Wend 

numnodes  =  numnodes  -  13 

Sheets("output").Range("Al:E65000").ClearContents 


'Declarate  the  Time-Zone' 

'Declarate  the  Source  Node' 

'Declarate  the  Auxiliar  Variable  "temp  1"' 
'Declarate  the  Auxiliar  Variable  "temp  2"' 
'Declarate  the  Auxiliar  Variable  "temptime"' 
'Declarate  the  Auxiliar  Variable  "Alt  Risk'" 
'Declarate  the  Auxiliar  Variable  "sink"' 
'Declarate  Row  Number' 

'Declarate  the  First  Node  to  be  Visited' 
'Declarate  the  Second  Node  to  be  Visited' 
'Declarate  Auxiliar  Variable  "temp"' 


Sheets("output").Range("J2:n5000").ClearContents 


'Redimension  variables' 

ReDim  Q(numnodes,  2)  As  Double 

ReDim  adj  (numnodes,  numnodes)  As  Integer 

ReDim  adj_risk(numnodes,  numnodes)  As  Double 

ReDim  adj_time(numnodes,  numnodes)  As  Double 

ReDim  previous(numnodes)  As  Long 
ReDim  risk(numnodes)  As  Double 
ReDim  time(numnodes)  As  Double 
ReDim  old_time(numnodes)  As  Double 
ReDim  point(numnodes,  10) 


'Redimension  of  number  of  nodes' 
'Redimension  of  the  Adjacent  Node' 
'Redimension  of  the  Risk  of  Adjacent 
Node' 

'Redimension  of  the  Time  related  to  the 
Adjacent  Node' 

'Redimension  of  the  Previous  Node' 
'Redimension  of  the  Risk' 

'Redimension  of  the  Time' 
'Redimension  of  the  Old  Time' 
'Redimension  of  the  Point  Number' 


'Get  variables' 

TZ  =  Worksheets("Input").Cells(2,  6).Value 
DS  =  Worksheets("Input").Cells(2,  8).Value 
hour  =  Worksheets("Input").Cells(3,  6).Value 
min  =  Worksheets("Input").Cells(4,  6).Value 
sec  =  Worksheets("Input").Cells(5,  6).Value 
day  =  Worksheets("Input").Cells(3,  4).Value 
month  =  Worksheets("Input").Cells(4,  4).Value 
year  =  Worksheets("Input").Cells(5,  4).Value 
boat  width  =  Worksheets("Input").Cells(7,  4).Value 
boat  speed  =  Worksheets("Input").Cells(8,  4).Value 
ratio  earth  =  Worksheets("Input").Cells(9,  4).Value 
pi  =  Application.WorksheetFunction.piQ 


'Get  Time  Zone' 

'Get  Daylight  Savings' 
'Get  Hour' 

'Get  Minute' 

'Get  Second' 

'Get  Day' 

'Get  Month' 

'Get  Year' 

'Get  the  Boat  Width' 

'Get  the  Boat  Speed' 

'Get  the  Ratio  of  the  Earth' 
'Get  the  pi  value' 
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For  i  =  1  To  numnodes 

'Calculate  Geocentric  Latitude  and  Longitude  of  Start  Point  in  Radians' 

point(i,  1)  =  Worksheets("Input").Cells(12  +  i,  4).Value  ’Get  Degrees  of  Latitude  of  Point  "i"' 

lat_d  =  point(i,  1) 

point(i,  2)  =  Worksheets("Input").Cells(12  +  i,  5).Value  ’Get  Minutes  of  Latitude  of  Point  "i"' 
latm  =  point(i,  2) 

point(i,  3)  =  Worksheets("Input”).Cells(12  +  i,  6).Value  'Get  Seconds  of  Latitude  of  Point  "i”' 
lat_s  =  point(i,  3) 

point(i,  4)  =  Worksheets("Input").Cells(12  +  i,  7).Value  'Get  Degrees  of  Longitude  of  Point  "i" 
lon  d  =  point(i,  4) 

point(i,  5)  =  Worksheets("Input").Cells(12  +  i,  8).Value  'Get  Minutes  of  Longitude  of  Point  "i" 
lon_m  =  point(i,  5) 

point(i,  6)  =  Worksheets("Input").Cells(12  +  i,  9).Value  'Get  Seconds  of  Longitude  of  Point  "i" 
Ion  s  =  point(i,  6) 

'Calculate  Latitude  of  Start  Point  in  Degrees' 

DD  lat  i  =  lat  d  +  lat  m  /  60  +  lat  s  /  3600 

'Calculate  Latitude  of  Start  Point  in  Radians' 

Radlati  =  WorksheetFunction.Radians(DDlati) 

'Calculate  Longitude  of  Start  Point  in  Degrees' 

DD  lon  i  =  lon  d  +  lon  m  /  60  +  Ion  s  /  3600 

'Calculate  Longitude  of  Start  Point  in  Radians' 

Radloni  =  WorksheetFunction.Radians(DDloni) 
point(i,  10)  =  Rad_lon_i 

'Calculate  Geocentric  Latitude  of  Start  Point  in  Radians' 
lat  =  Atn((l  -  (1  /  2983)  A  2)  *  Tan(Rad_lat_i)) 
point(i,  9)  =  lat 
Next  i 

infty  =  9999999 

’Get  adjacency  matrix  and  put  in  adj(IJ) 
count  =  81 

While  Sheets("Input").Cells(count,  2)  <>  "" 

adj(Sheets("Input").Cells(count,  2),  Sheets("Input").Cells(count,  3))  =  1 
adj(Sheets("Input").Cells(count,  3),  Sheets("Input").Cells(count,  2))  =  1 
count  =  count  +  1 
Wend 

'Initionalize  Variables' 

'Calculate  Universal  Time’ 

If  (hour  +  min  /  60  +  sec  /  3600  +  DS  -  TZ)  <  0  Then 
UT  =  (hour  +  min  /  60  +  sec  /  3600  +  DS  -  TZ)  +  24 
Elself  (hour  +  min  /  60  +  sec  /  3600  +  DS  -  TZ)  >  24  Then 
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UT  =  (hour  +  min  /  60  +  sec  /  3600  +  DS  -  TZ)  -  24 
Else 

UT  =  (hour  +  min  /  60  +  sec  /  3600  +  DS  -  TZ) 

End  If 


'Calculate  the  Julian  Date' 

JD  =  367  *  year  -  Int(7  *  (year  +  Int((month  +  9)  /  12))  /  4)  -  Int(3  *  (Int((year  +  (month  -  9)  /  7) 
/ 100)  +  1)  /  4)  +  Int(275  *  month  /  9)  +  day  +  1721028.5  +  UT  /  24 

start_time  =  JD 

Source  =  InputBox("Enter  number  of  node  to  start  at",  "Start  Node",  1)  'Get  start  node 

Sheets("input").Cells(2,  10)  =  (sunrise(radToDeg(point(Source,  9)),  radToDeg(point(Source, 
10)),  year,  month,  day,  TZ,  DS))  *  24 

Sheets("input").Cells(3,  10)  =  (solamoon(radToDeg(point(Source,  9)),  radToDeg(point(Source, 
10)),  year,  month,  day,  TZ,  DS))  *  24 

Sheets("input").Cells(4,  10)  =  (sunset(radToDeg(point(Source,  9)),  radToDeg(point(Source, 
10)),  year,  month,  day,  TZ,  DS))  *  24 


For  i  =  0  To  numnodes 

Q(i>  l)  =  i 

Q(i,  2)  =  infty 
risk(i)  =  infty 
previous(i)  =  infty 
time(i)  =  infty 
old_time(i)  =  start_time 
time(Source)  =  starttime 
time(0)  =  start_time 
Next  i 


’Initialization  Procedures' 

'Put  all  nodes  in  Q' 

'Unknown  distance  function  from  source  to  vertice  "v"’ 

'Total  risk  of  shortest  path  to  node  "i"' 

'Node  from  which  traveled  to  node  "i"  in  optimal  path  from  source' 
'Travel  time  to  node  "i"' 

’Travel  time  to  node  "i"' 

'Set  source  start  time  to  begingin  time' 

’Set  time  at  begining  to  begin  time' 


'Start  the  Algorithm' 

risk(Source)  =  0  'Distance  from  source  to  source’ 

Q(Source,  2)  =  0 
Q(0,  2)  =  -1  *  infty 
previous(Source)  =  0 


Do  'The  main  loop' 

'Sort  Q  by  risk' 
gap  =  (numnodes  -  1)  /  2 
While  (gap  >  0) 

For  k  =  gap  To  numnodes 
For  m  =  k  -  gap  To  0  Step  -gap 
If  (Q(m,  2)  >  Q(m  +  gap,  2))  Then 
temp2  =  Q(m,  2) 
tempi  =  Q(m,  1) 

Q(m,  2)  =  Q(m  +  gap,  2) 

Q(m  +  gap,  2)  =  temp2 
Q(m,  1)  =  Q(m  +  gap,  1) 
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Q(m  +  gap,  1)  =  tempi 
End  If 
Next  m 
Next  k 
gap  =  gap  /  2 
Wend 

If  Q(l,  2)  =  infty  Then 

Exit  Do  'All  remaining  vertices  are  inaccessible  from  source' 

End  If 

Q(l,  2)  =  infty  'Remove  "u"  from  Q  sets  dist  to  infmfity  so  wont  appear  in  sort' 

'Removes  “u”  from  adjacency  matrix  for  next  operation 
For  i  =  1  To  numnodes 
adj(Q(l,l),  i)  =  0 
Next  i 

'For  each  neighbor  ”v”  of  ”u”:  where  ”v"  has  not  yet  been  removed  from  Q.' 

For  i  =  1  To  numnodes 
If  adj(i,  Q(l,  1))=  1  Then 

'Compute  risk  for  time(q(l,l))  based  on  coming  from  node  Q(l,l)  to  node  ”i’". 

'Calculate  Cosine  of  Distance  between  Points  "i"  and  "j”' 

COSS  =  Sin(point(i,  9))  *  Sin(point(Q(l,  1),  9))  +  Cos(point(i,  9))  *  Cos(point(Q(l,  1), 

9))  *  Cos(point(i,  10)  -  point(Q(l,  1),  10)) 

'Calculate  Distance  between  Points  ”i”  and  "j”' 

ACOSS  =  Application.WorksheetFunctionAcos(COSS) 

'Calculate  Distance  between  Points  "i”  and  "j"  in  miles' 

DISTS  =  ACOSS  *  ratio  earth 

temptime  =  time(Q(l,  1))  +  (DISTS  /  boat_speed)  /  24 
alt_risk  =  risk(Q(l,  1))  +  travel_risk(Q(l,  1),  i,  temptime) 

If  temptime  <  time(i)  Then  'Relax  (u,v,a)' 

oldtime(i)  =  temptime 
time(i)  =  temptime 
previous(i)  =  Q(l,  1) 

For  j  =  1  To  numnodes 

If  Q(j,  1)  =  i  Then  Q(j,  2)  =  time(i) 

Next  j 

risk(i)  =  alt_risk 
End  If 
End  If 
Next  i 

Loop  Until  numnodes  <  0  '  loop  until  all  nodes  have  been  evaluated 
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'Get  apth  risk  and  time  traveled' 

Sheets("output").Cells(l,  1)  =  "Fastest  Path" 

Sheets("output").Cells(l,  2)  =  "Risk" 

Sheets("output").Cells(l,  3)  =  "Time" 

Sheets("output").Cells(l,  4)  =  "Survivability" 

Sheets("output").Cells(l,  5)  =  "Weighted  Survivability" 

sink  =  InputBox("Enter  node  number  for  drop  off  of  supplies",  "Supply  Drop  Off  Node", 
numnodes) 

rownum  =  2 

Sheets("output").Cells(rownum,  1)  =  sink 

Sheets(" output ").Cells(rownum,  2)  =  risk(sink)  -  risk(previous(sink)) 

Sheets(" output ").Cells(rownum,  3)  =  (time(sink)  -  time(previous(sink)))  *  24 
Sheets("output").Cells(rownum,  4)  =  1  -  Sheets("output").Cells(rownum,  2) 
Sheets("output").Cells(rownum,  5)  =  Sheets("output").Cells(rownum,  3)  * 
Sheets("output").Cells(rownum,  4) 

Sheets("output").Cells(rownum,  5)  =  Sheets("output").Cells(rownum,  3)  * 
Sheets("output").Cells(rownum,  4) 

rownum  =  rownum  +  1 
count  =  previous(sink) 

While  count  >  0 

Sheets("output").Cells(rownum,  1)  =  count 

Sheets("output").Cells(rownum,  2)  =  risk(count)  -  risk(previous(count)) 
Sheets("output").Cells(rownum,  3)  =  (time(count)  -  time(previous(count)))  *  24 
Sheets("output").Cells(rownum,  4)  =  1  -  Sheets("output").Cells(rownum,  2) 
Sheets("output").Cells(rownum,  5)  =  Sheets("output").Cells(rownum,  3)  * 
Sheets("output").Cells(rownum,  4) 
rownum  =  rownum  +  1 
count  =  previous(count) 

Wend 

Sheets("output").Cells(rownum  -  1,  2)  =  0 
Sheets("output").Cells(rownum  -  1,  4)  =  1 

rownum  =  rownum  +  1 

Sheets("output").Cells(rownum,  2)  =  "=sum(b2:"  &  "b"  &  rownum  -  1  &  ")" 
Sheets("output").Cells(rownum,  3)  =  "=sum(c2:"  &  "c"  &  rownum  -  1  &  ")" 
Sheets("output").Cells(rownum,  4)  =  "=product(d2:"  &  "d"  &  rownum  -  1  &  ")" 
Sheets("output").Cells(rownum,  5)  =  "=sum(e2:"  &  "e"  &  rownum  -  3  &  ")/sum(c2:"  &  "c"  & 
rownum  -  3  &  ")" 

rownum  =  rownum  +  2 

'Get  shortest  risk  paths' 

Call  Dijkstra(Source,  sink,  rownum) 

Call  getalldata 
End  Sub 
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Sub  Dijkstra(start_node  As  Integer,  end  node  As  Integer,  rownum2  As  Integer) 


Dim  Rad_lat_i  As  Double 
Dim  Rad_lat ;_j  As  Double 
Dim  Rad_lon_i  As  Double 
Dim  Rad_lonj  As  Double 
Dim  DD  lat  i  As  Double 
Dim  DD  lon  i  As  Double 
Dim  DD  lat  J  As  Double 
Dim  DD  lonj  As  Double 
Dim  COSS  As  Double 
Dim  ACOSS  As  Double 
Dim  DISTS  As  Double 
Dim  lat_d  As  Double 
Dim  lat_m  As  Double 
Dim  lat_s  As  Double 
Dim  lat  As  Double 
Dim  lon_d  As  Double 
Dim  lon_m  As  Double 
Dim  lon_s  As  Double 
Dim  Ion  As  Double 
Dim  pointQ  As  Double 


Dim  boat_speed  As  Double 
Dim  boat_width  As  Double 
Dim  river_speed  As  Double 
Dim  ratioearth  As  Double 
Dim  pi  As  Double 
Dim  hour  As  Double 
Dim  min  As  Double 
Dim  sec  As  Double 
Dim  day  As  Double 
Dim  month  As  Double 
Dim  year  As  Double 
Dim  JD  As  Double 
Dim  start_time  As  Double 
Dim  UT  As  Double 
Dim  DS  As  Integer 
Dim  t  As  Double 
Dim  numnodes  As  Long 


'Declarate  Latitude  of  Point  "i"  in  Radians' 

'Declarate  Latitude  of  Point  "j”  in  Radians' 

'Declarate  Longitude  of  Point  "i”  in  Radians' 

'Declarate  Longitude  of  Point  "j"  in  Radians' 

'Declarate  Latitude  of  Point  "i"  in  degrees' 

'Declarate  Longitude  of  Point  "i"  in  Degrees' 

'Declarate  Latitude  of  Point  ”j"  in  Degrees' 

'Declarate  Longitude  of  Point  "j"  in  Degrees' 

'Declarate  Cosine  of  Distance  between  Points  "i"  and  "j"’ 
'Declarate  ArcCosine  of  Distance  between  Points  "i”  and  "j"' 
'Declarate  Cosine  of  Distance  between  Points  "i"  and  "j"' 
'Declarate  Degrees  of  Latitude  of  Point  "i"’ 

'Declarate  Minutes  of  Latitude  of  Point  "i"' 

'Declarate  Seconds  of  Latitude  of  Point  "i"' 

'Declarate  Latitude  of  Point  "i"  in  Radians' 

'Declarate  Degrees  of  Latitude  of  Point  "i"' 

'Declarate  Minutes  of  Latitude  of  Point  "i"' 

'Declarate  Seconds  of  Latitude  of  Point  "i"' 

'Declarate  Longitude  of  Point  "i"  in  Radians' 

'Point[i][l]  Latitute  degrees  at  point  "i"' 

'Point[i][2]  Latitute  minutes  at  point  "i"' 

'Point[i][3]  Latitute  seconds  at  point  "i"' 

'Point[i][4]  Longitude  degrees  at  point  "i"' 

'Point[i][5]  Longitude  minutes  at  point  "i"' 

'Point[i][6]  Longitude  seconds  at  point  "i"' 

'Point[i][7]  Solar  Azimuth  at  point  "i"' 

'Point[i][8]  Solar  Altitude  at  point  "i"’ 

’Point[i][9]  Geocentric  Latitude  of  point  "i"' 

'Point[i][10]  Longitude  of  point  "i"' 

'Declarate  Boat  Speed' 

'Declarate  Boat  Width’ 

'Declarate  River  Speed' 

'Declarate  Ratio  Earth' 

'Declarate  the  Mathematical  Constant  Pi' 

'Declarate  Hours' 

'Declarate  Minutes' 

'Declarate  Seconds' 

'Declatate  Day' 

'Declarate  Month' 

'Declarate  Year' 

'Declarate  Julian  Date' 

'Declarate  Start  Time' 

'Declarate  Universal  Time 
'Declarate  Daylight  Saving' 

'Declarate  the  Auxiliar  Variable  “t”' 

'Declarate  the  Total  Number  of  Nodes' 
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Dim  Q()  As  Double 

Dim  adj()  As  Integer 

Dim  adj_risk()  As  Double 

Dim  adj_time()  As  Double 

Dim  previous()  As  Long 

Dim  risk()  As  Double 

Dim  time()  As  Double 

Dim  gap  As  Integer 

Dim  k  As  Integer 

Dim  m  As  Integer 

Dim  i  As  Integer 

Dim  j  As  Integer 

Dim  infty  As  Double 

Dim  count  As  Double 

Dim  TZ  As  Double 

Dim  Source  As  Integer 

Dim  tempi  As  Double 

Dim  temp2  As  Double 

Dim  temptime  As  Double 

Dim  alt_risk  As  Double 

Dim  sink  As  Integer 

Dim  rownum  As  Integer 

Dim  flrstnode  As  Integer 

Dim  secondnode  As  Integer 

Dim  temp  As  Integer 

Dim  timestep  As  Integer 

Dim  startsum  As  Integer 

Dim  minutes  between  As  Integer 


'Declarate  the  set  Q' 

'Declarate  the  Adjacent  Node' 

'Declarate  the  Risk  of  Adjacent  Node' 

'Declarate  the  Time  between  Adjacent  Nodes' 
'Declarate  the  Previous  node' 

'Declarate  the  Risk' 

'Declarate  the  time' 

'Declarate  the  Gap' 

'Declarate  the  Auxiliar  Variable  "k"' 

'Declarate  the  Auxiliar  Variable  "m" 

'Declarate  the  node  "i"' 

'Declarate  the  node  "j"' 

'Declarate  the  infinity  number' 

'Declarate  the  counter' 

'Declarate  the  Time-Zone' 

'Declarate  the  Source  Node' 

'Declarate  the  Auxiliar  Variable  "temp  1'" 
'Declarate  the  Auxiliar  Variable  "temp  2"' 
'Declarate  the  Auxiliar  Variable  “temp'" 
'Declarate  the  Auxiliar  Variable  "alt_risk" 
'Declarate  the  Auxiliar  Variable  "sink" 

'Declarate  the  Row  Number' 

'Declarate  the  First  Node  to  be  visited' 

'Declarate  the  Second  Node  to  be  visited' 
'Declarate  the  Auxiliar  Variable  "temp"' 
'Declarate  the  Auxiliar  Variable  "timestep"' 
'Declarate  the  Auxiliar  Variable  "startsum"' 
'Declare  the  Auxiliar  Variable  "minutes  between 


ff! 


minutes  between  =  CInt(InputBox("Enter  the  number  of  minutes  between  each  test  path", 
"Time  between  Paths",  15)) 

'Count  Number  of  Nodes' 
numnodes  =13 

While  Sheets("Input").Cells(numnodes,  2)  <>  "" 

numnodes  =  numnodes  +  1 

Wend 

numnodes  =  numnodes  -  13 


'Redimension  variables' 

ReDim  Q(numnodes,  2)  As  Double 

ReDim  adj  (numnodes,  numnodes)  As  Integer 

ReDim  adj_risk(numnodes,  numnodes)  As  Double 

ReDim  adj_time(numnodes,  numnodes)  As  Double 

ReDim  previous(numnodes)  As  Long 
ReDim  risk(numnodes)  As  Double 
ReDim  time(numnodes)  As  Double 
ReDim  point(numnodes,  10) 


'Redimension  of  the  number  of  Nodes' 
'Redimension  of  the  Adjacent  Node' 
'Redimension  of  the  Risk  of  the 
Adjancent  Nodes' 

'Redimension  of  the  Time  between 
Adjacent  Nodes' 

'Redimension  of  the  Previous  Node' 
'Redimension  of  the  Risk' 
'Redimension  of  the  Time' 
'Redimension  of  the  Node' 
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rownum  =  rownum2 


While  (Sheets("Input").Cells(3,  6)  +  (Sheets("Input").Cells(4,  6)  /  60))  <= 
Sheets("INput").Cells(5,  8) 

'Get  variables' 

TZ  =  Worksheets("Input").Cells(2,  6).Value 
DS  =  Worksheets("Input").Cells(2,  8).Value 
hour  =  Worksheets("Input").Cells(3,  6).Value 
min  =  Worksheets("Input").Cells(4,  6).Value 
sec  =  Worksheets("Input").Cells(5,  6).Value 
day  =  Worksheets("Input").Cells(3, 4).Value 
month  =  Worksheets("Input").Cells(4,  4).Value 
year  =  Worksheets("Input").Cells(5,  4).Value 
boat_width  =  Worksheets("Input").Cells(7, 4).Value 
boat_speed  =  Worksheets("Input").Cells(8, 4).Value 
ratio  earth  =  Worksheets("Input").Cells(9,  4).Value 
pi  =  Application.WorksheetFunction.pi() 

For  i  =  1  To  numnodes 

'Calculate  Geocentric  Latitude  and  Longitude  of  Start  Point  in  Radians' 

point(i,  1)  =  Worksheets("Input").Cells(12  +  i,  4).Value  'Get  Degrees  of  Latitude  of  Point  "i"' 

latd  =  point(i,  1) 

point(i,  2)  =  Worksheets("Input").Cells(12  +  i,  5).Value  Get  Minutes  of  Latitude  of  Point  "i"' 
lat_m  =  point(i,  2) 

point(i,  3)  =  Worksheets("Input").Cells(12  +  i,  6).Value  'Get  Seconds  of  Latitude  of  Point  "i"' 
lat_s  =  point(i,  3) 

point(i,  4)  =  Worksheets("Input").Cells(12  +  i,  7).Value  'Get  Degrees  of  Longitude  of  Point  "i"' 
lon  d  =  point(i,  4) 

point(i,  5)  =  Worksheets("Input").Cells(12  +  i,  8).Value  'Get  Minutes  of  Longitude  of  Point  "i"' 
lon_m  =  point(i,  5) 

point(i,  6)  =  Worksheets("Input").Cells(12  +  i,  9).Value  'Get  Seconds  of  Longitude  of  Point  "i"' 
Ion  s  =  point(i,  6) 

'Calculate  Latitude  of  Start  Point  in  Degrees' 

DD  lat  i  =  lat  d  +  latm  /  60  +  lat  s  /  3600 

'Calculate  Latitude  of  Start  Point  in  Radians' 

Radlati  =  WorksheetFunction.Radians(DDlati) 

'Calculate  Longitude  of  Start  Point  in  Degrees' 

DD_lon_i  =  lon_d  +  lon_m  /  60  +  lon_s  /  3600 

'Calculate  Longitude  of  Start  Point  in  Radians' 

Radloni  =  WorksheetFunction.Radians(DDloni) 
point(i,  10)  =  Rad_lon_i 


'Get  Time  Zone' 

'Get  Daylight  Savings' 

'Get  Hours' 

'Get  Minutes' 

'Get  Seconds' 

'Get  Days' 

'Get  Months' 

'Get  Years' 

'Get  the  Boath  Width' 

'Get  the  Boat  Speed' 

'Get  the  Ratio  of  the  Earth' 

'Get  the  the  Mathematical  constant  Pi ' 
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'Calculate  Geocentric  Latitude  of  Start  Point  in  Radians' 
lat  =  Atn((l  -  (1  /  2983)  A  2)  *  Tan(Rad_lat_i)) 
point(i,  9)  =  lat 
Next  i 

infty  =  9999999 

'Get  adjacency  matrix  and  put  in  adj(I,j) ' 
count  =  81 

While  Sheets("Input").Cells(count,  2)  <>  "" 

adj(Sheets("Input").Cells(count,  2),  Sheets("Input").Cells(count,  3))  =  1 
adj(Sheets("Input").Cells(count,  3),  Sheets("Input").Cells(count,  2))  =  1 
count  =  count  +  1 
Wend 

'Initionalize  Variables' 


'Calculate  Universal  Time  (UT)' 

If  (hour  +  min  /  60  +  sec  /  3600  +  DS  -  TZ)  <  0  Then 
UT  =  (hour  +  min  /  60  +  sec  /  3600  +  DS  -  TZ)  +  24 
Elself  (hour  +  min  /  60  +  sec  /  3600  +  DS  -  TZ)  >  24  Then 
UT  =  (hour  +  min  /  60  +  sec  /  3600  +  DS  -  TZ)  -  24 
Else 

UT  =  (hour  +  min  /  60  +  sec  /  3600  +  DS  -  TZ) 

End  If 


'Calculate  the  Julian  Date' 

JD  =  367  *  year  -  Int(7  *  (year  +  Int((month  +  9)  /  12))  /  4)  -  Int(3  *  (Int((year  +  (month  -  9)  / 
7)  / 100)  +  1)  /  4)  +  Int(275  *  month  /  9)  +  day  +  1721028.5  +  UT  /  24 

starttime  =  JD 

Source  =  start  node  'get  start  node 


For  i  =  0  To  numnodes 

Q(i,  l)  =  i 

Q(i,  2)  =  infty 
risk(i)  =  infty 
previous(i)  =  infty 

time(i)  =  starttime 
Next  i 

’Start  the  Algorithm 
risk(Source)  =  0 
Q(Source,  2)  =  0 
Q(0,  2)  =  -1  *  infty 
previous(Source)  =  0 

Do 


'Initialization  Procedures' 

'Put  all  nodes  in  Q’ 

’Unknown  distance  function  from  source  to  "v"' 

’Total  risk  of  shortest  path  to  node  ”i” 

'Node  from  which  traveled  to  node  ”i”  in  optimal  path  from 
source' 

'Travel  time  to  node  "i"' 


'Distance  from  source  to  source’ 


'The  main  loop’ 
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'Sort  Q  by  risk' 

gap  =  (numnodes  -  1)  /  2 
While  (gap  >  0) 

For  k  =  gap  To  numnodes 
For  m  =  k  -  gap  To  0  Step  -gap 
If  (Q(m,  2)  >  Q(m  +  gap,  2))  Then 
temp2  =  Q(m,  2) 
tempi  =  Q(m,  1) 

Q(m,  2)  =  Q(m  +  gap,  2) 

Q(m  +  gap,  2)  =  temp2 
Q(m,  1)  =  Q(m  +  gap,  1) 

Q(m  +  gap,  1)  =  tempi 
End  If 
Next  m 
Next  k 
gap  =  gap  /  2 
Wend 

If  Q(l,  2)  =  infty  Then 

Exit  Do  'All  remaining  vertices  are  inaccessible  from  source' 

End  If 

Q(l,  2)  =  infty  'Remove  u  from  Q  sets  dist  to  infmfity  so  wont  appear  in  sort' 

’Removes  "u"  from  adjacency  matrix  for  next  operation’ 

For  i  =  1  To  numnodes 
adj(Q(l,l),  i)  =  0 
Next  i 

’For  each  neighbor  ”v”  of  "u”:  where  "v"  has  not  yet  been  removed  from  Q’ 

For  i  =  1  To  numnodes 
If  adj(i,  Q(l,  1))  =  1  Then 

'Compute  risk  for  time(q(l,l))  based  on  coming  from  node  Q(l,l)  to  node  "i"' 

’Calculate  Cosine  of  Distance  between  Points  "i"  and  "j"’ 

COSS  =  Sin(point(i,  9))  *  Sin(point(Q(l,  1),  9))  +  Cos(point(i,  9))  *  Cos(point(Q(l,  1),  9))  * 
Cos(point(i,  10)  -point(Q(l,  1),  10)) 

'Calculate  Distance  between  Points  "i”  and  "j”’ 

ACOSS  =  Application.WorksheetFunctionAcos(COSS) 

'Calculate  Distance  between  Points  "i"  and  "j"  in  miles' 

DISTS  =  ACOSS  *  ratio  earth 

temptime  =  time(Q(l,  1))  +  (DISTS  /  boat_speed)  /  24 

altrisk  =  risk(Q(l,  1))  +  travel_risk(Q(l,  1),  i,  temptime) 

If  alt  risk  <  risk(i)  Then  'Relax  (u,v,a)' 
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risk(i)  =  alt_risk 
previous(i)  =  Q(l,  1) 

For  j  =  1  To  numnodes 

If  Q(j,  1)  =  i  Then  Q(j,  2)  =  risk(i) 

Next  j 

time(i)  =  temptime 

Elself  alt_risk  =  risk(i)  And  temptime  <  time(i)  Then 
risk(i)  =  alt_risk 
previous(i)  =  Q(l,  1) 

For  j  =  1  To  numnodes 

If  Q(j,  1)  =  i  Then  Q(j,  2)  =  risk(i) 

Next  j 

time(i)  =  temptime 

End  If 
End  If 
Next  i 

Loop  Until  numnodes  <  0  'Loop  until  all  nodes  have  been  evaluated' 

'Get  apth  risk  and  time  traveled' 

Sheets(" output ").Cells(rownum,  1)  =  "Path  "  &  Sheets("Input").Cells(3,  6)  &  & 

Format(Sheets("Input").Cells(4,  6),  "00") 

Sheets("output").Cells(rownum,  2)  =  "Risk" 

Sheets(" output ").Cells(rownum,  3)  =  "Time" 

Sheets("  output")  Uells(rownum,  4)  =  "Survivability" 

Sheets("output").Cells(rownum,  5)  =  "Weighted  Survivability" 
Sheets("output").Cells(rownum,  5)  =  "Weighted  Survivability" 

sink  =  end_node 
rownum  =  rownum  +  1 
startsum  =  rownum 

Sheets("  output  ").Cells(rownum,  1)  =  sink 

Sheets("output").Cells(rownum,  2)  =  risk(sink)  -  risk(previous(sink)) 

Sheets("  output")  Uells(rownum,  3)  =  (time(sink)  -  time(previous(sink)))  *  24 
Sheets(" output ").Cells(rownum,  4)  =  1  -  Sheets("output").Cells(rownum,  2) 
Sheets("output").Cells(rownum,  5)  =  Sheets("output").Cells(rownum,  3)  * 
Sheets("output").Cells(rownum,  4) 

rownum  =  rownum  +  1 
count  =  previous(sink) 

While  count  >  0 

Sheets("output").Cells(rownum,  1)  =  count 

Sheets("output").Cells(rownum,  2)  =  risk(count)  -  risk(previous(count)) 
Sheets("output").Cells(rownum,  3)  =  (time(count)  -  time(previous(count)))  *  24 
Sheets("output").Cells(rownum,  4)  =  1  -  Sheets("output").Cells(rownum,  2) 
Sheets("output").Cells(rownum,  5)  =  Sheets("output").Cells(rownum,  3)  * 
Sheets("output").Cells(rownum,  4) 
rownum  =  rownum  +  1 
count  =  previous(count) 
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Wend 

Sheets("output").Cells(rownum  -  1,  2)  =  0 
Sheets("output").Cells(rownum  -  1,  4)  =  1 

rownum  =  rownum  +  1 

Sheets("output").Cells(rownum,  2)  =  "=sum(b"  &  startsum  &  &  "b"  &  rownum  -  1  &  ")" 
Sheets("output").Cells(rownum,  3)  =  "=sum(c"  &  startsum  &  &  "c"  &  rownum  -  1  &  ")" 
Sheets("output").Cells(rownum,  4)  =  "=product(d"  &  startsum  &  &  "d"  &  rownum  -  1  & 

fiyi 

Sheets("output").Cells(rownum,  5)  =  "=sum(e"  &  startsum  &  &  "e"  &  rownum  -  3  & 

")/sum(c"  &  startsum  &  &  "c"  &  rownum  -  3  &  ")" 

rownum  =  rownum  +  2 

If  Sheets("Input").Cells(4,  6)  =  60  -  minutes_between  Then 
Sheets("Input").Cells(4,  6)  =  0 

Sheets("Input").Cells(3,  6)  =  Sheets("Input").Cells(3,  6)  +  1 
Else 

Sheets("Input").Cells(4,  6)  =  Sheets("Input").Cells(4,  6)  +  minutes_between 
End  If 
Wend 
End  Sub 


Sub  getalldata() 

'Create  coluns  on  the  "Output  Spreadsheet"  with  the  Total  Calculated  Outcomes' 

On  Error  GoTo  100 
Dim  count  As  Integer 
Dim  rownum  As  Integer 

count  =  2 
rownum  =  2 

While  count  <  50000 

While  Left(Sheets("Output").Cells(count,  1),  4)  <>  "Path" 

count  =  count  +  1 

Wend 

Sheets("output").Cells(rownum,  10)  =  LTrim(Right(Sheets("output").Cells(count,  1),  5)) 
While  Sheets("output").Cells(count,  1)  <>  "" 
count  =  count  +  1 
Wend 

Sheets("output").Cells(rownum,  11)  =  Sheets("output").Cells(count  +1,2) 
Sheets("output").Cells(rownum,  12)  =  Sheets("output").Cells(count  +  1,  3) 
Sheets("output").Cells(rownum,  13)  =  Sheets("output").Cells(count  +1,4) 
Sheets("output").Cells(rownum,  14)  =  Sheets("output").Cells(count  +1,5) 
rownum  =  rownum  +  1 
count  =  count  +  1 
Wend 
100 

End  Sub 
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Appendix  C:  Blue  Dart 


Fabio  Ayres  Cardoso,  Captain,  Brazilian  Air  Force 
(55-21)  2453-0604,  cardosoaf@yahoo.com.br 
word  count: 

Stealthy  River  Navigation  in  Jungle  Combat  Conditions 

It’s  common  sense  that  logistics  is  fundamental  for  conducting  war  at  the 
operational  level.  To  keep  fighting  in  the  inherently  chaotic  and  frictional  environment, 
combat  units  need  to  be  constantly  supported  and  resupplied  of  items  that  are  consumed 
during  the  battle.  In  certain  places,  like  the  jungle  and  tropical  rainforests,  terrain 
conditions  and  climate  offer  an  extra  dimension  to  military  operations  which  needs  to  be 
carefully  considered  in  addition  to  enemy  actions.  However,  if  local  features  are 
exploited  they  can  become  a  notable  ally  to  sustain  and  prolong  military  campaigns. 

For  example,  under  a  logistic  perspective,  jungle  surroundings  provide  excellent 
camouflage  for  stationary  entities  like  storage  depots,  airfields,  and  camps.  The  jungle 
camouflage  also  provides  a  means  to  effectively  hide  and  cover  mobile  units  such  as 
trucks,  airplanes,  and  boats.  These  deny  and  deceive  strategies  are  especially  important  to 
protect  the  supply  chains  set  up  between  depots  and  combat  forces  which  opposing  forces 
are  constantly  trying  to  disrupt.  A  key  aspect  of  military  logistics  is  to  preserve  lines  of 
communication  and  supply  between  suppliers  and  deployed  troops  under  enemy  threat. 

The  above  mentioned  covert  conduct  of  logistics,  also  called  stealthy  navigation, 
uses  local  attributes  to  conceal  movement.  In  jungle  areas,  the  most  obvious  method  to 
avoid  enemy  detection  is  to  take  advantage  of  shadows  or  dark  areas  provided  by  large 
and  dense  vegetation.  However,  it  requires  both  good  knowledge  of  the  actual 
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environment  as  well  as  skills  and  judgment  to  maximize  benefits.  Nowadays,  this  kind  of 
discernment  can  be  obtained  using  computational  models  which  combine  local 
characteristics  with  optimization  tools.  This  research  combines  a  solar  position  calculator 
to  estimate  shade  areas  along  a  river  and  a  graph  algorithms  to  compute  combinations  of 
available  paths.  When  combined,  it’s  possible  to  schedule  sorties  over  a  river  network  in 
a  jungle  context  reducing  the  likelihood  of  detection  by  enemy  forces  and  the 
confiscation  or  destruction  of  supplies. 

That  approach  offers  different  and  complementary  advantages: 

a)  It  presents  a  way  to  save  military  assets  in  combat  zones  by  reducing  the  demand  for 
armed  escort  across  watercourses.  These  forces  are  now  available  for  offensive 
actions. 

b)  It  matches  the  strategy  employed  by  jungle  warriors  in  terms  of  furtive  actions,  covert 
movements,  and  denied  existence. 

c)  It  provides  a  means  to  greatly  reduce  supply  chain  costs.  By  using  the  local 
environment  to  conceal  movement,  boats  do  not  have  to  be  made  of  technologically 
advanced  materials  that  provide  very  weak  radar  returns. 

Key  words:  jungle  warfare,  travel  risk  minimization,  river  network,  stealthy  navigation. 

“Captain  Fabio  Ayres  Cardoso  is  Master’s  degree  student  from  the  Air  Force  Institute  of 
Technology.” 
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Appendix  D:  Quad  Chart 


Stealthy  River  Navigation  in 
Jungle  Combat  Conditions 


Plotted  Network  Diagram  used  in  Model  Development 
(Satellite  photo  of  Amazon  River  basin) 


Capt.  Fabio  Ayres  Cardoso 

AFIT/ENS 

ADVISOR 
Dr.Alan  Johnson 

READER 

Major  Daniel  Mattioda.  PhD 
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unmanned  aeml  vehcle s)  Plannee  mist  select  among  the  abatable  waeicoirses  ones 
that  provide  stealthy  navigation  to  the  oomba  forte  locaion  where  the  reqtested  simples 
are  safety  unleaded 


LIMITATIONS  AND  ASSUMPTIONS 


b)  The  dynamic  nattire  of  the  waecouise  navi  ark  isn’t  modeled  The  envronmeital  effect 
of  fog  at  sunrise  and  sunsa  due  to  the  high  tumitty  of  topical  ram  forests,  as  well  as 
the  etcessuecondensaron  on  diy  seasons  e  net  modeled  Addtonaly.  the  major  river 
attributes  (Row  speed  and  depth)  as  w el  as  the  forest  rvenne  feature  varatians  throu^i 
the  path  are  also  comptetfty  neglect  el 

c)  The  nsk  cafculann  whch  cons ders  an  amnpresert  enemy  exacerbaes  foe  patertia 
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Boat  width  =  1.64  Tea  Boa  speed  =  1.86  nWhr 
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n  Survivability  (probability  of  mission  success) 
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a)  Research  identifies  benefits  of  concealment  schemes  ised  to 
avoid  detection  and  movement  tracking  in  terms  of  postion  and 
speed 

b)  Changes  in  risk  quantifier  parameters  don't  affect  proportionally 
and  equally  the  boa  vulnerability  through  the  route  Increasing  boa 
speed  demonstrates  low  impact  in  reducing  boa  vulnerability.  On 
the  other  hand,  small  vanaons  on  boat  width  produces  large 
consequences  for  risk  travel.  Thus,  boa  width  is  an  important 
consideration  during  the  stealthy  supply  watercraft  design  process 

c)  Different  analytic  approaches  to  the  same  operational  scenario 
provides  military  planners  the  abSty  to  attain  simultaneously 
tactical  and  logistics  objectives  The  outcomes  ate  compaibiy 
displayed  with  their  goals 


Boa  width  =  €.56  feet 


□ 


a)  Increase  the  number  of  units  deployed  on  the  battlefield  to  evahiae 
the  optima  start-time  and  boa  dimensions 

b)  Change  environment  feaures  such  as  river  flow  speed  (dvection  and 
mtensity),  and  river  depth  (include  seasonal  alterations}  which  can 
affect  boa  transit  capabifty. 

c)  Evaluate  the  optima  depot  location  to  mmrmee  total  travel  risk 
between  the  depot  and  the  deployed  fighting  forces  given  a  network 
and  theater  operations  boundaries 

d)  Include  the  moon  lum mosey  influence  on  exposure  leva  calculators 
during  network  night  crossing 

e)  Develop  ether  comba  logistics  effectiveness  measures  for  assessing 
trava  network  vulnerability  to  include  fua  consumption  rate,  defense 
posture,  and  detection  avoidance  from  waer  and  aerial  searches 


Department  of  Operational  Sciences 
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